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Resumo 



A produgao de pioes em colisoes nucleao-nucleao, junto ao limiar, tern sido urn desafio 
nas ultimas decadas. A reacgao pp — *> pptr em particular, e muito sensivel a mecanismos 
de curto alcance porque a conservagao do isospin suprime o termo de troca de pi5es 
que de outra forma seria dominante. Assim, tern sido muito dificil estabelecer a relativa 
importancia dos varios processos de reacgao. 

Apos rever o estado-da-arte da teoria, abordamos a validade da distorted-wave Born 
approximation (DWBA), atraves da sua ligagao aos diagramas da teoria das perturbagoes 
ordenadas no tempo (TOPT). Analisamos igualmente as escolhas possiveis para a energia 
do piao trocado, inevitaveis no formalismo nao-relativista subjacente a DWBA. 

O operador de re-dispersao resultante de TOPT e comparado com o obtido pela abor- 
dagem mais simples de matriz-S, que tern sido usada eficazmente abaixo do limiar. A 
tecnica de matriz-S, tendo reproduzido os resultados de TOPT para a re-dispersao em 
pp — > ppn° e aplicada a produgao de pioes neutros e carregados, descrevendo-se com 
sucesso as secgoes eficazes nestes diferentes canais. Os principals mecanismos de produgao 
e ondas parciais correspondentes a momento angular mais elevado sao incluidos. Final- 
mente, discutimos o efeito na secgao eficaz das aproximagoes usuais para a energia do 
piao trocado. 



PALAVRAS-CHAVE: 

Produgao de pi5es, prescrigoes de energia, 
ordenadas no tempo, re-dispersao. 



matriz-S, DWBA, teoria das perturbagoes 



Abstract 



Understanding pion production in nucleon-nucleon collisions near threshold has been a 
challenge for the last decades. In particular, the reaction pp — > pp-ir is highly sensitive to 
short-range mechanisms, because isospin conservation suppresses the otherwise dominant 
pion exchange term. However, the relative importance of the various reaction processes 
has been very difficult to establish. 

After reviewing the state-of-the-art of the theoretical approaches, we address the va- 
lidity of the distorted-wave Born approximation (DWBA) through its link to the time- 
ordered perturbation theory (TOPT) diagrams. As the energy of the exchanged pion is 
not determined unambiguously within the non-relativistic formalism underlying DWBA, 
we analyse several options to determine which one is closer to TOPT. 

The S-matrix technique, successfully used below threshold, is shown to reproduce 
the results of TOPT for the re-scattering mechanism in 7r° production. It is afterwards 
applied to full calculations of both charged and neutral pion production reactions, the 
cross sections of which are described successfully. The main production mechanisms and 
partial waves corresponding to high angular momentum are included in the calculations. 
Finally we discuss the effect on the cross section of the frequent prescriptions for the 
energy of the exchanged pion. 



KEYWORDS: 

Pion production, energy-prescriptions, S-matrix, DWBA, time-ordered perturbation the- 
ory, re-scattering. 
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Preface 



The study of pion production processes close to threshold was originally initiated to 
explore the application of fundamental symmetries to near-threshold phenomena. These 
investigations on meson production reactions from hadron-hadron scattering began in 
the fifties, when high energy beams of protons became available. A strong interdepen- 
dence between developments in accelerator physics, detector performance and theoretical 
understanding led to an unique vivid field of physics. 

Triggered by the unprecedented high precision data for proton-proton induced reac- 
tions (in new cooler rings), the interest on pion production studies was revitalised in the 
last decade. The (large) deviations from the predictions of one-meson exchange models 
controlled by the available phase-space, are indications of new and exciting physics. 

The reaction NN — > NNn is the basic inelastic process related to the nucleon-nucleon 
(NN) interaction. It sheds light on the NN and irN interactions and is key to under- 
standing pion production in more complex systems. Close to the threshold the process is 
simpler because it is characterised only by a small number of combinations of initial and 
exit channels. Moreover, at these reduced energies, meson production occurs with large 
momentum transfers, making it a powerful tool to study short range phenomena. 

Pion production occurs when the mutual interaction between the two nucleons causes a 
real pion to be emitted. The other contribution comes from a virtual pion being produced 
by one nucleon and knocked on to its mass shell by an interaction with the second nucleon. 
This is the so-called re-scattering diagram, which is found to be highly sensitive to the 
details of the calculations, namely the treatment of the exchanged pion energy. 

This pion re-scattering mechanism is suppressed in pp — > ppir due to isospin conserva- 
tion. The transition amplitude then results from a delicate interference between various 
additional contributions of shorter range. A treatment of these mechanisms consistent 
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with the NN interactions employed in the distortion of the initial and final state is essen- 
tial to clarify this large model dependence. In this work, a consistent description of not 
only neutral, but also charged pion production, is shown to be possible. 
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In this thesis we will address the problem of charged and neutral pion production in 
nucleon-nucleon collisions. The main steps of this investigations are: 

1. Starting with relativistic field theory, time-ordered perturbation theory (TOPT) is 
used to justify the distorted- wave Born approximation (DWBA) approach for pion 
production; 

2. Using the DWBA amplitude fixed by TOPT as a reference result, the effect of the 
traditional prescriptions for the energy of the exchanged pion in the re-scattering 
operator is analysed; 

3. Defining a single effective production operator within the S-matrix technique, its 
relation to the DWBA amplitude yielded by TOPT is established; 

4. Employing the S-matrix approach, charged and neutral pion production reactions 
are described consistently 

This work is based on the following publications: 

• V. Malafaia and M. T. Pena, Pion re-scattering in ir° production, Phys. Rev. C 69 
(2004) 024001 [nucl-th/0312017]. 

• V. Malafaia, J. Adam and M. T. Pena, Pion re-scattering operator in the S-matrix 
approach, Phys. Rev. C 71 (2005) 034002 [nucl-th/0411049]. 

• V. Malafaia, M. T. Pena, Ch. Elster and J. Adam, Charged and neutral pion pro- 
duction in the S-matrix approach, [nucl-th/0511038], submitted to Phys. Lett. B. 

• V. Malafaia, M. T. Pena, Ch. Elster and J. Adam, Neutral and charged pion pro- 
duction with realistic NN interactions, to submit to Phys. Rev. C. 

Chapter ^ introduces the physics accessible through the study of the meson production 
reactions. It focuses on the specific aspects of hadronic meson production reactions close 
to threshold, namely the rapidly varying phase-space, the (high) initial and (low) final 
relative momenta, and the general energy dependence of the production operator. 

Chapter |21 is a historical review of the main theoretical approaches to pion production 
developed so far. The first part is dedicated to the distorted-wave Born approximation 
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(DWBA), the frequently used approach in all the calculations near the threshold energy, 
in which the nucleon-nucleon interaction is treated non-perturbatively, whereas the tran- 
sition amplitude iViV — > NNtt is treated perturbatively. The second part focuses on 
the coupled-channel phenomenological approaches and the third part aims to present the 
actual status of chiral perturbation theory (xPT) calculations. 

DWBA calculations apply a three-dimensional formulation for the initial- and final- 
NN distortion, which is not obtained from the Feynman diagrams. As a consequence, 
in calculations performed so far within the DWBA approximation, the energy of the 
exchanged pion has been treated approximately and under different prescriptions. A 
clarification of these formal issues is thus needed before one can draw conclusions about 
the physics of the pion production processes. The re-scattering mechanism, being highly 
sensitive to these energy prescriptions, is the ideal starting point for this clarification. 

Chapter El aims to obtain a three-dimensional formulation from the general Feynman 
diagrams. This chapter discusses the validity of the DWBA approach by linking it to the 
time-ordered perturbation theory (TOPT) diagrams which result from the decomposition 
of the corresponding Feynman diagram. Since in the time-ordered diagrams energy is not 
conserved at individual vertices, each of the re-scattering diagrams for the initial- and 
final-state distortion defines a different off-energy shell extension of the pion re-scattering 
amplitude. This imposes the evaluation of the matrix elements between quantum me- 
chanical wave functions, of two different operators. Henceforth, in Chapter 0] we are lead 
to an alternative approach to the TOPT formalism (the S-matrix approach), which in 
contrast to it avoids this problem. 

The S-matrix approach provides a consistent theoretical framework for the two di- 
agrams, as well as for the NN distortion. Besides, from the practical point of view, 
it simplifies tremendously the numerical effort demanded in TOPT by the presence of 
logarithmic singularities in the pion propagator for ISI. The S-matrix construction has 
already been successfully used below pion production threshold and in particular for NN 
interactions and electroweak meson exchange currents. In Chapter 0] also, the effective 
DWBA amplitude obtained in Chapter El is employed to re-examine the nature and ex- 
tent of the uncertainty resulting from the approximations made in the evaluation of the 
effective operators. 

4 



Preface 



In Chapter the S-matrix technique, shown in Chapter 0] to reproduce the results 
of time-ordered perturbation theory for the re-scattering mechanism in ir° production, is 
applied to charged and neutral pion production reactions. The major production mech- 
anisms, namely the contributions from the direct-production, re-scattering, Z-diagrams 
and A-isobar excitation are considered. Higher angular momentum partial waves, which 
are not included in traditional calculations, are also considered. The last part is dedicated 
to the effect on the cross section of the usual prescriptions for the energy of the exchanged 
pion. For all the charge channels, the S-matrix approach for the description of the pion 
production operators reproduces well the DWBA result coming from TOPT. Importantly, 
the effect of some approximations usually employed is also assessed. The ir + reaction is 
seen to be especially sensitive to those. Previous failures in its description are overcome 
and clarified. 

Finally, in Chapter El we outline and summarise the most relevant aspects discussed 
in this thesis, and mention the future prospects of pion production in nucleon-nucleon 
collisions. 
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Meson production close to threshold 
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Abstract: Meson production reactions in nucleonic collisions near threshold are a pow- 
erful tool to investigate short-range phenomena. Pion production reactions play a very 
special role, since they yield the lowest hadronic inelasticity for the nucleon-nucleon inter- 
action and thus they are an important test of the phenomenology of the nucleon-nucleon 
interaction at intermediate energies. 
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Meson production close to threshold 



1.1 What physics can we learn from meson produc- 
tion reactions? 

Quantum Chromodynamics (QCD), the fundamental theory for strong interactions, 
unfolds an impressive predictive power mainly at high energies. However, at low energies 
the perturbative expansion no longer converges. 

Although a large amount of data on hadronic structure and dynamics is available 
from measurements with electromagnetic probes (for instance, from MAMI at Mainz, 
ELSA at Bonn and JLAB at Newport News), there is still much to be learned about 
the physics with hadronic probes at intermediate energies, comprising the investigation of 
production, decay and interaction of hadrons. In particular, meson production reactions 
(close to threshold) in nucleon-nucleon, nucleon-nucleus and nucleus-nucleus collisions 
constitute a very important class of experiments in this field. 

With the advent of strong focusing synchrotrons having high-quality beams (for in- 
stance, the IUCF Cooler in Bloomington, CELSIUS in Upsala and COSY in Jiilich), a 
new class of experiments could be performed in the last decade, differing from the previous 
ones with respect to the unprecedented quality of the data (polarised as well as unpo- 
larised) for several NN — > NNx reaction channels (a recent review on the experimental 
and theoretical aspects of meson production can be found in Refs. [H Ej and in Ref . jB] , 
respectively) . 

The study of meson production close to threshold has several attractive key features, 
in particular concerning, 

(i) Large momentum transfer in the entrance channels 

Meson production near threshold occurs at large momentum transfers and therefore 
is a powerful tool to study short range phenomena in the entrance channel; 

(ii) Simplicity of the entrance and exit channels 

The analysis of the reaction data is straightforward allowing one to study the un- 
derlying reaction mechanisms; 

(iii) Small phase-space 

Only a very limited part of the phase-space is available for the reaction products 
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and hence, only a few partial waves contribute to the observables, especially very 
near the threshold. 

Although a small number of partial waves is clearly an advantage for calculations, on the 
other hand, as a result of the small phase-space, the cross sections are also small. There 
is a delicate interplay between the several possible reaction mechanisms and thus it is 
essential that all the technical aspects are under control for a meaningful interpretation 
of the data. 

The large momentum transfer can also look like an disadvantage since it is difficult 
to reliably construct the production operator. However, in the regime of small invariant 
masses, the production operator is largely independent of the relative energy of a partic- 
ular particle pair in the final state 1 . Consequently, dispersion relations can be used to 
extract low-energy scattering parameters of the final state interaction and at the same 
time, to estimate the error in a model independent way. 

In the overall, meson production reactions in nucleonic collisions have a huge poten- 
tial to give insight into the strong interaction physics at intermediate energies, namely 
concerning the following aspects: 

• Final state interactions 

Scattering of unstable particles off nucleons is experimentally very problematic since 
it is difficult to prepare intense beams of these particles with the required accuracy. 
Production reactions where such particles emerge in the final state are an attractive 
alternative. 

• Baryon resonances in a nuclear environment 

The systems studied in nucleon-nucleon and nucleon-nucleus collisions allow to in- 
vestigate particular resonances in the presence of other baryons and excited by 
various exchanged particles. One example is the N* (1535), which is clearly visible 
as a bump in any r) production cross sectionjH]. 

• Charge symmetry breaking 

The existence of available several possible initial isospin states (for instance pp, pn 
: If there are resonances near by, this statement is no longer true. 
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and dd), with different possible spin states, allows experiments which enable the 
study of isospin symmetry breaking. 

• Effective field theory in large momentum transfer reactions 

The investigation of hadronic processes at low and intermediate momenta are es- 
sential to test the convergence of the low-energy expansion of chiral perturbation 
theory (%PT). 

• Three-nucleon forces 

The information on the short-range mechanisms that can be deduced from pion 
production in nucleon-nucleon (NN) collisions is relevant for constraining the three- 
nucleon forces jl] . 

This work will focus on pion production reactions in nucleon-nucleon collisions. It is 
the lowest hadronic inelasticity for the nucleon-nucleon interaction and thus an impor- 
tant test of the understanding of the phenomenology of the iViV interaction. Secondly, 
as pions are the Goldstone bosons of chiral symmetry, it is possible to study pion pro- 
duction also using xPT. This provides the opportunity to improve the phenomenological 
approaches via matching to the chiral expansion, as well as to constrain the chiral contact 
terms via resonance saturation. As mentioned before, a large number of (un)polarised 
data is available 2 to be used as constraints. Moreover, meson-exchange models of the 
nucleon-nucleon interaction above the pion threshold rely on detailed information about 
the strongly coupled inelastic channels which must be treated together with the elastic 
interaction. Also, information on pion production in the NN system is required in models 
of pion production or absorption in nuclei. 

2 From the experimental point of view, it is important to notice that since the vector mesons have 
much larger widths compared with the pscudoscalar mesons, their detection is very difficult on top of 
a large physical background of multi-pion production events. Secondly, as there seems to be a general 
trend that the larger the mass, the smaller the cross section, the generally heavier vector mesons have 
smaller production cross sections, are thus much harder to investigate than the lighter onesPJ. 
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1.2 Specific aspects of hadronic meson production 
close to threshold 

1.2.1 Rapidly varying phase-space 

Threshold production reactions are characterised by excess energies which are small 
compared to the produced masses. In the near threshold regime the available phase- 
space changes very quickly (although remaining small). Therefore, to compare different 
reactions, one needs an appropriate measure of the energy relative to threshold. For pion 
production, the traditionally used variable is the maximum pion momentum t] (in units 
of the pion mass). For all heavier mesons, the so-called excess energy Q, defined as 

Q = y/s - V 's threshold , (1.1) 

is used instead. In Appendix |XJ a compilation of useful kinematic relations is presented, 
and the importance of relativistic kinematics for the near threshold reactions is stressed. 

If Q gives the available energy for the final state, the interpretation of rj is some- 
what more involved. In a non-relativistic, semiclassical picture, the maximum angular 
momentum allowed can be estimated via 

Imax ^ Rq', (1.2) 

where R is a measure of the force range and q' is the typical momentum of the corre- 
sponding particle. Identifying R with the Compton wavelength of the meson of mass m x , 
T) can be interpreted as the maximum angular momentum allowed [5]: 

q< 

Imax- -V (with ft = C = 1). (1.3) 

m x 

To compare the cross sections for reactions with different final states in order to extract 
information about the reaction mechanisms, one has to choose carefully the variable that 
is used to represent the energy. Indeed, as it is shown in Fig. Il.l|. the total cross sections 
for pp — > ppir , pp — > pprj and pp — > pprj' are different when compared at equal rj (panel 
(b)) or at equal Q (panel (a)). When the dominant final state interaction is the pp 
interaction, which is the case for those reactions, it appears thus to be more appropriate 
to compare the cross sections at equal Q, since then at any given energy, the impact of 
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Figure 1.1: Total cross sections for the reactions pp — > ppir , pp — » pprj and pp — > pprj' as 
a function of a)the excess energy Q and b)the maximum pion momentum rj in units of the 
mass of the produced meson. Figure taken from Ref. 



the final state interaction is equal for all the reactions. This is not the case for equal 
values of rj, as rj depends on the mass of the produced meson. 

1.2.2 Initial and final relative momenta 

Meson production in nucleon-nucleon collisions requires that the kinetic energy of the 
initial particles is sufficiently high to put the outgoing meson on its mass shell. In other 
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words, the relative momentum of the initial nucleons must exceed the threshold value, 



m 



2 



&iai=]/rn x M+^, (1.4) 

where M is the nucleoli mass and m x is the mass of the produced meson. For a close- 
to-threshold regime, the particles in the final state have small momenta and thus plnluai 
also sets the scale for the typical momentum transfer. In a non-relativistic picture, this 
large momentum transfer translates into a small reaction volume, characterised by a size 
parameter 3 , 

R ~ -rr — 0.5fm for pion production. (1.5) 

Pinitial 

The two nucleons in the initial state have thus to approach each other very closely before 
the production of a meson can happen. For this reason it is important to understand not 
only the elastic but also the inelastic NN interaction to obtain quantitative predictions. 



1.2.3 Remarks on the production operator and nucleonic dis- 
tortions 

In the near threshold regime, all the particles in the final state have low relative 
momenta and thus can potentially undergo strong final state interactions (FSI) that can 
induce strong energy dependencies. On the other hand, close to the threshold, the initial 
energy is significantly larger than the excess energy Q, and consequently the initial state 
interaction (ISI) should at most mildly influence the energy dependence. The dependence 
of the production operator on the excess energy should also be weak, since it is controlled 
by the typical momentum transfer, which is significantly larger than the typical outgoing 
momenta. Fig. II .21 illustrates the energy dependence of the production operator both FSI 
and ISI cases. 

The large momentum transfer characteristic of meson production reactions at thresh- 
old leads to a large momentum mismatch for any one-body operator that might contribute 
to the production reaction. 

3 As in Eq. I|1.3fl . h = c = 1 is assumed. 
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Figure 1.2: Sketch of the excess energy Q-dependence of the meson production reaction. A 
is the production operator. The left diagram shows the complexity possible, with interaction 
between the three final particles, whereas the right diagram shows the first and potentially 
dominant term (with FSI only for the two final nucleons). 



1.3 Theoretical considerations on AW — ► NNx reac- 
tions 

Most of the theoretical models 4 for the reactions NN —>■ NNx can be grouped in two 
distinct classes: 

• Distorted wave Born approximation (DWBA) 

A production operator is constructed within some perturbative scheme approxima- 
tion and then is convoluted with the nucleon wave functions. 

• (Truly) Non-perturbative approaches 

Integral equations are solved for the full (NN, NNx) coupled- channel problem, de- 
scribing multiple re-scattering and preserving three-body unitarity. 

The great majority of theoretical studies of pion production in nucleon-nucleon collisions 

in the threshold region have been done within the DWBA formalism jOj. This approach is 

motivated by the fact that close to pion production threshold, where the kinetic energy 

of the particles in the final state is practically zero, the forces between the nucleons are 

4 The development of theoretical approaches for the reactions NN — > NNx has a long history. A 
review of earlier works can be found in Rcfs. [7|[H]. 
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a) b) c) 

Figure 1.3: Contributions to the production operator in pion production in iViV collisions 
which in the DWBA formalism are convoluted with NN wave functions: a) direct-production, 
b) re-scattering and c) short-range contributions. The solid and dotted lines represent the 
nucleons and pions, respectively. 

much stronger than the interaction between the pion and the nucleon. Consequently, only 
the interaction between the nucleons is taken into account up to all orders, for instance, 
by employing wave functions that are solutions of a scattering (Lippmann-Schwinger) 
equation, whereas the pion production process is treated perturbatively and the pion is 
assumed to propagate freely after its production. Typically, the diagrams that contribute 
are of the type of those of Fig. 11.31 

On the other hand, all calculations performed so far for pion production within the full 
(NN, NNx) coupled-channel approach were within the framework of time-ordered per- 
turbation theory [TP] (TOPT), or its extension to the iVA sector done by the Helsinki [TT| 
H2HE], Argonne [H | ITH | \TR \ ITT] and Hannover [TH ] ITH ] l2*0j groups, which have a reasonable 
predictive power at higher energies but cannot describe the physics very near threshold, 
as discussed in Sec. 12.2.11 

1.4 The cross section 

In Fig. 11.41 there is an overview of the total cross sections in pp interactions below 
4GeV beam momentum. In Table 11.11 we list the threshold momenta p thr ' and threshold 
laboratory energies , for the AW — > NNn reactions considered in this work. The 
right column is a compilation of references with the experimental determination of a to t- 

The first threshold which opens with increasing beam energy is 7r° production followed 
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Figure 1.4: Cross sections for pp interactions as a function of the beam momentum. Figure 
taken from Ref. 

very soon after by 7r =l= production (see Table fT"Tj) . Pion production exhausts all inelasticity 
in this momentum range (see Fig. I1.4JI and therefore is fundamental to understand the 
nucleon-nucleon interaction. 
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reaction p thr - [MeV] T^- [MeV] 



pp - 


-* ppn 


724.4 


279.7 


[211 E21 E3 EH ESI E 


pp - 


-> pnir + 


737.3 


289.5 


|2Zl EH1 EHl ED] 


pn - 


-> pp7r~ 


737.3 


289.5 


[2HE21E21 



Table 1.1: Threshold momenta and threshold laboratory energies for pion production reactions 
in the NN collisions considered in this work. The right column refers to the experimental 
determination of the corresponding cross sections. 
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a) 



b) 



c) 



d) 



Figure 2.1: The main mechanisms considered in the works on pion production: a) direct 
production, b) re-scattering, c) pion re-scattering via a A and d) short-range contributions. 
The solid and dashed lines represent the nucleon and pion fields, respectively. Heavy mesons 
(cr, lo,...) are represented by the wavy line. 



2.1 DWBA in the Meson exchange approach 

2.1.1 Problems of the earlier calculations 
The work of Koltun and Reitan 

Pioneering work on pion production was done in the 1960s by WoodrufjHj and by 
Koltun and ReitanJS^. All later investigations of meson production, including the very 
recent efforts to analyse the high-precision data from the new generation of accelerators, 
have followed basically the same approach, if one excludes the Hannover model [TH | ITT H l2*U] 
for coupled NN, NA-NNir channels. 

These works focused on the reactions pp — > ppir and pp — > dir + . The processes 
considered were direct production by either nucleon (diagram (a) of Fig. 123} , the so-called 
impulse approximation, and production from pion-nucleon scattering (diagram (b)), the 
so-called re-scattering term. The nN — > ttN transition amplitudes were parameterised in 
terms of scattering lengths, through the Hamiltonian 1 7i = 7i\ + TI2, with 



rii ■ fnNN -, 

Hi = % cr 



V.T • 7T + — (pT ■ 7T + T ■ -hp) 



{2.1) 



H2 = 47T— !-7T ■ 7T + 47T— ^-T • 7T X 7T (2.2) 

where a and r are the usual nucleon spin and isospin operators, and p is the nucleon 



1 Actually, Eq. 112. Ill and Eq. H2.2J) are the Lagrangian, but we kept here the term Hamiltonian for 
historical reasons. 
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momentum operator. The pion(nucleon) mass is (M) and the pion field is 7r. The 
ttNN pseudovector coupling constant is f n NN- 

Tii is obtained from a non relativistic reduction of the pseudovector NNw vertex 2 and 
gives diagram (a) in Fig. 12.11 The first term of Eq. (|2.1|) represents p-wave irN coupling, 
while the second term ( "galilean" term) accounts for the nucleon recoil effect [SB] • For s- 
wave pion production, only the second term contributes. Since this second term is smaller 
than the first term by a factor of ~ the contribution of the Born term to s-wave pion 
production is intrinsically suppressed, and as a consequence the process becomes sensitive 
to two-body contributions, Fig. 12.1T b) and (d). 

7^2 is a phenomenological effective Hamiltonian describing pion re-scattering. The 
isoscalar and isovector parts of the nN s-wave scattering amplitude, Ai and A2, were 
obtained from the s-wave phase shifts S\ and £3 for the pion-nucleon scattering, through 
the Born-approximation relations |34j. 

Ax = ~ (5 1 + 25 3 ) = 0.005 A 2 = ~ (5 1 - 5 3 ) = 0.045, (2.3) 

or] or] 

where r] is the maximum pion momentum (in units of m n ). Although A2 ^> Ai, the isospin 
structure of the A2 term is such that it cannot contribute to it production. 

We note that without initial- or final-state distortions, the diagram (a) of Fig. 12.11 
vanishes because of four-momentum conservation. The calculations of Ref. |B~3] were 
performed using the Hamada- Johnston phenomenological potential for the NN distortion, 
and neglecting the Coulomb interaction between the two protons. The calculated cross 
section for pp — > ppn of YJrfiib was found to be consistent with the measured cross 
section near threshold, which was, however, not well determined |37j by that time. 

The role of final-state and Coulomb interactions 

When the first high precision data [21] on the reaction pp — > ppn appeared, they con- 
tradicted the predicted rf dependence near threshold. According to the work of Refs. [2H 
13*%] the energy dependence of the s-wave cross section followed from the phase space and 
2 In the chiral limit for vanishing momenta, the interaction of pions with nucleons has to vanish. Thus 
the coupling of pions naturally occurs either as derivative or as an even power of the pion mass. In general, 
the pseudovector coupling for the irNN vertex is preferred. The pseudovector coupling automatically 
incorporates a strong(weak) attractive p-wave(s- wave) interaction between pions and nucleons|35j. 

22 



2.1 DWBA in the Meson exchange approach 



23 



so 



40 



* 30 



3 20 



1 1 F*-l 1 |—T 1 1 1 1 1 1 TT-| 1 1 I 



pp — ppir 



without 

Coulomb 



10 7 /& 



with 
Coulomb 




L_i i i i_l ■ i i i I i i i i 1 i ■ i 

O 0.2 0.4 0.6 



10 



Figure 2.2: Energy dependence of the total cross section for the reaction pp — > ppir . The 
solid dots are data from Ref. (25 . The squares, cross and diamonds are older data. The solid 
curve is the full calculation and the dashed line shows the effect of omitting the Coulomb 
interaction. The data are to be evaluated using the scale on the left, while the theory used 
the scale on the right. Note that the calculations underestimate the data by a factor of ~ 5. 
Figure taken from Ref. |38J. 



a simple treatment of the final-state interaction |39[ 14*0] between the two (charged) pro- 
tons. This was found sufficient to reproduce the shape of the measured cross section up 
to 7] ~ 0.5 (see Fig. 12. 2|) . where higher partial waves start to contribute [22]. Also, the 
inclusion of the Coulomb interaction was found to be essential to describe the energy 
dependence of the total cross section in particular for energies close to threshold. The 
validity of the effective-range approximation employed in Ref. for the energy depen- 
dence of the final state turned out to be limited to energies rather close to threshold 
(V < 0-4). 

However, the theory failed in describing quantitatively the cross section for pp — > ppir 
by a factor of 5, which was in contrast to the reaction pp —* pmr + , where the discrepancy 
was less than a factor of 2, as reported in Ref. [22] . Ref. jHH] suggested that the problem 
arose from the use of an over-simplified pion nucleon interaction, namely in considering the 
exchanged pion to be on-shell. The on-shell s-wave pion nucleon interaction is constrained 
to be small by the requirements of chiral symmetry but, for the production reaction to 
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Figure 2.3: Cross section for pp — > ppir . The solid line is the full calculation. The dotted 
line is the direct production due only to second term of Eq. (12.1)) and going beyond the static 
approximation, the dash-dotted line includes also re-scattering through the A. The dashed 
line is the full purely nucleonic production. All the calculations were multiplied by a common 
arbitrary factor of 3.6. The data points are from Ref. [21] ■ Figure taken from Ref. [41] . 



proceed, either the re-scattered pion or a nucleon must be off-shell. This means that the 
ttN amplitude relevant for pion production must be larger than the theoretical one for 
on-shell particles, whose investigation was to be subsequently pursued. 

The inclusion of the A 

The work of Ref. jH] considered the p-wave re-scattering through a A (1232) resonance 
(diagram c) of Fig. 12.1)1 by introducing finite-range NA coupled-channel admixtures to 
the nucleonic wave functions. The transition potential NN — > NA included both it and 
p exchanges. The pion production vertex was taken from Eq. ()2.2)) and including the 
relativistic effects arising from the use of the pion total energy E n instead of the pion 
mass (static approximation). The isoscalar and isovector parameters Ai and A2 of the 
phenomenological hamiltonian of Eq. (|2.2|) were allowed an energy dependence through 
the momentum of the pion. The results are shown in Fig. 12.31 

The small difference between the theoretical models was attributed to the different 
value used for the nN coupling constant and to the relativistic kinematics. The inclusion 
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of the non-galilean term and the s-wave scattering gave an enhancement of over 60% 
(dashed curve). Another 25% enhancement arose from the inclusion of the re-scattering 
through the A. However, the cross section was still missed by almost a factor of ~ 4. 



Charged pion production 

The first calculations on n + production where those of Schillaci, Silbar and Young[42j. 
General isospin and phase space arguments HO] were employed to predict the spin, 
isospin and energy dependence of the total cross section, including all partial waves for 
the NN amplitudes, but accounting only for s-wave pion-nucleon states. It fails if con- 
tributions from the A resonance are significant. 

The second prediction was made by Lee and Matsuvama|16j with a coupled channel 
formalism that focused on the effects of the A intermediate state. In Ref. the A 

process is handled rigorously while the non resonant pion production process is introduced 
as a perturbation. The estimated A contribution was roughly 15% of the total cross 
section. Both calculations were not able to describe the experimental data for pp — > pnn + 
as it is shown in Fig. 12.41 (solid and dotted lines, respectively). 

Calculations based on a relativistically covariant one-boson exchange model [13], also 
suggested that the contribution of a A-isobar is not important at energies below 350MeV 
(due to the fact that at lower energies pions are predominantly produced in a ttN relative 
s-state and thus the possibility of forming a A-isobar is greatly reduced), but dominate 
at higher beam energies. The cross sections for pp — > pnir + near threshold were under- 
predicted a factor of 2-4. 

The work of Ref. [H] applied the Watson theoremjSH] 3 to the final state interaction 

3 In 1952 Watson showed that for a short-range strong (attractive) NN interaction and in the regime 
of low relative energies of the interacting particles, the energy dependence of the total NN — > NNx cross 
section is determined only by the phase space and by the on-shell NN T-matrix, 

Vnn^nnx (v) / dp(q)\T (k 0} k )\ 2 <x / dp(q) 
Jo JO 

Here, the momentum of the outgoing meson is q and dp (q) denotes the phase space. The relative 
momentum on the final nucleons is fco, and S (ko) are the corresponding phase-shifts of the final NN 
subsystem (restricted to s- waves). Recently, the work of Ref. [321 concluded that Watson's requirement 
of an attractive FSI is unnecessary to obtain the energy dependence of the cross section given by Eq. (|2.4|l . 
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Figure 2.4: Comparison of measured total pp — > pnix + cross section with the theoretical 
predictions of Refs. [HI E2|. The solid lines are the calculations of Schillacci, Silbar and 
Young|42j. The dotted lines are the calculations of Lee and Matsyuama[T7]. Figure taken 
from Ref. (22] 

of Ref. jlH]- The calculated cross sections were within 25% of the pp — ► pnii + data[28 , 
but were also under-predicted near threshold. 

Also, in the work of Ref. [IB], the isoscalar heavy meson exchange found to dominate 
in pp — > ppn was however shown to be less significant in pp — > pmr + , where the re- 
scattering diagram was the most important one. Further theoretical studies on these 
issues were then clearly needed to clarify the role of the different production mechanisms. 

2.1.2 The first quantitative understandings 

The first quantitative understanding of the pp — > ppir data was reported by Lee and 
Riskaj!?] and later confirmed by Horowitz et al. |48] . where it was demonstrated that 
short range mechanisms (diagram (d) of Fig. 12. lj) can give a sizeable contribution. In 
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these works, the difficulty in describing the cross section for pp — * ppn was overcame 
by considering the pair terms, positive and negative energy components of the nucleon 
spinors, connected to the isoscalar part of the NN interaction. 



The importance of short-range mechanisms 

In the work of Ref. the short range two-nucleon mechanisms that are implied by 
the nucleon-nucleon interaction were taken into account by describing the pion-nucleus 
interaction by the extension of Weinberg's effective pion-nucleon interaction to nuclei: 



C = -^-A» ■ <Vr, (2.5) 



where is the isovector axial current of the nuclear system 4 . This formulation reduces 
the calculation of matrix elements for nuclear pion production to the construction of 
the axial current operator, which is formed of single-nucleon and two-nucleon (exchange) 
current operators. The single-nucleon contribution is 5 



one— body ^ / j 



i=l,2 



. P'i+Pi (j) 

2M 



(2.6) 



where pi (pi) is the initial(final) nucleon momenta. When the nucleon-nucleon interaction 
is expressed in terms of Fermi invariants (scalar, vector, tensor and axial-vector) there is 
an unique axial exchange charge operator that corresponds to each invariant. The general 
two-body-exchange charge operator is then 

AL-bod y = 7^73 [A° (S) + A (V) + A (T) + A (A)] . (2.7) 
(2tt) 

The axial exchange charge operators associated with the scalar (S) and vector (V) com- 
ponents of the NN interaction are the most important. Dropping terms that involve 



4 The relationship between the current and the amplitude M. is Q^A 11 = f^Ai, where Q = (E^^q) 
is the four-momentum of the emitted pion. Near threshold, the interactions that involve s-pions should 
dominate, and thus the amplitude is simply given by M. = —j-A^E^, which coincides with the second 
term of Eq. (|2.1|l . 

5 The conventional single-nucleon pion production operator (first term of Tii of Eq. (|2.1|l ) is recovered 
by using the Goldbcrgcr-Trciman relation M- = . 
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isospin flip and therefore do not contribute to it production, A (S) and A (V) readjlHj 
A ° {S) = K^ rW + v s ^ r(2) ] • ^ + 0- " 2 ) > ( 2 - 8 ) 



A°(V) = |{[4Wr (1) +^(fcV 2) 



(2.9) 



+^zuy(jfe) (r (1) x r (2) ) 8® • jfe J + (1 ^ 2) . 

The momentum operators are defined as Pj = and = P2 ~ P2 = Pi — p{ is 

the momentum transfer. The momentum dependent potential functions are isospin 
independent (+) and isospin dependent (— ) functions associated with the corresponding 
Fermi invariants. These functions can be constructed from the components of complete 
phenomenological potential models [23? or alternatively, by employing phenomenological 
meson exchange models [1H]. 

The contribution of the s-wave pion re-scattering to the pp — > ppir amplitude can be 
included by adding to A° nc „ body + A° wo „ body the following two-body axial charge operator 

AOw ^7^^ i ^^i /( * 2) + (1 " 2) - (2 - 10) 

(2ir) m w m 2 _j_ ^ 

where f(k 2 ) is a monopole form factor (to be consistent with the Bonn boson exchange 
model for the nucleon-nucleon potential) and E n the energy of the produced pion. Note 
that A (it) of Eq. (j2.1()j) does not include the dependence on the energy of the exchanged 
pion (i. e., the static approximation is employed here). 

The results for the pp — >• pptr are shown in Fig. 12.51 As already found in previous 
works [HH1 141 j. the impulse and re-scattering mechanism were not enough to describe the 
data (dot-dashed line of Fig. The short-range axial charge operators enhance largely 
the cross section and remove most of the under prediction (solid and dotted lines of 
Fig. 12.5(1 . However, it was also found that for both iViV potentials employed (Paris and 
Bonn), the energy dependence of the data was not reproduced in detail. According to 
Ref. [Hj, this might be due to the neglect of the energy dependence of the parameter Ai 
in the effective amplitude of Eq. (|2.10j) . or to p-wave contributions and nNN three-body 
scattering distortions in the final state. Although these corrections were expected to be 
small, they could in principle lead to significant contributions through the interference 
with large amplitudes and thus have large significant effects on the predicted energy 
dependence. 

28 



2.1 DWBA in the Meson exchange approach 



29 




Figure 2.5: Cross section for pp — > ppn°. The so I id (dotted) line is the full calculation using 
the Bonn(Paris) potential to construct the axial exchange charge operator. The dot-dashed 
line is obtained by keeping only the one-body term of Eq. (12. bp and the pion re-scattering 
term of Eq. (12.101) using the Paris potential. The data points are from Refs. (2HI23- Figure 
taken from Ref. [47J. 

Shortly after the work of Ref. jTT] , where the meson-exchange contributions were cal- 
culated from phenomenological potentials, Ref. [IB] used an explicit one-boson exchange 
model for the NN interaction and for the calculation of the MEC's. As in Ref. the 
irN re-scattering vertex was restricted however to the on-shell matrix element and to s- 
wave pion production. The largest contribution was found to come from the Z-diagrams 
mediated by a-exchange, which was of the order of the one-body (impulse) term. The 
next important contribution was from the u meson Z-diagrams (35 — 45% of the one-body 
term contribution). 

The importance of off-shell effects 

Shortly after the discovery of the importance of the short-range mechanisms, Hernandez 
and OsetjHTJ demonstrated, using various parameterisations for the ttN — > ttN transition 
amplitude, that its strong off-shell dependence could also be sufficient to remove the 
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Figure 2.6: Cross section for — > pp7r° using the Bonn B potential for the initial and final 
state interaction between the two protons. All calculations include both impulse and pion 
re-scattering diagrams. Figure taken from Ref. [53J. 



discrepancy between the Koltun and Reitan model [33] and the data. 

The importance of the off-shell amplitudes was also seen in a relativistic one boson ex- 
change model|221- I n the work of Ref. [53], the model independent off-shell ttN amplitude 
obtained by current algebra (and used previously in the Tucson-Melbourne three-nucleon 
force) was also considered as input for the pion re-scattering contribution to pp — > ppir 
near threshold. It was found that this pion re-scattering contribution, together with the 
direct-production term, provided a good description of the ir° production data, when 
the current algebra nN amplitude parameters were updated with the phenomenological 
information obtained from the new meson factory irN scattering data (see Fig. 12. 6|) . 

Z-diagrams: perturbative vs. non-perturbative 

In the succeeding years many theoretical efforts were made for the calculation of the 
pp — > ppn cross section. In Refs. |43j l5lj covariant one boson exchange models were 
used in combination with an approximate treatment of the nucleon-nucleon interaction. 
Both models turned out to be dominated by heavy meson exchanges, thus giving further 
support to the picture proposed in Refs. [HJ EH] • 
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Figure 2.7: Total cross section for pp — > ppn as a function of the proton laboratory energy. 
The net effect of the resonances is small. Figure taken from Ref. |5Zj . 



However, in Ref. jHS] the negative energy nucleons were re-examined using the covari- 
ant spectator description, for both the production mechanism and for the initial and final 
state pp interaction. This approach differs crucially from earlier ones by including non 
perturbatively the intermediate negative-energy states of the nucleons 6 . 

The perturbative result for the direct-production diagram was found to be about 3 
times larger that the non-perturbative one. Although the calculation in Ref. [S3] did not 
include the re-scattering diagram contribution, it showed that the sensitive cross section 
for 7T° production seemed to be an ideal place to look for effects of relativistic dynamics. 

The role of the nucleon resonances 

Additional short-range contributions were also suggested, namely the p — uj meson 
exchange current [5^] . resonance contributions (SHI E3 EE] (see Fig. I2.7|) and loops that 
contain resonances |58j. All those, however, turned out to be smaller when compared to 
the heavy meson exchanges and the off-shell pion re-scattering. 

6 As mentioned before, in perturbative approaches these contributions (often called Z- diagrams) are 
simulated by the inclusion of effective meson-exchange operators acting in two-nucleon initial and final 
states. 
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2.2 Coupled-channel phenomenological calculations 

The starting point of the NN-NA approach is the recognition that the nucleon is a 
composite system. Since the A isobar is the most important mode of nucleonic excitation 
at intermediate energies, a possible process contributing (in second order) to nucleon- 
nucleon elastic scattering is the transition from a pure nucleonic state into a state into 
a nucleon plus a A (or two A's) with the inverse transition taking the system back to 
a two-nucleon state again. From this point of view, the nucleon-nucleon problem is a 
coupled-channel system involving at least the NN and the NA channels [H]- 

2.2.1 The Hannover model 

The Hannover model[T%| ITHj l2Ti] for the NN system considers the A isobar and pion 
degrees of freedom in addition to the nucleonic one. The model is based on a hamiltonian 
approach within the framework of non covariant quantum mechanics. In isospin-triplet 
partial waves, it extends the traditional approach with purely nucleonic potentials. It 
is constructed to remain valid up to 500MeV CM energy. The Hilbert space considered 
comprises NN and NA basis states, connected by transition potentials. Pion production 
and pion absorption are mediated by the A isobar excited by ir and p exchange. The model 
accounts with satisfactory accuracy for the experimental data of elastic nucleon-nucleon 
scattering, of the inelastic reactions pp — > ir + d and of elastic pion-deuteron scattering. 

Lee and Matsuyama also performed calculations |15 [ ITH j IT7j for pion production within 
a coupled- channel approach, which differed from that of the Hannover group mainly by 
the treatment of the energy in the A propagator. Both the theoretical predictions |1 61 IT§] 
from the Hannover and the Lee and Matsuyama models for the pp — > pmr + differential 
cross section were found to be quite sensitive to the inclusion of the NA potential, but 
the under-prediction of the data could not be completely removed. The calculations 
of Refs. [T6t IT9] considered an energy region well above pion production threshold (~ 
580MeV — 800MeV), since pion production is assumed to occur only via an intermediate 
A excitation and thus the details of the irN amplitude (related to chiral symmetry and 
the chiral limit), which are important close to threshold, were not included. 
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2.2.2 The Jiilich model 

The Jiilich model [23 EH attempts to treat consistently the NN and the irN 
interaction for meson production close to threshold, taking both from microscopic models. 
Although not all parameters and approximations used in the two systems are the same, 
the same effective Lagrangians consistent with the symmetries of the strong interaction 
underlie the potentials to be used in the Lippmann-Schwinger equations for NN and irN, 
independently. All single pion production channels including higher partial waves are 
considered. 

The model used for the NN distortions in the initial and final states is based on the 
Bonn potential jH2j. The A-isobar is treated in equal footing with the nucleons through a 
coupled-channel framework including the A" A" as well as the iVA and AA channels. The 
model parameters were adjusted to the phase shifts below the pion production threshold. 

Since all the short range mechanisms suggested in literature to contribute to pion 
production in iViV collisions mainly influence the production of s-wave pions, in the Jiilich 
model only a single diagram was included (heavy meson through the uj as in diagram (d) 
of Fig. 12. lj) to parameterise these various effects. The strength of this contribution was 
adjusted to reproduce the total cross section of the reaction pp —> ppn° close to threshold. 
After this is done, all the parameters are fixed. 

The model describes qualitatively the data, as shown in Fig. 12.81 The most striking 
differences appear for double polarisation observables in the neutral pion channel. As a 
general pattern the amplitudes seem to be of the right order of magnitude, but show a 
wrong interference pattern. For charged pion production most of observables are described 
satisfactorily. In contrast to the neutral channel, the charged pion production was found to 
be completely commanded by two transitions, namely, 3 Pi — > 3 S±s, which is dominated 
by the isovector pion re-scattering, and 1 D2 ^ 3 Sip, which governs the cross section 
especially in the regime of the A resonance. 

2.3 Chiral perturbation theory 

In the late 90 's there was the hope that might resolve the true ratio of re- 

scattering and short-range contributions in pion production. It came as a big surprise, 
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Figure 2.8: Comparison of the Julich phenomenological model of jHUEIIllHI] t° the data. The 
solid lines show the result of the full model and the dashed lines show the results without the 
A contributions. Figure taken from Ref. [3J. 



however, that the first results for the reaction pp —>■ ppir |3T)| EBj showed that the %PT 
irN scattering amplitude interfered destructively with the direct contribution, making 
the discrepancy with the data even more severe, and thus suggesting a significant role 
for heavy meson exchanges in n° production [HU EE]- hi addition, the same isoscalar 
re-scattering amplitude also worsened the discrepancy in the ir + channel [66 . 

In low-energy pion physics, the constrains to an effective field theory (%PT) come 
from chiral symmetry, since it forces not only the mass of the pion to be low, but also 
the interactions to be weak: the pion needs to be free of interactions in the chiral limit 
for vanishing momenta. The first success of (xPt) was the application to meson-meson 
scattering jHZj. Treating baryons as heavy allowed straightforward extension of the scheme 
to meson-baryon|68. as well as baryon-baryon systems[nHl[I0J[7TJ[72]. The A isobar could 
also be included consistently in the effective field theory [75] . It was shown recently [741 
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I7H] that when the new scale induced by the initial momentum p ~ y/m n M 7 for meson 
production in nucleon-nucleon collisions is taken properly into account, the series indeed 
converges. 

The starting point for the derivation of the amplitude is an appropriated Lagrangian 
density, constructed to be consistent with the symmetries of QCD and ordered according 
to a particular counting scheme. The leading order Lagrangian is [HI EHl EE] 
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where /„. is the pion decay constant in the chiral limit, is the axial- vector coupling of 
the nucleon and Ha is the AA^-zr coupling. The the nucleon, pion and A field are N, ir and 
^Aj respectively. The isobar-nucleon mass difference is A, and the quark mass difference 
contribution to the neutron-proton mass difference is 5. The 1/2 — > 3/2 spin transition 
matrix operator is S and the 1/2 — > 3/2 isospin transition matrix operator is T. The are 
normalised such that 



1 



SiSj — — (25ij — ieijkdk) 



(2.13) 



The constants q, not constrained by chiral symmetry, depend on the details of QCD 
dynamics. They are at present unknown functions of the fundamental QCD parameters, 



7 Since \JMra v is smaller than the characteristic mass scale of QCD (Mqcd ~ lGeV), at least in the 
chiral limit, the contribution of other states (the Roper, the p meson, etc) can be buried in short-range 
interactions. 
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tree 
L i 


loop 


cf ee (/A) 


cl oop (A) 


1 


-0.64 


-0.93 


-0.64 


-0.93 


2 


1.78 


3.34 


0.92 


0.64 


3 


-3.90 


-5.29 


-1.20 


-2.59 


4 


2.25 


3.63 


0.90 


2.28 



Table 2.1: Low-energy constants a for and of Eq. (JZUJ-Eq. (ET^ (in GeV" 1 ). 
The first two columns are from Ref. jZZ] and Ref. [76J. The last two columns are obtained 
subtracting a chosen A contribution of = — = 2cf = 2.7GeV _1 [3|- 

with A without A 

-SK 

Figure 2.9: Illustration of resonance saturation. Figure taken from Ref. [3J. 

and can be extracted from a fit to elastic irN scattering 8 . The corresponding values are 
in Table l2~Tl In a theory without explicit A's, their effect is absorbed in the low energy 
constants (see Fig. 12. 9|) . which is called resonance saturation (infinitely heavy or static 
limit). Thus, in the case of a theory which considers explicit A's, the A contribution 9 
needs to be subtracted from the values given in the first two columns of Table 12.11 



8 In the standard case of processes involving momenta of order m v , the predictive power is not lost, 

because at any given order in the power counting only a finite number of unknown parameters appear. 

After these unknown parameters (low energy constants) are fitted to a finite set of data, all else can be 

predicted at that order |Hrj|. 

9 Note that there is some sizeable uncertainty in the A contribution[HHj- The empirical values of the 

low energy constants ci, C2, C3, C4 can be understood from resonance exchange. In particular, assuming 

that ci is saturated completely by scalar meson exchange (which is in agreement with indications from 

the NN force), the values for C2-C4 can be understood from a combination of A, p and scalar meson 

exchange EU: = -cf = 2cf = -2.54 - 3.10GcV _1 . 
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XX 

a) b) c) d) 

Figure 2.10: Diagrams contributing to the impulse term that are irreducible in the context 
of chiral power counting for the pp — > ppit reaction. The solid and dotted lines denote 
the nucleons and the pions, respectively. Leading(sub-leading) vertices are represented by 
solid(open) circles. They are originated by the three-momentum(energy) dependence of the 
vertex. 



2.3.1 Power counting for the impulse term 

Within the framework of x?T, only pion exchange is considered. Also, the impulse 
term (Fig. 12. ll (a)) is included in the class of irreducible diagrams defined in Weinberg's 
sense. A sub-diagram is considered reducible in Weinberg's sense if it includes a small en- 
ergy denominator of the order of ~ m^/M and irreducible otherwise |63j. In the following, 
we present the power counting for the impulse term. 

Since near threshold the pion carries an energy of the order of the pion mass, at 
least one of the nucleon intermediate states, before or after pion emission, must be off 
mass-shell by ~ m w . The transition to an off- mass-shell state is induced by a relatively 
high-momentum (~ y/m^Mj meson exchange mechanism. Therefore, the irreducible sub- 
diagrams of Fig. 12.11 (a) are in lowest order to be drawn as in Fig. 12.101 
The corresponding power counting is done as follows: 

- Close to threshold, the a ■ q term of the ttNN vertex is suppressed, and the pion- 
nucleon interaction proceeds via the Galilean term N^t -iz a • VN of Eq. (J2.12I) . 
This yields a factor of 

1 pE n ^Mm^m n mj 2 

U M ~ UM ~ UM 1 ' 2 ' 1 ' ' 

- Since the nucleon-nucleon interaction originates from virtual (static) pion exchange, 
as in Fig. 12.101 (a) and (b), each three- momentum dependent ttNN vertex con- 
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tributes with a factor of y/m n M / f n and the pion propagator contributes with a 
factor of (m^-M) -1 , thus giving the overall factor of f~ 2 from the nucleon-nucleon 
interaction. If the nucleon-nucleon interaction arises from exchange of a heavier 
meson, as in Fig. 12.101 (c) and (b), the overall factor is also ( i 



2 



Near threshold, the contribution of the non-relativistic two-nucleon propagator is 

(-^intermediate — -^initial) ~ m 7r 1 '■ 

Finally, the overall contribution of the impulse term is then 

rnj 2 1 1 1 fm^ 

x ^— =«A/TT- ( 2 - 15 ) 



f % M 1 / 2 f*m w /3y M 

In conclusion, the impulse term is not irreducible in the Weinberg's sense. Within %PT 
the impulse term is calculated as given by the (distorted) diagrams in Fig. 12.101 

2.3.2 Why 7r° is problematic 

For neutral pion production there is no meson exchange operator at leading order and 
the nucleonic operator gets suppressed by the poor overlap of the initial and final state 
wave functions (an effect not captured by the power counting) and interferes destructively 
with the direct production of the A. Thus the first significant contributions appear at 
NNLO. As there is a large number of diagrams at NNLO, the different short-range mech- 
anisms found in the literature and discussed before are of similar importance and capable 
of removing the discrepancy between the Koltun and Reitan result and the data. Charged 
pion production is expected to be significantly better under control, since there is a meson 
exchange current at leading order and there are non- vanishing loop contributions [66 . 

The importance of the pion loops 

The most prominent diagram for neutral pion production close to threshold is the 
pion re-scattering via the isoscalar T- matrix that, for the kinematics given, is dominated 
by one-sigma exchange [2]. Within the effective field theory the isoscalar potential is built 
up perturbatively. The leading piece of the one-sigma exchange gets cancelled by other 
loops that cannot be interpreted as a re-scattering diagram and therefore are not included 
in the phenomenological approaches. This is an indication that in order to improve the 
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Figure 2.11: Reduced cross section for pp — > pnn + with (left panel) and without (right panel) 
the A contribution. The dotted and solid lines correspond, respectively to the Reid 93 and 
Argonne 18 NN potentials. Figure taken from Ref. [66J. 

phenomenological approaches, at least in case of neutral pion production, pion loops 
should be considered as well. 

2.3.3 Charged pion production in xPT 

The main calculations on pp — > pnir + in the framework of xPT are those of Ref. |66j . 
which also included the mechanism proposed in Refs. jUlEE], where the short-range in- 
teraction is supposed to be originated from Z-diagrams mediated by a and u changes. 
The Coulomb interaction is disregarded in the pnir + final state. The channels considered 
were 3 Pi — » ( 3 S±)p and 3 Pq — > CSq) s. For the 3 Si final state, the Weinberg- Tomozawa 
term contribution is found to be the largest. Most of the other contributions were much 
smaller and tended to cancel each other to some extent. The exception was the A con- 
tribution which had a significant destructive interference with the Weinberg- Tomozawa 
term. Although the theory produces the correct shape for the r\ dependence, it fails in 
magnitude by a factor of ~ 5 (see Fig. 12.11]) . 

Ref. [EH! suggested that the A contribution may be overestimated by a factor of 50% 
or more. Actually, the cross sections with and without the A contribution (respectively, 
left and right panel of Fig. I2.11|) . differed by almost a factor of 2. This could arise from 
the uncertainties on the 7rA^ A coupling constant, or from the neglecting of the A- A" mass 
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difference in energies. The kinetic energy of the A was also neglected, account of which 
would further decrease its amplitude. 

2.3.4 The Mow k approach 

As there are no xPT NN potentials available yet for pion production calculations, in 
practices one usually uses a hybrid (not much consistent) %PT approach, in which the 
transition operators are derived from %PT but the nuclear wave functions are generated 
from high-precision phenomenological NN potentials. However, a conceptual problem 
underlying these hybrid xPT calculations is that, whereas the transition operators are 
derived assuming that relevant momenta are sufficient small compared with the chiral 
scale x ~ lGeV, the wave functions generated by a phenomenological potential can in 
principle contain momenta of any magnitude |78j. 

A systematic method based on the renormalisation group approach was recently 
developed |79j to construct from a phenomenological "bare" NN potential an effective 
NN potential, Vi w-k, by integrating out momentum components above a specified cutoff 
scale A. For a A = 2.1fm _1 , it was found[?n] that the low momentum behaviour of the 
NN wave functions calculated from V[ ow _k is essentially model independent. 

The very recent first calculation [7H] for pp — > ppw using the VJ w-k approach led 
to a cross section closer to the experimental data than the one calculated with "bare" 
potentials. 

2.4 Energy prescription for the exchanged pion 

In most of the calculations on pion production, several approximations to the pion 
production vertex and the kinematics have been tacitly assumed from the very first investi- 
gation by Koltun and Reitan[33]. However, these prescriptions were found to have a signif- 
icant effect both in the magnitude and energy dependence of the cross section |5? ] l59 | 180] . 
as it is illustrated by Fig. 12.121 and by Fig. 12.1 HI 

Fig. 12.121 refers to the calculations [HE] for pp — > ppir based on the Jiilich model 
presented in Sec. 12.21 It differs from the calculation of Ref. [SJ since the momentum of 
the exchanged pion is not fixed to the value corresponding to the particular kinematic 
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Figure 2.12: Cross section for pp — > ppir . The solid curve is the full calculation of the re- 
scattering and direct production diagrams whereas the dashed line corresponds to the direct 
production diagram only. The dashed-dotted curve is the full result scaled by a factor of ~ 4. 
The dotted curve is the full calculation with the approximate treatment of the pion production 
vertex of Ref. j5TJ (frozen kinematics). Figure taken from Ref. [59J. 
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Figure 2.13: Cross section for pp — > ppir , including both impulse and re-scattering ampli- 
tudes. The solid line is the calculation assuming four-momentum conservation and that the 
intermediate particles are on-mass-shell. The dotted and dashed-dotted lines correspond to 
two different fixed-threshold prescriptions for the energy of the exchanged pion. Figure taken 
from Ref. (ED]. 

situation at the pion production threshold when there are no distortions in the initial and 
final pp states. 

The full treatment of the pion production vertex (solid line of Fig. I2.12|) reduced the 
7T° production rate by a factor of 2. Therefore, the work of Ref. concluded that the 
enhancement of the re-scattering amplitude due to offshellness falls a bit too short to 
explain the scale of the pp — *> pptr cross section. 

In Fig. 12.131 we show a calculation for the pp — > ppn° reaction within %PT, which also 
aimed to investigate the effect of the simplifying assumptions on the energy-momentum 
flow in the re-scattering diagramJSU]. For this case, the discrepancy between the approx- 
imations was found to be of a factor of 3. 

A clarification of these formal issues is thus necessary before one can draw conclusions 
about the physics of the pion production processes. More specifically, one needs to obtain 
a three-dimensional formulation from the more general Feynman procedure. This will be 
the subject of the next Chapter, and was the starting point of this thesis work. 
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Abstract: DWBA calculations, which are currently applied to pion production, are 
based on a quantum-mechanical three-dimensional formulation for the initial- and final- 
A^A^ distortion, which is not obtained from the four-dimensional field-theoretical Feynman 
diagrams. As a consequence, the energy of the exchanged pion has been treated approxi- 
mately in the calculations performed so far. This Chapter will discuss the validity of the 
DWBA approach through the link with the time-ordered perturbation theory diagrams 
which result from the decomposition of the corresponding Feynman diagram. 
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All the meson production mechanisms described in the previous Chapter are derived 
from relativistic Feynman diagrams. Nevertheless, within the framework of DWBA, the 
evaluation of the corresponding matrix elements for the cross section proceeds through 
non-relativistic initial and final nucleonic wave functions. Therefore the calculations apply 
a three-dimensional formulation in the loop integrals for the nucleonic distortion, which 
is not obtained from the underlying field-theoretical Feynman diagrams. Namely, the 
energy of the exchanged pion in the re-scattering operator (both in the 7riV amplitude 
and in the exchanged pion propagator) has been treated approximately and under different 
prescriptions in calculations performed till now. In other words, the usual approach, as 
in Refs. jSHl EH EH E3 EH EJ] , is to use an educated "guess" for the energy- dependence 
of the virtual pion-nucleon interaction and then to use a Klein-Gordon propagator for 
the exchanged pion propagator. In particular, most of the calculations assume a typical 
threshold kinematics situation even when the exchanged pion is expected to be largely 
off-shell. 

A theoretical control of the energy for the loop integration embedding the non-relativis- 
tic reduction of the Feynman 7r-exchange diagram in the non-relativistic nucleon-nucleon 
wave functions became then an important issue, as pointed out in Refs. [5^1 183j. 

In this Chapter, following the work of Ref. |84j . we deal with the isoscalar re-scattering 
term for the pp — > ppir reaction (near threshold). Although for tt° production, as men- 
tioned in Chapter the re-scattering mechanism is indeed small, the amount of its 
interference with the (also small) impulse term depends quantitatively on the calculation 
method. We note, furthermore, that the pion isoscalar re-scattering term, which is energy 
dependent, increases away from threshold, and that for the charged pion production re- 
actions the isovector term is important, and also depends on the exchanged pion energy. 
Thus the knowledge gained from the application discussed here to the pp —* ppir reaction 
near threshold is useful for other applications. 

Specifically, this Chapter will focus on the investigation of 

i) the validity of the traditionally employed DWBA approximation. We will use as 
a reference the result obtained from the decomposition of the Feynman diagram 
into TOPT diagrams, and realise how the last ones link naturally to an appropri- 
ate quantum- mechanical DWBA matrix element. This study for a realistic irNN 
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coupling was not done before; 

ii) the (numerical) importance of the three-body logarithmic singularities of the exact 
propagator of the exchanged pion, which are not present when the usual approxi- 
mations in DWBA are considered; this study was not done before. 

The work of Ref. [HI] generalised the work of Refs. [H21 EH] which considered a solvable 
toy model for scalar particles and interactions and treated the nucleons as distinguishable 
and therefore pion emission to proceed only from one nucleon. It is described by the 
following Lagrangian 

i=l,2 ^ ' 

9 fN\N 2 Ti + g a J2 N}N,a + £ £ iV/iV, (<V) 2 (3.1) 

J* i=l,2 ^ i=l,2 

where Ni and ir are the nucleon and pion field, respectively. Eq. (jH.lj) further as- 
sumes a Yukawa coupling for the pion-nucleon coupling, ^jl an d describes the pion re- 
scattering through a irN seagull vertex inspired by the chiral irN interaction Lagrangian 
of Eq. (12.121) . jsQnQo, where c is an arbitrary constant. In this toy model, the nuclear 

J TV 

interactions are described through u-exchange, which also couples to the nucleons via 
a Yukawa coupling, g a , and the nucleons are treated non-relativistically (in particular, 
contributions from nucleon negative-energy states are not considered). 

There are some nonrealistic features of this model that also motivated our work and the 
one of Ref. [HI] . First of all, the model does not satisfy chiral symmetry, which would have 
required a derivative coupling of the pion to the nucleon instead of the simpler Yukawa 
coupling. Furthermore, in this toy model the nucleons, interacting via a exchange, have 
a stronger overlap than in a more realistic model, since it does not include short-range 
repulsive nucleon-nucleon interactions that keep the nucleons apart. Note also that near 
threshold kinematics, the scalar particle is produced in a s-wave state, as is the final 
nucleon-nucleon pair (and the initial nucleon -nucleon pair too, by angular momentum 
conservation requirements). This is not the case for the production of a real pion, since it 
is a pseudo-scalar particle and, therefore, near threshold the production of s-wave pions 
calls for a initial P-wave nucleon state. 

Our calculation employs a physical model for nucleons and pions and investigates how 
much of the features found in Refs. jH2 EH] survive in a more realistic calculation which 
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uses a pseudo- vector coupling for the 7r iViV vertex, the xPT amplitude (SHI and the 
Bonn B potential jH2j for the nucleon-nucleon interaction. 



3.1 Extraction of the effective production operator 

3.1.1 Final-state interaction diagram 
The amplitude 

The Feynman diagram for the reaction pp — > ppir where the NN final-state interaction 
(FSI) proceeds through sigma exchange is represented in part a of Fig. 13.31 After the 
nucleon negative-energy states are neglected, it corresponds to the amplitude 

Mra= i3i JSf fm Q?-l l+ ieQ<> 2 -l + is K$ -I'+ie (3 - 2) 

1 



Q'o 2 -ul + ie 

where the exchanged pion has four- momentum Q' = (Q' ,q'). The four-momenta of the 
intermediate nucleons are Q\ and Q 2 and the four-momentum of the exchanged o is 
K. The on-mass-shell energies for the intermediate nucleons, exchanged pion and sigma 
mesons are ou^2, u-k, u a , respectively. In Appendix El details on these functions are given. 
All quantities are referred to the three-body centre-of-mass frame of the 7rNN final state. 

In Eq. (13. 2|) / (Q' ) is a short-hand notation for both the M. w nn and M^nn- We note 
that / in fact may depend on the three-momenta and energies of both the exchanged and 
produced pion, and also on the nucleon spin. However, for simplicity, only the dependence 
on the energy of the exchanged pion, Q' Q , is emphasised, because this is the important 
variable for the main considerations of these Chapter. For the toy model described by the 
lagrangian of Eq. (|3.1|) , it is 

/ (Q'o) = (t) x (^Q'oK) , (3.3) 



\ fit J \ fn 

but for the more realistic case of Ref. [Hlj which uses a pseudo-vector coupling for the 
V-kNN vertex and the xPT iiN amplitude (HH], it reads 



f (Q'o) = (jf ■ X ^ [ 2Cl + [c 2 + C 3 



c 3 - 



8M J ml ml 
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MQ'o} 



E% — U>2 

o 



E 2 -F 2 - u a 
o 



Re{Q' } 



o 

E2-F2 + UJ a 



o 

E % + uj x - E\ 



Figure 3.1: Schematic position (scaled for a kinematics close to pion production threshold) 
of the poles in Q' of the amplitude of Eq. ()3.5|) . The diamonds, open circles and bullets 
represent, the nucleon poles, the sigma poles and the pion poles, respectively. 



In Eq. (J3.3J) and Eq. ()3.4|) . E^q^) is the energy (three- momentum) of the produced pion. 
As mentioned before, f n is the pion decay constant and qa is the nucleon axial vector 
coupling constant, which is related to the pseudovector coupling constant J'^nn by the 
Goldberger-Treiman relation M- = ^ NN . The low-energy constants c's are defined on 
Table O 



Using four-momentum conservation, the amplitude of Eq. ()3.2}) can be written as 

Mfsi = ~ %gl \ ^ f m f To f u + T F a f\u iF (3 - 5) 

1 111 

F 2 - E 2 + Q' - uj a + ie F 2 - E 2 + Q' + uj a - ie Q' - uj n + ie Q' Q +u 7r — ie 

where the functions E\ 2 {F\p) stand for the energy of the initial (final) nucleons. 

The integrand in Eq. (J3.5|) has three poles in the upper half-plan and three poles in 
lower half-plan. They are schematically reproduced on Fig. 13.11 

To integrate Eq. ()3.5|) over the energy variable Q' Q we can close the contour on one 
of the half-planes and pick each of the three poles enclosed. However, to get a more 
straightforward connection with the DWBA formalism through time-ordered perturbation 
theory, it is better to perform a partial fraction decomposition to isolate the pion poles 
before integrating over Q' . All details of this decomposition are given in Appendix IU1 
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h 

k 



□ 




Figure 3.2: Schematic position of the poles in Q' resulting from the partial fraction decom- 
position of the amplitude of Eq. (13. 5j) . The meaning of the symbols is the same of Fig. 13.11 



After the Q' integration, one obtains for Eq. (|3.5J) : 

a F SI r) f d Q 1 
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with E tot = 2E = Fi + F 2 + E n and E = E 1 = E 2 . 

This equation evidences that there are six contributions to the amplitude. These six 
terms, originated by the four propagators of the loop, can be interpreted as time-ordered 
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a CL\ CL2 CI3 04 <2 5 0,6 

Figure 3.3: Decomposition of the Feynman diagram in terms of six time-ordered diagrams 
for the final-state interaction. The pion(sigma) field is represented by a dashed(solid double) 
line. The nucleons are represented by solid lines. The DWBA amplitude may be identified to 
the first four time-ordered diagrams (g^ to a 4 ). The last two diagrams (a 5 to a 6 ) are usually 
called stretched boxes. 

diagrams. They are represented by diagrams a\ to a§ in Fig. 13.31 This interpretation 
justifies the extra subscript label TOPT for the M.fsi amplitude in Eq. ()3.6|) . 

Obviously, the result of direct integration is the same of the result obtained doing a 
partial fraction decomposition before integrating. For the / (Q' Q ) functions considered, 
since they have a simple 1 (linear) dependence in Q' Q , the integrals in which the two or 
three poles are all in the same half-plan vanish (first four diagrams of Fig. 13 .2jl . and we 
end up with only six terms, corresponding to four distributions of the poles (four last 
terms of Fig. 13. 2|) . 

Because the partial fraction decomposition of the propagators in Eq. (|3.5|) was done 
prior to the integration over the variable Q' Q , we have the following outcome which is 
independent of the choice of the contour of this integration: the only terms from the 
decomposition which do contribute to the integral correspond to the ones with only one 
pole, which happens to be the Q' = lu w = a/ m 2 + q 12 or the Q' = —LU n = — a/ m 2 + q 12 
pion poles. The other terms, with nucleon poles and/or sigma poles, together or not with 
pion poles, have all these poles located on the same half-plane and consequently their 
contribution vanish (upper part of Fig. 13. 2[) . 

We stress at this point that this method for the energy integration implies effectively 
that the irN re-scattering amplitude is evaluated only for on-mass-shell pion energies. In 
this way, off-shell extrapolations which are not yet solidly constrained are avoided, which is 

1 The linear dependence in Q' Q of the functions / (Q' ) guarantees that the integral over the curve of 
radius R which closes the contour vanishes when R — > 00. 
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Figure 3.4: Schematic position of the poles in Q' resulting from the partial fraction decompo- 
sition of the amplitude of Eq. (13. 5 j) for the DWBA and stretched boxes terms. The meaning 
of the symbols is the same of Fig. 13.11 



an advantage. Other methods may need the contribution of the off-shell amplitude in the 
integrand with the form shown in Eq. ()3.5j) . But the net result is the same, provided that 
all the contributions from all (nucleon, sigma and pion) propagator poles are considered. 
So far, calculations [SHI EH1 EH EH EM EJJ did not consider the pion propagator poles, 
since they approximate that propagator by a form free of any singularity. Consequently, 
they exhibit a strong dependence on the irN amplitude at off-mass-shell energies of the 
incoming pion. 

The effective pion propagator for FSI distortion 

From the six terms in Eq. (jH.fij) . the first four terms (corresponding to diagrams a\ to 
of Fig. 13. 3|) have the special feature that any cut through the intermediate state intersects 
only nucleon legs. Thus, they may be identified to the traditional DWBA amplitude for 
the final-state distortion. In contrast, in the last two diagrams as and a§ of Fig. 13.31 
any cut through the intermediate state cuts not only the nucleon legs, but also the two 
exchanged particles in flight simultaneously. They are called the stretched boxes [52]. The 
corresponding pole distribution is illustrated in Fig. 13.41 

Because of the identification of diagrams a\ to with DWBA, we may collect the 
four first terms of Eq. ()3.6|) and obtain what we may call the reference expression for the 
DWBA amplitude: 



M 



fsi _ 1 f A' 



DWBA 



1 



/ M j^—^ — b : 7T T nn- (3-7) 



2 J (2tt) 3 \- v ' J {E 1 + E 2 -E 7T -oo 1 - uj 2 ) 
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Here T stands for the transition-matrix of the final-state interaction. We have used both 
sigma exchange, (with m a = 550MeV) as in Ref. jH2j, which makes Eq. ()3.7j) coincide 
exactly with Eq. (j3.6|) . and also the T-matrix calculated from the Bonn B potential j 
For the cx-exchange case, one gets 

1 1 



n 2 

yDWBA _ Ha 



(3. 



. (E to t — Fx — — U2 — LOa) {E t ot — F 2 — E w — cox — u a 
The details of the T-matrix calculation are in Appendix [DJ I n the derivation of the 
integrand of Eq. ()3.7)1 the propagators for the two nucleons in the intermediate state 
fused into only one overall propagator with the non-relativistic form, 

Gnn = (TP I p w ~, — TTv ( 3 - 9 ) 

The function / (uj^) includes the contribution of the two pion poles u n and — u n corre- 
sponding to two different time-ordered diagrams, 

~ _ f (ay) (Ei - E n - uh - ay) + / (-ay) (E 2 - u 2 - ay) (3 10) 

UJn 

where the function / (ay) is the product of the irN amplitude with the wNN vertex. The 
multiplicative kinematic factors // = ^-^-^-^ anc i ff = e 2 -loi-uj 2 ma y ^ e treated as 
form factors (see Fig. 13.61) . 

In the derivation of Eq. ()3.7|) from Eq. ()3.6|) the function G n for the pion propagator 
turns to be exactly 

Gn ~ [™ + f] 2 - [(£ to4 -£-f)-^-oy] 2 ' (3 ' U) 
which gives the form of the effective pion propagator appropriate for a DWBA final-state 
calculation, and can also be written as 

Gy — - — — , : -. (3.12) 

[ui + E n - Ei + ay J [E 2 -oj 2 - UJ n \ 

The same approach will next be applied to the ISI case. 

3.1.2 Initial- state interaction diagram 
The amplitude 

The corresponding Feynman diagram (b in Fig. I3.5|) . for the NN initial-state interaction 
(ISI) when it proceeds through sigma exchange, after the nucleon negative-energy states 
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Figure 3.5: Decomposition of the Feynman diagram in terms of six time-ordered diagrams for 
the initial-state interaction. The pion(sigma) field is represented by a dashed(solid double) 
line. The nucleons are represented by solid lines. The DWBA amplitude corresponds to the 
first four time-ordered diagrams (fei to 64) and the stretched boxes to the last two (6 5 to b e ). 



are neglected, generates the amplitude: 

M ISI = -€/^/W;)77-^ ~- F4 _ F ly T- ( 3 - 13 ) 

J (27r) Qq + F 2 - u 2 + is F x + E n - Q' -u x + ie 

1 111 

Q' Q - E 2 + F 2 - u a + ie Q' - E 2 + F 2 + uj c - ie Q' - u n + ie Q' Q + u„ - ie ' 

where we have used a notation analogous to one used for the final-state amplitude of 
Eq. JH. 

As before, in order to perform the integration over the exchanged pion energy Q' Q , a 
partial fraction decomposition to isolate the poles of the pion propagator was done. By 
closing the contour such that only the residues of the Q' = ±u n = ±^/g' 2 + ml poles 
contribute, in an entirely similar way to the FSI case, one obtains the amplitude: 
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The six terms in Eq. ()3.14jl are interpreted as contributions from time-ordered diagrams 
represented in Fig. 13. 5) b\ to b§. This interpretation justifies the extra subscript label 
TOPT for the Ai 1 amplitude. Although Eq. ()3.6|) and Eq. ()3.14|) are formally alike, for 
the initial state an extra pole is present (besides that from the nucleons propagator), since 
it is energetically allowed for the exchanged pion to be on-mass-shell. All the singularities 
are decisive for the real and imaginary parts of Eq. (|3.fi|) and Eq. (|3.14j) . 

The effective pion propagator for the ISI distortion 

Analogously to the final-state case, the first four terms in Eq. (|3.14jl (diagrams b\ to 64 
of Fig. 13.5) where any cut of the intermediate state intersects only nucleon lines) are 
identified with the DWBA amplitude for the initial-state distortion, and the last two 
(diagrams 65 and b§ of Fig. 13.51 with two exchanged particles in flight in any cut of the 
intermediate state) to the stretched boxes. In other words, the decomposition obtained 
in Eq. (|3.14|) allows to write the exact or reference expression for the DWBA amplitude 
for the initial-state distortion, by collecting the four first terms, which have intermediate 
states without exchanged particles in flight, i.e., 



M I8J __WiV 

M DWBA - 2 J s 



fM°* 7^—^ (3-15) 

J [hi + t, 2 — ui — u 2 ) 



where 



J Kr) = • (3-16) 



Again, / (u n ) stands for the product of the ttN amplitude with the 7iNN vertex. As 
before, the multiplicative factors f[ = E ^+ F ^~^~^ anc l J* — f 2 -u2-^ ^ ave a form-factor 
like behaviour (see Fig. 13.6)1 . cutting the high momentum tail. 

In Eq. (13.15)1 we made also the replacement of the exchanged particle potential by the 
A^A^ scattering transition-matrix. The two-nucleon propagator 

Gnn = 7^ — — ^ r (3-17) 

is the non-relativistic global NN propagator. As we did for the final-state interaction, we 
extracted from Eq. (j3.15j) the effective pion propagator to be included in the initial state 
distortion in a DWBA-type calculation. It reads, 

Gw = [^ + ^ + f] 2 -[(f -^ + ^)-^] 2 ' (3 ' 18) 
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Figure 3.6: Form-factor like behaviour of the kinematic multiplicative factors for the simple 
case when / (u n ) = —f (—u) w ), as a function of the exchanged pion momentum. The dotted 
line corresponds to threshold and the solid line to the maximum momentum of the emitted 
pion for T iah = 440MeV. 



or, 



Q — 

(F 2 -u 2 - uj n ) (Fi + E^-ujx- u n ) ' 
For the a-exchange case, one gets 

,2 



DWBA 9c 



2cu n 



1 



1 



(Etot -E x -u 2 - v a ) (E tot - E 2 - uji - cv a ) 



(3.19) 



(3.20) 
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Figure 3.7: Importance of the stretched boxes relatively to the DWBA FSI amplitude M-dwba 
of Eq. (13. 7 j) . as a function of the mass of the scalar particle for the final NN interaction. The 
7iN amplitude is contact re-scattering vertex as in Eq. (13. 3|) (solid line) and the xPT amplitude 
as in Eq. (13. 4j) (dashed line). The energy is taken at the pion production threshold. Absolute 
values of the amplitudes are considered. 

3.2 Stretched Boxes vs. DWBA 

In this section, the iViV and (NN)w channels considered in all calculations refer to the 
transition 3 P — > ( 1 *S'o) s, which is the dominant one for pp — > ppir (see Chapter |3J). The 
details of the partial wave decomposition analytical formulae are given in Appendix lEl 

The DWBA amplitude was found to be clearly dominant over the stretched boxes in 
the realistic model considered, as the dashed line in Fig. 13 . 71 documents for the FSI case. 

It is interesting to compare this result with the one of Ref . [H2] • As shown in Fig. 13 .7\ 
the stretched boxes amplitude is less than 1%, of the total amplitude and therefore is 
about 6 times more suppressed than in the dynamics of the toy model used in that 
reference. Replacing the irN amplitude from %PT by a simple contact amplitude, the 
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Figure 3.8: Importance of the stretched boxes compared to the DWBA (FSI) for the cross 
section as a function of the mass of the scalar particle for the final NN interaction. The nN 
amplitude is the %PT amplitude (left) and a contact re-scattering vertex (right). 

stretched boxes amplitude becomes slightly more important, but they still do not exceed 
4% of the DWBA amplitude (solid line in Fig. ET7J) . 

In terms of the cross section, the weight of the stretched boxes relatively to DWBA is 
even smaller, of the order of 0.02% at most. This is seen in Fig. IH.8I where we compare, 
for three different values of the laboratory energy, the cross section obtained with only 
the DWBA contribution {(Jdwba), with the one obtained with only the stretched boxes 
terms (a stretched) ■ In both cases considered, V^ N + V*^N an< ^ V-kNN + contact (left and 
right panel in Fig. K18[ respectively) the cross section with the DWBA terms is clearly 
dominant, and in a more pronounced way when compared to the less realistic case of 
Ref. jH2], n ot only at threshold but even for higher energies as 440MeV. 

However, we note that the ratios and energy dependence of the amplitudes and of the 
cross sections are significantly influenced by the irN amplitude used in the calculation. 
The stretched boxes are seen to be amplified by the contact irN re-scattering amplitude, 
due to an interplay between the ttN and the A^A" amplitudes. Relatively to the more 
realistic %Pt amplitude, the contact ttN amplitude gives a larger weight to the low- 
momentum transfer. The realistic vrA" amplitude satisfies chiral symmetry. This implies 
cutting small momenta and giving more weight to the region of large momentum transfer, 
which however in turn is cut by the nucleonic interactions. The difference between the 
two nN amplitudes is clearly seen by comparing the behaviour of each curve on the left 
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Figure 3.9: The same as Fig. O but for the DWBA ISI amplitude M^ba of Eq. (I3~T5|) . 

panel of Fig. 13.81 with the corresponding curves on the right panel, for small values of the 
NN interaction cut-off. 

In the case of the initial-state amplitude, the stretched boxes amplitude is also much 
smaller than the DWBA amplitude for the two cases shown in Fig. 13.91 

The cross sections obtained with only the stretched boxes terms were found to be less 
than 1.2% of the DWBA cross sections, even for laboratory energies as high as ~ 
440MeV (see Fig. EHUl) . 

The results presented in Fig. 13.71 to Fig. 13.101 justify the DWBA treatment for pion 
production. 
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Figure 3.10: The same as Fig. 13.81 but for the ISI case. 



3.3 The logarithmic singularity in the pion propaga- 
tor (ISI) 

From the evaluation of the effective pion propagator for the ISI case of Eq. (I3.19J1 . at 
threshold, 

°* = (l^) 2 _ (^-^) 2 = u> n (rn*-u„y (3 ' 21) 
it is straightforward to conclude that G n has a pole, since it is energetically allowed for 
the exchanged pion to be on- mass-shell (u n = m n ). The treatment of the pole leads to 
a moving singularity which was not considered in the traditional DWBA calculations of 
Refs. |2H1 E3 EH HUH ETJ • It requires a careful numerical treatment which was included 
in our calculations. 

This is done by rewriting G n of Eq. ()3.19|) as 

1/1 1 



with 



Cei — Ct2 \Otl — UJ W Ct 2 — UJtt 

oli = F 2 — to 2 

a 2 = Fi + - u\. 



(3.22) 

(3.23) 
(3.24) 



With the definitions Fi j2 , loi >2 an d ^-k introduced in Sec. 13.11 (see Appendix |Bj) , then 
Eq. flEZZ) reads 



G-k 



ai — a 2 



a>i + u> n 1 



a 2 + u n 1 



Qnh (a? 2 ) Vi - X\ qj 2 (x 2 ) y 2 - x x 
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where 7 2 (#2) = \0u ~ Qk\ — 2q u q^x 2 + q\- The angles involved in the partial wave 

decomposition performed for the amplitudes are x\ = A. (q u — qt, qt) and 22 = A- (qu, qt) 
(see Appendix |Ej) . Here, y^, roots of the denominator of G n , are given by: 

2 o2 2 

Vi = a \~ with {3 2 = ml + % + I 2 {x 2 f and z = 1, 2 (3.26) 

(x 2 ) 4 

Once is written in this form, for the calculation of the partial wave decomposition, 
we have to deal with its poles at X\ = y^. We used the subtraction technique (see 
Appendix |FJ) exposing the three-body logarithmic singularity 2 : 

r 1 _mqL PL(Xi)dXi = pv f'/MdM Pi(ii)+2/fci!lftU _ 

J-i Vi-xi + is J -1 Vi - x 1 

-iirP L {y i )f(y i ,x 2 ) (3.27) 

where Pl and are the Legendre polynomials of order L and Ql are the Legendre functions 
of the second kind of order L[H3]. These last functions exhibit logarithmic singularities 
which are given by the condition yi = ±1, which are determined analytically. Then, the 
integral over the moving logarithmic singularities is handled by a variable mesh. When the 
number of singularities (always between zero and two) is n s , the interval of integration 
is divided into 2n s + 1 regions including as breaking points those given by the found 
singularities. In each one of the regions we considered a Gaussian mesh. 



3.4 Conclusions 

In this Chapter we have generalised the work done in Refs. |H2*| 183] on the study of 
the DWBA approaches for the irreducible pion re-scattering mechanism of the reaction 
pp — > ppn . We considered a physical model for nucleons and pions, combining a pseudo- 
vector coupling for the 7rNN vertex and the x?T 7riV re-scattering amplitude. We treated 
the initial and final state distortion exactly, taking into account in the first case the three- 
body moving singularities in the pion propagator. 

Both for the final- and initial-state interaction, our results show that the DWBA 

formalism is quite adequate at threshold and even at higher energies, since this part 

of the full amplitude is clearly seen to be dominant over the stretched boxes. This is 

independent of the model for the ttN re-scattering amplitude. Nevertheless, relatively to 
2 An alternative treatment may be found in Rcf. |83| . 
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less realistic models, a chirally constrained amplitude reinforces even more the relative 
importance of the DWBA amplitude in the total amplitude. 

The amplitudes given by Eqs. ([3.7[) and (j3.15|) obtained as non-relativistic reductions of 
the corresponding Feynman diagrams for the final- and initial-state interactions, can now 
be employed to study the commonly used prescriptions for the energy of the exchanged 
pion, both in the propagator and at the re-scattering vertex. This energy, not fixed 
in a non-relativistic quantum-mechanical calculation, was fixed by TOPT to be u^, its 
on-mass-shell energy. 

Since in the time-ordered diagrams energy is not conserved at individual vertices, each 
of the re-scattering diagrams for the initial- and final-state distortion defines a different off- 
energy shell extension of the pion re-scattering amplitude, as calculated in Sec. 13.1. II and 
Sec. 13.1.21 As a consequence, the two operators for FSI and ISI are different, demanding 
the evaluation of two different matrix elements of such operators between the quantum- 
mechanical NN wave functions. In the next Chapter, an alternative approach will be 
presented, where this problem does not arise. 
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Chapter 4 

The S-matrix approach for quantum 
mechanical calculations 



Contents 







4.1.1 Factorisation of the effective re-scattering ppera^pi 


4.1.2 The S-matrix technique 






4.2.1 ExDansion of the effective oion Dropaeator 


4.3 Conclusions 



Abstract: The pion re-scattering process is certainly part of the pion production 
mechanism, but its importance relative to other contributions varies considerably depend- 
ing on the approximations made in evaluation of the effective operators. In this Chapter, 
we will re-examine the nature and extent of this uncertainty. The pion re-scattering 
operator originated from the S-matrix construction, previously applied only below pion 
production threshold, is seen here to approximate well the exact DWBA result consistent 
with time-ordered perturbation theory. 
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The S-matrix approach for quantum mechanical calculations 



4.1 The pion re-scattering operator in the S-matrix 
technique 

To derive the quantum mechanical effective operators for pion production, and other 
effective nuclear operators in general, one starts from the relativistic (effective) Lagrangian 
written in terms of hadronic fields. The interactions, mediated by meson exchanges before 
and after the production reaction takes place, are included in the effective NN (and 
nucleon-meson) interaction, while from the irreducible parts connected to the reaction 
mechanism (e.g., pion production) one obtains effective operators, whose expectation 
values are to be evaluated between the initial and final nucleonic wave functions. One 
aims to get such effective operators consistent with the realistic description of the NN 
interaction, which can then be used in studies of the corresponding reactions not only in 
the simplest (one or two-nucleon) systems, but preferably also in heavier nuclei. 

The covariant techniques based on the Bethe-Salpeter equation or its quasipoten- 
tial re- arrangements are these days practically manageable only below meson production 
threshold. However, above the threshold the dressings of the single hadron propagators 
and interaction vertices via the meson loops have to be included explicitly. For this reason 
the construction of the production operator is so far realised mostly in the Hamiltonian 
quantum-mechanical framework (usually non-relativistic, or, with leading relativistic ef- 
fects included perturbatively by means of the decomposition in powers of p/m, where p 
is the typical hadronic momentum and m is the nucleon mass). 

The derivation of the nuclear effective operators below the meson production thresh- 
old within the Hamiltonian framework - leading to hermitian and energy independent 
NN and 3N potentials, and conserved electromagnetic and partially conserved weak cur- 
rent operators - can be done in many different ways (see discussion in Refs. |Ho^ |H3 
IHHl IHHj and references therein). At the non-relativistic order the results are determined 
uniquely. As for the leading order relativistic contributions, they were shown to be uni- 
tarily equivalent |89j . The unitary freedom allows to choose the NN potentials to be static 
(in the CM frame of the two-nucleon system) and identify them with the successful static 
semi-phenomenological potentials. 

Also, above the threshold energy the static limit is commonly employed, since more 
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elaborate descriptions which include the mesonic retardation and loop effects are techni- 
cally considerably more complex [001 EU IH2j, especially for systems of more than two 
nucleons. Both the static approaches and the ones including "retardation" typically 
consider contributions of several one-meson exchanges, and the potentials are fitted to 
describe the data. It is therefore difficult to assess how well do they approximate the 
covariant amplitudes (corresponding to the same values of physical masses and coupling 
constants) which are so far outside the scope of existing calculational schemes, but which 
we believe do provide in principle the consistent description of the considered reactions. 

Thus, the ultimately exact approach to the description of pion production (and in par- 
ticular of the pion re-scattering contribution) would be either the covariant Bethe-Salpeter 
or quasipotential frameworks (extended above the pion threshold) or the quantum me- 
chanical coupled-channel technique including retardation. In these approaches one has to 
treat consistently the non-hermitian energy- dependent NN interaction (fitted to the data 
also above pion production threshold) and consider besides the effects of renormalisation 
of vertices, masses and wave functions via meson loops. 

In this Chapter (following Refs. jS2 | 183 } l84" | I93j) we rather numerically estimate the 
range of the predictions from several frequently used simplifying approximations, and 
compare them to the result obtained from the reduction of the corresponding covariant 
Feynman diagrams for the pion re-scattering operator. This reduction coincides with the 
time-ordered perturbation theory |H4"j. To this end we deal with retardation effects in 
the exchanged pion propagator, i.e., its energy dependence, as well as with the energy 
dependence of the 7rN scattering amplitude in the vertex, from which the produced pion 
is emitted. 

4.1.1 Factorisation of the effective re-scattering operator 

In the previous Chapter, we made the connection to the usual representation of the 
pion re-scattering operator for non-relativistic calculations, following the work of Ref. |84j . 
We started from the covariant two-nucleon Feynman amplitudes including final and initial 
state interaction (FSI and ISI, respectively), shown in Fig. 14. lb ) and 14.1b )). 

To obtain the effective re-scattering operator, the negative energy contributions in 
the nucleon propagators (to be included in complete calculations) were neglected. By 
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The S-matrix approach for quantum mechanical calculations 



O 

a b) 

Figure 4.1: Feynman diagrams for pion re-scattering. The pion field is represented by a dashed 
line, the NN interaction by solid double line and the nucleons by solid lines. 
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V 
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Figure 4.2: The two time-ordered diagrams for FSI considered here. The additional stretched 
box diagrams are neglected. 



integrating subsequently over the energy of the exchanged pion the resulting Feynman 
amplitudes were transformed into those following from time-ordered perturbation theory. 

The irreducible "stretched box diagrams" (i.e., those with more than one meson in 
flight in the intermediate states) give a very small contribution and can be therefore also 
neglected (as we showed in the previous Chapter in Fig. 13.71 -Fig. 13.10)) . Thus, the full 
covariant amplitude is in the lowest order Born approximation well approximated by the 
product of the NN potential and the effective pion re-scattering operator, which can be 
extracted from these time-ordered diagrams (Fig. 14.2)) . Summarising, what was obtained 
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in Chapter El for the final state distortion amplitude, was 

M ™ = /tT^IFXp " ~-0 rs (4.1) 

J (2iry Fi + F 2 - oJi - u) 2 + is 



2u„ 



E 2 -u 2 - uj n Ei - lji - E n - u n 



\fG*. (4.2) 



The effective pion re-scattering operator was factorised into an effective pion re-scattering 
vertex / and an effective pion propagator G n . For the diagram with FSI (Fig. 14. 2j) this 
factorisation reads (Eq. flUD), Eq. (ETHH) and Eq. lETTSt ): 

/ = — [(E 1 -uj 1 -E 7r -uj 7r )f(uj 7T ) + (E 2 -uj 2 -uj 7T )f(-u 7T )} (4.3) 

G * = "TP W ^ V ^ 

[Ei - Ui - En - u 7T )(E 2 -u 2 - u) n ) 



2Un 



1 1 

+ 



oj 2 — uj a Fi — uoi — uj a 



(4.5) 



As before, q' is the momentum of the exchanged pion, Lu n = a/ vn\ + q 12 is its on-mass-shell 
energy and f(uj n ) is the product of the amplitude with the nNN vertex. The variables 
Ei, Ui, Fi are the on-mass-shell energies of the i-th nucleon in the initial, intermediate and 
final state, respectively, E n is the energy of the produced pion, Ei + E 2 = Fi + F 2 + E n . 

Again, we note that / (u n ) in fact depends on the three-momenta and energies of both 
(the exchanged and the produced) pions and on the nucleon spin. We indicate explicitly 
only the dependence on the exchanged pion energy, as the only one important for the 
considerations below. 

The inclusion of some pieces of the integrand of Eq. (j4.1j) into the propagator 
and of others in the modified vertex / is somewhat arbitrary. The appearance of the 
unusual effective propagator G n and the effective vertex / is the result of combining two 
time-ordered diagrams with different energy dependence into a single effective operator. 

The NN interaction is in Fig. 14.21 and in Eq. (|4.1jl simulated by a simple a-exchange 
potential. Though not realistic, this interaction suffices for model studies of approxi- 
mations employed in derivations of the effective pion re-scattering operator, as done in 
references Refs. jHS |HH] • Since some results do depend on the behaviour of the NN scat- 
tering wave function, in particular in the region of higher relative momenta, we perform 
our calculations (as in Ref. jSlj) also with V a replaced by a full iViV T- matrix, generated 
from the realistic Bonn B potential |62j. 
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The S-matrix approach for quantum mechanical calculations 



We note that the meson poles are not neglected in the integration over the energy 
Q'q of the exchanged pion, which generates Eq. (|4.1|) (see Sec. 13. lj) . A result similar to 
Eq. (|4.1|) can also be obtained for the amplitude with the initial state interaction (ISI). 



Mi si 
f 



d 3 q' 

(2?r) 3 ' CT E x + E 2 - ui - u 2 + ie 



Yr 



O, 



Ek + Fi — U)\ — UJ n F 2 — Ll>2 — LOjt 



\tG. 



[{F 2 -u 2 - w v )f(v v ) + (K + F 1 -u 1 - u v )f(-u v )] 



(E w + F 1 -u 1 - u n ) (F 2 -UJ2- 



2u a 



E 2 — u 2 — u a Ei — uj\ — uj a 



(4.6) 
(4.7) 
(4.8) 
(4.9) 

(4.10) 



The two amplitudes differ however in the contribution from the pion poles to the remaining 
integration over the three-momentum. For the amplitude with FSI there are no such 
poles. However, for the ISI case there are values of the exchanged pion three-momentum 
for which the propagator G n has poles. These poles have been considered in all our 
numerical calculations for the cross section (see Sec. 13.3)1 . As we will see in this Chapter, 
they are one of the main reasons for deviations between several approximations and the 
reference results calculated from Eqs. (|4. 1114.5)1 and Eqs. (|4. 6114 .10]) . 

It is worth mentioning that although the FSI and ISI diagrams graphically separate the 
A^V interaction and the pion re-scattering part (when the stretched boxes are neglected), 
they do not define a single effective operator (as a function of nucleon three-momenta and 
the energy of emitted pion) . Since in these time-ordered diagrams energy is not conserved 
at individual vertices, each of these diagrams defines a different off-energy shell extension 
of the pion re-scattering amplitude. This is an unpleasant feature, since one would have 
to make an analogous construction for diagrams with both FSI and ISI. Moreover, one 
would have to repeat the whole analysis for systems of more than two nucleons. Only 
after the on-shell approximation is made consistently (in the next subsection), the pion 
re-scattering parts of FSI and ISI diagrams coincide and one can identify them with a 
single effective re-scattering operator. 
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4.1.2 The S-matrix technique 

The S-matrix technique is a simple prescription to derive effective nuclear operators 
from the corresponding covariant Feynman diagrams [HE1 EZ1 EH1 EH]- F° r electromag- 
netic operators and also for NN and 3N potentials at energies below the first nucleonic 
inelasticities, the S-matrix approach analytically reproduces the results of more labori- 
ous constructions, based on time-ordered or non-relativistic diagram techniques. We will 
prove here that the same is numerically true for energies above pion production threshold. 

This approach is well defined and understood below the meson production threshold. 
As a simple tool it was employed also above the threshold in Refs. |H1 EZ|, but only 
to derive the Z-diagram operators. The two-nucleon effective operators are by definition 
identified with the diagrams describing the irreducible mechanism of the corresponding re- 
action. The only exception are the nucleon Born diagrams from which the iteration of the 
one-nucleon operator has to be subtracted. The operators of the nuclear electromagnetic 
and weak currents, as well as the pion absorption operators and nuclear potentials, are 
obtained by a straightforward non-relativistic reduction of the corresponding Feynman 
diagrams, in which the intermediate particles are off-mass-shell but energy is conserved 
at each vertex: therefore the derived effective operators are also defined on-energy-shell. 
The nuclear currents and other transition operators are defined to be consistent with a 
hermitian energy independent NN potential, which has the usual one boson exchange 
form employed in realistic models of the NN interaction, and can be used also in systems 
of more than two nucleons. 

The (j-exchange potential 

For the cx-exchange potential the S-matrix technique in the lowest order of non-relativistic 
reduction yields 

A 2 - (ml + q*) 

where A = A x = A 2 , A x = e[ — ei and A 2 = e 2 — e 2 , with e[ and 6j being the on-mass- 
shell energies of the i-th nucleon after and before the meson exchange, respectively. As 
pointed out above, this defines the potential only on-energy-shell (where actually A = 0). 
However, the Lippmann-Schwinger equation and even the first order Born approximation 
to the wave function require the potential off-energy-shell. 

67 



68 



The S-matrix approach for quantum mechanical calculations 



The S-matrix prescription takes the symmetric combination 

A = ^ + ^= Al + A2 . (4.12) 
2 2 2 v ; 

Most realistic NN potentials, namely those fitted to the data below pion threshold, in 
particular the Bonn B[H2j potential used in this work, are energy-independent and static in 
the nucleon CM frame and are therefore consistent with this construction. The extended 
S-matrix approach jSHl EH EH EH] given by Eq. (j4.11|) defines the most general off-energy- 
shell continuation of V a as a class of unitarily equivalent potentials parameterised by the 
"retardation parameter" v. 

We show in Fig. 14.31 the comparison of the cross section for pp — ► ppir° with the 
approximation of Eq. (j4.11j) and the static approximation (A = 0) for the potential, to 
the "exact" reference result of Eq. (|4.1jl and Eq. (|4.6jl . From Fig. 14 .3\ one concludes 
that the effect of the approximations for V a is small till Q ~ O.lm^. As expected, the 
deviations of the described approximations (dashed and dotted line, respectively) from 
the reference result (solid line) increase with energy, with a maximum deviation of about 
20%. 

The pion re-scattering operator 

For the pion re-scattering diagram in the FSI case (Eq. ()4.2jl ). according to the S- 
matrix approach one assumes energy conservation at each vertex. Then the exchanged 
pion is no longer on-mass-shell and we have to replace f{uJ n ) —> f{E 2 — uj 2 ) and f{—uo n ) — > 
f(E 2 — 002)'- in the first time-ordered diagram the virtual pion is entering the re-scattering 
vertex and in the second one it is emitted from this vertex (as defined on Fig. 14. 2j) . 
Therefore, one obtains 



Ors - O s rs = -±^-f^— , (4.13) 



/ (E 2 - uj 2 



(E 2 - uj 2 ) - ujI 

which is of the S-matrix formjOS]- The nucleon energies before and after pion emission 
are E 2 and u 2 , respectively (see Fig. 14.41) . 

For three-momentum values corresponding to the on-mass-shell conditions of the nu- 
cleons and the emitted pion, the characterisation ()4.13j) of the S-matrix technique suf- 
fices. However, whenever the operators are used in convolution integrals for the nucleon 
states distortion, the momentum variable is free. The on-mass shell condition has then 
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Figure 4.3: Comparison of the cross section for pp — > ppn calculated with the static and 
on-shell approximations for the cx-exchange potential (dashed and dotted lines, respectively) 
to the reference result of Eq. (14. 1J) and Eq. (14. 6 j) (solid line), as a function of the laboratory 
energy Ti ab . The values of the cross sections are scaled by requiring a reference result of 1 for 

T lab = 440MeV. 
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±l— 2 -=r 2 
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Figure 4.4: Illustration of the S-matrix prescription. 
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to be reinforced explicitly by demanding the conservation of energy u w = E 2 — lu 2 = 
— (Ei — iiJi — E n ). This is implemented by setting additionally 

n = + + (4.i4) 

2 2 

as it is illustrated on Fig. 14.41 The general form of the S-matrix prescription (HSj for the 
pion re-scattering operator becomes henceforth 

6 s = ® . (4.15) 

Note that the S-matrix approach defines a single effective operator (j4.15|) both for the 
FSI and ISI diagrams. In Eq. (|4.15j) . the energy of the exchanged pion is determined by 
applying consistently energy conservation in both vertices of each diagram and guarantee- 
ing that intermediate particles are on-mass shell. The prescription (|4.14|) is also consistent 
with the most realistic NN potentials (see Eq. (|4.12Jl ). 

4.2 Energy prescriptions for the exchanged pion 

The on-shell approximation (also labelled the U E — E'" approximation as intro- 
duced in Refs. jHDllHSIHS] actually coincides with the approximation defined by Eq. (|4.13j) . 
The re-scattering operator of Eq. (|4.1H|) can indeed be obtained directly from Eq. ()4.3j1 
and Eq. (|4.4|) by the substitutions following from the on-energy-shell prescription and the 
energy conservation in individual vertices Q = E 2 — u 2 = — (Ei — u>i — E n ) as follows: 

A 1 -2 (E 2 -u 2 )xf (E 2 -u 2 ) + 0xf(uj 2 - E 2 ) ~ 

rs ^ 2(E 2 -u 2 ) (E 2 -u 2 -uj 7r )(u J2 -E 2 -^) rs - { - ' 

We note that this form of O rs only coincides with Eq. ()4.15|) for on-mass shell nucleons, 
which is not the case in intermediate states, where the momentum integration variable 
does not in general satisfy the condition of Eq. (j4.14j) . 

In the work of Ref. |H1] the extra kinematic factors in Eq. (|4.H|) multiplying the function 
/ (u n ) were interpreted as form factors and kept unaltered, i.e., the substitution f((J n ) — > 
f(E 2 — lo 2 ), f(—Ld n ) — > f(E 2 — uj 2 ) was made only in G n and f(uj n ) of Eq. (|4.2|) . not in 
the kinematic factors included in the function /. 

In the following we will analyse in detail the S-matrix approach as well as the effects 
of other approximations frequently used for the energy of the exchanged pion, namely 
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the on- shell, fixed threshold-kinematics and static approximations. The fixed kinematics 
approximation, besides considering on-mass-shell particles and energy conservation at the 
vertices, further assumes that the energy of the emitted pion is m n , its threshold value 
(Q — > in Eq. (|4.15)l ). On the other hand, the static approximation assumes that no 

energy is exchanged (Jl — > in Eq. (|4.15ll ). 

Table 14.11 lists the different approximations considered here and in other works for the 
full re-scattering operator. The corresponding approximations for the pion propagator 
G n and the re-scattering vertex / are in Table 14.21 and Table 14 .3[ respectively. 
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Or 



TOPT 



on- shell 



fixed-kin. 



static 



S-matrix 



FSI 



ISI 



1 

2ui w 



f(E 2 ~^2) 



1 

2oj t 



/(^) _|_ /(-^) 



'( 



(£2-^2) -w 2 - 



/(o) 



(E 2 -^2) + (^l+g,r--El) 

2 



F2— ^2— O'tt Fi+Ejt— <^1— Wtt 



/(^2-^2) 



'( 



/(o) 



(ui2-F 2 )+{F 1 +E^-uj 1 ) 
2 



(E 2 -^2) + (^l+-g7T-gl) 

2 



- ui 



{w 1 -Eir-F 1 ) + (F 2 -W 2 ) 

2 



Table 4.1: Energy prescriptions for the exchanged pion in the full re-scattering operator. The 
first line is the reference result coming from TOPT of Eq. (14. 2 j) and Eq. (14.71) . The second, 
third and fourth lines are for the frequently used prescriptions for the energy of the exchanged 
pion (on-shell[57J |Hni E21 IH3. fixed-kinematics[2U EH ESI and static approximation, re- 
spectively). The last line presents the explicit expressions for the full production operator in 
the S-matrix approach of Eq. (I4.15jl . 

Note that the S-matrix approach defines a single effective operator: although they may look 
different because the labels for the energies in the FSI and ISI cases are different, the operators 
in the last line have the same dependence on the on-shell energies. 
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Table 4.2: Frequently used prescriptions (on-shell(HZI IHII1 E2 IHHJ , fixed-kinematics |3o1 l55l l5o1 
155] and static approximation) for the energy of the exchanged pion in the pion propagator. 
The first line is the reference result coming from TOPT of Eq. (14. 4|) and Eq. ()4.9|) . The last 
line presents the explicit expressions for G n within the S-matrix approach. 
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FSI 


ISI 


TOPT 


/ K) // + / (-u x ) fi 


fMfi + f{-u v )f* 


on-shell 


f(E 2 -u 2 )f{ + f(u 2 -E 2 )f{ 


f(u2-F 2 )fl + f(F 2 -u 2 )f* 


fixed-kin. 


f if) fi + f (-^) fi 




static 


f(0)f{ + f(0)fl 


/(o)/i+/(o)/j 



Table 4.3: The same of Table l4~21 but for the re-scattering vertex / alone. The first line is the 
reference result of Eq. fl4.3jt and Eq. (14. 8j) . The functions f[ and f l k (k = 1,2) refer to the 
multiplicative factors mentioned in Sec. 13.11 

As in Chapter El for numerical calculations we considered the NN — > (NN)ir transi- 
tion in partial waves 3 Po — > ( 1 So)s for the pp — > ppn° reaction. The amplitudes and cross 
sections are evaluated both with the simple interaction V a of Eq. (|4.5J) and Eq. ()4.10j) and 
with the Bonn B potential |62j. 

We tested the S-matrix prescription for the re-scattering operator (|4.15jl and the ap- 
proximations discussed in the previous section (see Table UJJJ . We include also the results 
for the on-shell, fixed threshold-kinematics, static and the S-matrix approximations for 
the effective pion propagator G n , as listed in Table l4~2l 

In Fig. 14.51 we show that the amplitudes with the S-matrix operator s (dotted line 
with +'s on the upper panels) are the closest to the reference result (solid line). Us- 
ing the same approach both for the operator and for V a increases slightly the gap from 
the reference result (line with x's versus solid line on the upper left panel). The fixed 
threshold- kinematics version of O rs , denoted as 0^ fc , works well for small values of the 
excess energy Q, but starts to deviate rapidly with increasing Q (dashed-dotted line in 
the upper panels of Fig. I4.5J1 . The on-shell approximation for the operator (dotted-line 
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in the upper panels), although deviates largely from the reference result when the NN 
interaction is described through u-exchange, is close to the reference result for the Bonn 
B potential (right panel) 1 . The static approximation for the re-scattering operator (O st ) 
overestimates significantly the amplitude (dashed versus solid lines in the upper panels). 
The same conclusions hold individually for the FSI and ISI amplitudes (Fig. 14.61 and 
Fig. 14.71 respectively). 

The investigation of how much of these features survive below pion production thresh- 
old is on Fig. 14.81 The fixed kinematics approximation (dashed-dotted line) is now found 
to disagree the most with the reference result. The deviation of the static approxima- 
tion (dashed line) from the reference result (solid line) is much smaller when compared 
to the situation above threshold of Fig. 14.51 This justifies the traditional approximation 
of static exchange below threshold. As before, the several approximations for the pion 
propagator yield results very close to each other (lower panels of Fig. 14. 8|) . Also, the S- 
matrix approach (dotted line with +'s in the upper panels) describes quite satisfactorily 
the reference amplitude. 



1 A more detailed analysis of these differences can be found in Sec. 14.2.11 
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Figure 4.5: Absolute values of the FSI + ISI amplitude as a function of the excess energy 
Q = 2E — 2M — E w (in units of m^). The right(left) panels correspond to the amplitudes 
with cr-exchange(Bonn B potential) for the NN interaction. The amplitudes are taken at the 
maximum pion momentum q™ ax , determined by Q. The upper panels correspond to approxi- 
mations for the whole operator O rs ; the lower panels to approximations for the pion propagator 
G w only. The solid line (O ex ) is the reference calculation. The dashed, dashed-dotted, dot- 
ted lines and dotted-line with +'s correspond to the static, fixed threshold-kinematics and 
on-shell approximations and the S-matrix approach, respectively. The corresponding operators 
are O st , O fk , O on and s , and G st , G fk , G on and G s . The dotted line with x's corresponds 
to both the pion production operator and cr-exchange potential in the S-matrix approach. All 
amplitudes were normalised by a factor defined by the maximum value of the reference result. 
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Figure 4.7: The same of Fig. 14.51 but for the ISI amplitude. 
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Figure 4.8: Absolute values of the FSI + ISI amplitude as a function of the symmetric of 
excess energy Q (in units of m n ), below pion production threshold. The meaning of the lines 
is the same of Fig. 14.51 
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4.2.1 Expansion of the effective pion propagator 

The lower panels on Fig. 14.51 show that all considered approximations taken only for 
the effective pion propagator do not differ much from each other, as already found on 
Ref. [El]. It means that the choices for the energy of the exchanged pion in the effective 
propagator alone are not very decisive (solid line versus dotted, dashed, dashed-dotted 
and dotted line with x's). We notice however that for cx-exchange there is a considerable 
deviation of all these approximations from the reference result. 

In order to understand this we considered the expansion of the effective pion propa- 
gator G n in Eq. f)4.4j) in terms of an "off-mass-shell" dimensionless parameter y: 

2E - E n - u x - u 2 

y = — p — > 4 - 17 

u 1 -u 2 + E w 

which measures the deviation of the total energy from the energy of the intermediate state 
with all three particles on-mass-shell|93j . As before, E = E\ = E 2 . This Taylor series 
expansion can give an insight on the small effect of retardation effects in the propagator, 
and it reads 



G n -- 



EV+^i-i^A 2 _ , .2 

2 / TV 



-2E + E w + u 1 + u 2 ) 



(4.18) 





^Tay 



where Gj) ay has the form of the usual Klein-Gordon propagator. 

We notice here that in the case of the ISI amplitude, the representation of the pion 
propagator G w by its Taylor series, the first term of which is G^iy, fails due to the presence 
of a pole in the propagator. 

Figure E31 compares the first four terms Gj) (i = 1, ...,4) of this expansion (dashed, 
dotted and dashed-dotted line and bullets, respectively) with the full effective propagator 
in Eq. (|4.4jl (solid line), as a function of the two-nucleon relative momentum for two 
different values of the excess energy Q. The convergence of the series demands at least 4 
terms. Besides, as expected, this convergence is momentum-dependent. 

We have also compared the first term of this expansion with the already considered 
on-shell, fixed threshold-kinematics and static approximations for the pion propagator. 
These results are shown in Fig. 14.101 We realise that all these approximations are very 
near to the 1st order term of the Taylor series. The corrections arising from higher order 
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Figure 4.9: Convergence of the Taylor expansion ()4.18p of the pion propagator G w in the 
FSI amplitude as a function of the two-nucleon relative momentum. The solid line is the full 
effective propagator of Eq. Q4.4I) . The dashed, dotted and dashed-dotted lines correspond to 
the first three terms of the Taylor expansion. The bullets represent the fourth term. Left 
panel: at threshold; right panel: above threshold at maximum pion momentum for an excess 
energy Q = O.Sm^. The bottom panels zoom into the region of low relative momentum 
{q k < 500MeV). 

terms in the expansion are negligible only for low momentum transfer, more precisely in 
the range qt < lOOMeV. 

The deviations of G st , G^ k and G on from the effective propagator G^ given by Eq. ()4.4j) 
cannot explain the relatively large deviations obtained on the bottom-left panel of Fig. I4.5I 
between considered approximations and the reference result. These deviations follow from 
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FSI (Q=0m ,q =0) FSI (Q=0.5m ,q =q max ) 




20 
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Figure 4.10: Approximations for G w in the FSI diagram and the first term of the Taylor series. 
Left and right panels have the same meaning as on Fig. 14.91 



the ISI contribution. 

When we compare the results for the final state interaction (lower panels of Fig. 14 .6|) 
with the results obtained for the initial state interaction (lower panels of Fig. 14.7)1 . we 
recognise however that the deviation between the approximate results and the exact one 
is much more significant for the ISI amplitude. The on-shell G on approximation (and 
also the fixed kinematics approximation) gives a larger amplitude. Its effect is much more 
pronounced than it is for the final-state distortion. This difference can be understood since 
the initial state interaction induces high off-shell energies in the intermediate nucleons. 
This enlarges the gap between the G on or G^ k and the G exact calculations, compared to 
what happens for the FSI case, where the nucleons emit the outcoming pion before their 
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Figure 4.11: Imaginary part of the ISI amplitude M.dwba of Eq. ()4.6jl as a function of 
the excess energy. The exact result (full line) and fixed kinematics (dashed line), on-shell 
(short-dashed line) and static (dotted line) approximations for the pion propagator are shown. 



interaction happens. Also, what is behind the approximations G on , G^ k or G st being 
worse representations of the exact amplitude for the ISI case than for the FSI case, is 
the physics related to the logarithmic singularities of the exact pion propagator. These 
singularities are present only for the ISI amplitude, but are absent in the case of the 
approximate forms for that propagator. Their effect is most important for the imaginary 
part of the amplitude, which is shown in Fig. 14.111 

Also, the weight of the ISI term depends on the NN interaction employed. It is 
comparable to the FSI term for V a (for which the deviations are large, as seen on the 
bottom- left panel of Fig. 14.5)1 . but it is less important for the full Bonn B potential (and 
therefore the corresponding deviations on the bottom-right panel of Fig. 14.51 are indeed 
much smaller). 

All the findings for the amplitudes manifest themselves also in the results for the cross 
section. We show in Fig. 14.121 the effects of the considered approximations on the cross 
section, first taking only the FSI contribution. On the left panel the amplitude includes 
V a for the NN interaction, and on the right panel the Bonn B T-matrix is used. The 
curves compare the reference result (solid line in all panels) with the S-matrix results 
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(upper panels) and their fixed threshold-kinematics version (lower panels). The S-matrix 
approach (dashed line) is the closest to the reference result (upper panels of Fig. I4.12|) . 

For the case of the A" A" interaction described by V a we also show the result following 
from the S-matrix prescription applied to the iVTV interaction (dotted line on left panels 
in Fig. 14. 12)1 . For the fixed threshold- kinematics versions (bottom panel) the deviations 
from the reference result increase pronouncedly with the excess energy Q, as expected. 
The approximations for the energy of the exchanged pion taken in the pion propagator 

and in the re-scattering vertex f(u n ), but not in the kinematic factors of Eq. (j4.3j) 
(see Table l4~3*)l . overestimate the cross section by a factor of 5 (solid line with bullets). 

Finally, we present on Fig. 14.131 the comparison between the approximated total cross 
sections with both FSI and ISI included. The approximation dictated by the S-matrix 
approach (dashed and dotted lines on the upper panels of Fig. 14.13)) is clearly seen as the 
best one. For the Bonn potential calculation, it practically coincides with the reference 
result. As shown in the previous section, this procedure amounts to extend the on- 
shell approximation, used in Ref. [EI] for G n and f{uj n ) alone, also to the multiplicative 
kinematic factors showing up in the operator / (see Eqs. (j4.2M.4jl ). 

To conclude we notice, moreover, that for the realistic A^A" interaction, the difference 
between the S-matrix approach (upper right panel on Fig. 14.13)1 and its fixed threshold- 
kinematics version (lower right panel of the same figure) is not very important near 
threshold, provided the excess energy does not exceed ~ 30MeV (Q/m n ~ 0.2). 

4.3 Conclusions 

This Chapter discussed the choices for the exchanged pion energy which are unavoid- 
able in the three-dimensional non-relativistic formalism underlying DWBA. The main 
conclusions are: 

1) The usual approximations to the effective pion propagator E3 E31 EDI EI], 
obtained from a quantum-mechanical reduction of the Feynman diagram describing the 
pion re-scattering process, are rather close to the first order term of a Taylor series in a 
parameter measuring off-mass-shell effects in the intermediate states. The series converges 
rapidly for the FSI amplitude near threshold. As a consequence, retardation effects are not 
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Figure 4.12: Effects of the approximations for the re-scattering operator O rs and for the effec- 
tive pion propagator G w as a function of the excess energy Q. The cross section curves shown 
correspond to the FSI amplitude alone. The upper panels correspond to the re-scattering 
operator O rs given by Eq. ()4.15|) and the lower panels to fixed threshold-kinematics approx- 
imation. The solid line is the reference calculation (|4.1j) . The dashed line is the S-matrix 
calculation for the re-scattering operator O rs given by Eq. ()4.15j) (upper panels) and the fixed 
threshold-kinematics approximation for ()4.15|) (lower panels). The dotted line corresponds to 
take the S-matrix approximation not only for O rs , but also for the a-exchange interaction 
()4.1ip . The dashed-doted line corresponds to the on-shell (upper panels) and fixed threshold- 
kinematics (lower panels) prescriptions only for G n . The solid lines with bullets refer to the 
energy prescriptions taken for the propagator G w and for /(o» T ) in Eq. (14. 3|) . as in [84J, but not 
for extra kinematic factors in /. All the cross sections were normalised with a factor defined 
by the maximum value of the reference result. 
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Figure 4.13: The same of Fig. 14.121 but for the total (FSI+ISI) cross-section and considering 
only the approximations for O and G n . 



decisive in the pion re-scattering mechanism near the threshold energy for pion production. 

2) As for the pion energy in the nN re-scattering amplitude, the on-shell approach 
when used only in f(u n ) overestimates significantly the reference result. 

3) Nevertheless, and this is the key point of this Chapter, this deviation is dramatically 
reduced if the approximation coming from the S-matrix approach is used consistently in 
the whole effective operator (HB]- This procedure amounts to extend the on-shell approx- 
imation used in Ref. [EI] for G n and f(u) n ) to the full operator O rs , including kinematic 
factors which differently weight the two dominant time-ordered diagrams. The amplitudes 
and cross sections obtained with the S-matrix effective operator are very close to those 
obtained with the time-ordered one in the considered kinematic region. Nevertheless, the 
static approximation works well below threshold. 
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The re-scattering operator for the neutral pion production in the isoscalar 7rN channel 
indeed seems to be relatively unimportant: its enhancement reported in previous papers 
followed from inconsistent or too crude (static or fixed threshold-kinematics) treatment 
of the energy dependence of the effective operator. Our findings explain why the 
calculation of Ref . [HO] , where the on-shell approximation is used, artificially enhances the 
contribution of the isoscalar re-scattering term. On the other hand, our results indicate 
that the fixed kinematics choice done in Refs. |3T)1 EH1 EH] for the different production 
operators considered, is valid in the restricted kinematic region where Q < 0.2m,,-. 

The re-scattering mechanism is filtered differently by other spin/isospin channels in 
pion production reactions. For charged pion production reactions the general irreducible 
re-scattering operator comprises also the dominant isovector Weinberg- Tomozawa term of 
the 7rN amplitude, and its importance is therefore enhanced. Next Chapter investigates 
these reactions within the S-matrix approach. 
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to charged and neutral pion production reactions. The contributions from the direct- 
production, pion re-scattering, Z-diagrams and explicit A-isobar excitation terms are 
considered. High angular momentum partial wave channels are included and the conver- 
gence of the partial wave series is investigated. 
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5.1 Pion production operators 

Chapter HJ following the work of Ref. has shown that the S-matrix approach re- 
produces the energy dependent re-scattering operator for neutral pion production derived 
from time-ordered diagram techniques. This conclusion provides the motivation to apply 
this approach further to total cross sections for not only neutral, but also charged pion 
production. As discussed before, the S-matrix technique has the advantage to define a sin- 
gle effective operator, contrarily to what happens when deriving the quantum-mechanical 
production operator in the framework of non-relativistic, time-ordered perturbation the- 
ory. 

Another usual problem of pion production calculations concerns the need of a realistic 
nucleon-nucleon potential valid for the energies above the pion production threshold, 
which necessarily are needed for the initial nucleonic states. In this Chapter, we apply for 
the first time an NN interaction including two boson exchange potentials which describes 
well not only NN phase-shifts but also inelasticities. 

5.1.1 The mechanisms and their operators 

We considered the contributions of the impulse, re-scattering, A-isobar mechanisms 
and Z-diagrams as represented in Fig. 15. IK 

The ttN amplitudes are those of chiral perturbation theory from Ref. [HE], resulting 
from the leading-order Lagrangian of Eq. (|2.11j) and the next-to-leading order Lagrangian 



where pi and q\ are the nucleon momentum before and after pion emission, is the 
momentum of the emitted pion and a is an isospin index. The second term of Eq. (|5.1j) 
accounts for the nucleon recoil effect. For the re-scattering ttN-ttN vertex (with no isospin 

^^Our calculations do not include contributions to the diagrams on Fig. 15.11 involving the two T- 
matriccs. Wc note however that all the contributions with at least one T-matrix are considered, together 
with the undistorted operator. 
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a) b) c) d) 

Figure 5.1: Mechanisms considered for charged and neutral pion production, a) direct pro- 
duction, b) re-scattering, c) Z-diagrams and d) A-isobar contribution. T is calculated from 
realistic NN interactions. 
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where Q' is the zeroth component of the exchanged pion four-momentum and e a b c the 
Levi-Civita tensor (as before, a, b and c are isospin indices). 

The A-isobar contribution was calculated from Eq. (j2.11|) and reads 
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A - E„ tt 2 



(5.4) 



where Q is given by Eq. (j4.14|) . Since the A-resonance is here included explicitly, its 
contribution needs to be subtracted from the values of the low energy constants q's of 
Eq. (J5.2j) . This was done taking the heavy baryon limit of Eq. ()5.4j) . giving the contribution 
cf ~ — = — 2.78GeV _1 which was subtracted from the — 5.29GeV _1 value of the c 3 
constant appearing in the (q n ■ q') term of Eq. (15.21) . 

The heavy-meson mediated Z-diagrams contributions given in Ref. [UJ and used in 
Ref. [27] for the 7r° production calculation, were also included: 
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for cr-exchange, and 



M 




(5.6) 



for oj exchange, where V§ and Vy are respectively, the scalar component and isospin 
independent vector component of the nucleon-nucleon interaction. The i-th nucleon final 
momentum is qi, and k = q\ — f>2- The kinematics conventions for the diagrams may be 
found in Appendix El Details for the partial wave decomposition of the amplitudes are 
given in Appendix |E| 

5.1.2 Nucleon-nucleon potentials at intermediate energies 

Since the energies considered for the nucleonic initial state interaction are necessarily 
above the pion production threshold, we took a recently developed NN interaction fitted 
to the nucleon-nucleon scattering data|9() [ I95j. This nucleon-nucleon interaction was 
developed by the Ohio group as an extension of the Bonn family [B^j of realistic nucleon- 
nucleon potentials well into the intermediate and high energies region. 

The free NN T-matrices are obtained from a meson-exchange model including A- 
isobar and N*(1440)/N*(1535) degrees of freedom to account for pion-production above 
the pion thresholdPHllIIIJIHIi]- To properly generate inelasticities, the model incorporates 
the particle data book nucleon resonances as intermediate excitations within two meson 
exchanges loops. As an extension of the family of the Bonn potentials, it describes well 
the NN phase-shifts and inelasticities up to 1 GeV. 

The model is based on the solution of the relativistic Thomson equation. This equa- 
tion was favoured to time-ordered field theory since it gives a good description of the 
nucleon inelasticities, without the numerical complexities from the one-pion exchange cut 
singularities. In time-ordered perturbation theory the last ones are formally unavoidable, 
but actually, the main contribution to the nucleon singularities arise from the nucleon 
resonances, and the effect from the one-pion exchange cut is indeed quite small. 

The Thomson equation for NN scattering was solved in order to obtain the fit of 
the NN interaction to both phase-shifts and inelasticities. Table l5~Tl lists the parameters 
employed by the Ohio group NN potential and consistently used in all the calculations 
in this work. Fig. 15.21 shows the quality of the phase shifts in the relevant energy region 
for the initial and final nucleonic states. 
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Figure 5.2: Phase shifts calculated with the Ohio group NN potential. 
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Table 5.1: iViV-meson and iVA-meson vertices of the Ohio group NN model. The meson 
mass is m, the coupling constants are g and /, and A is the form-factor cutoff. All of the 
corresponding form factors are of the monopole type y j^g^ J |53|- 

5.2 Calculational details 

5.2.1 Three-body kinematics 

The phase-space of the reaction NN — > B1B2X is five- dimensional: the three parti- 
cles in the final state introduce 3x3 = 9 degrees of freedom, but the four-momentum 
conservation reduces this number to 5. For energies close to threshold, the final state 
can be treated non-relativistically. The natural coordinate system is therefore given by 
the Jacobi coordinates (HHl where first the relative momentum of one pair of particles is 
constructed, and then the momentum of the third particle is calculated relatively to the 
two-body system. Obviously, there are three equivalent sets of variables possible. 

For reactions of the type NN — > BiB 2 x it is common to work with the relative 
momentum of the two-nucleon system and to treat the meson x separately, as the spectator 
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particle. This choice is theoretically the most convenient, since one is working already with 
the relative momentum of the dominant final state interaction. Close to threshold, for the 
final-state distortion (but not for the initial-state distortion), the nucleons can be treated 
as non-relativistic, and thus for any given relative energy of the outgoing two-nucleon 

system, 

The modulus of the meson momentum \q x \ is fixed by energy conservation. The phase- 
space is then 

t = {e,n k ,n n } (5.8) 

where and Q n are the solid angle for the relative NN momentum and for the pion 
momentum, respectively (see Appendix [HJ)- initial state will be denoted by p. Explicit 
expressions for the vectors appearing are given in Appendix [XJ As long as the initial 
state is unpolarised, the system has azimuthal rotation symmetry, reducing the number 
of degrees of freedom from 5 to 4. 

5.2.2 The role of the different isospin channels 

Meson production reactions in nucleon-nucleon collisions give also the possibility to 
manipulate another internal degree of freedom: the isospin. Nucleons are isospin | par- 
ticles in a SU (2) doublet, 



\p) = 



t = 2>*3 = +2 ) 



t = ^3 = -^), (5.9) 



T 2 

where \p) and \n) represent the proton and neutron wave functions, respectively, and 
is the third component of the isospin operator t. 

A two-nucleon pair can be either in an isotriplet state (total isospin T = 1 with the 
three possible projections of the total isospin T 3 = +1 for a pp state, T 3 = for a pn state 
and T3 = —1 for a nn state) or in an isosinglet state (T = with T3 = for a pn state). 

Even in an isospin conserving situation, various transitions are possible, depending on 
whether an isovector or an isoscalar particle is produced. In particular, pion production 
in NN collisions can be expressed in terms of only three independent reactions with 
cross sections denoted by cr 10 , er 01 and er n (here, the reference of Ref. jH] is adopted, 
where Tj and Tf in o^t, are the NN pair initial and final total isospin, respectively). 
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Isospin conservation forbids the existence of <7oo- 111 total, there are seven relations for 
NN -> NNtt reactions 0Ej: 



a (pp — > pj»7T ) = an (5.10) 

a (pp ^ pnir + ) = (7io + 0ii (5.11) 

a (rap -> 72£>7T°) = - [010 + o- i] (5.12) 

(np -> p^TT") = - [011 + O i] (5.13) 

(np -> rari7r + ) = - [011 + O i] (5.14) 

(nn — > npir~) = a w + o"n (5.15) 

(nn — > nn7r°) = <Tn (5.16) 



Note that none of the reactions involves more than two of the fundamental cross sections 
and that the last two reactions are related with the first two by charge symmetry. The 
relative strength of the different transition amplitudes can thus provide significant infor- 
mation on the production operator. At medium energies, the cross sections an and 010 
are dominated by the excitation of the intermediate A33 resonance and are well measured 
even close to threshold. In contrast, the isoscalar cross section 001, which has to be ex- 
tracted from pion production data in neutron-proton and proton-proton collisions is not 
well known. Actually, the determination of 001 includes a principal model dependence, as 

Cnp->ppiT- — 2 + °oi) (5-17) 

only holds in the case of exact isospin invariance. However, due to the different particle 
masses entering in the reactions np — > pprr' , np — > nnir + and pp — > ppir , isospin invari- 
ance is only an approximate symmetry, and thus the comparison of the cross sections can 
not be performed at the same beam energy, as mentioned before 

Meson production data close to threshold in np collisions suffer from the fact that 
only relative cross section measurements have been performed, as a consequence of the 
(poor) quality of the neutron beam. Unfortunately, the isospin cross section 001 can only 
be accessed via these np interactions. It is obtained from the reaction cross section for 
np — > ppn~ by subtracting the ern part, which results in larger uncertainties than in other 
isospin cross sections. 
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5.2.3 Partial wave analysis 

Some of the complexity inherent to three-body final state reactions can be dealt using 
partial wave analysis. For energies close to threshold, since the total kinetic energy in the 
final state is small, the reactions are dominated by transitions to states with low orbital 
angular momentum values. Therefore, an expansion of the outgoing spherical waves will 
involve mainly low-order spherical harmonics, and thus the observables should have a 
relatively simple angular dependence. 

The amplitude of the outgoing wave in spin state o\OiOx obtained when the nucleons 
in the initial state have spin projections a a and cr& j^Hl EH] is given by 

M™" = iKi \/jJTT ( s « s >>Wb \ s i°i) (LiSiOvi \JM) (5.18) 

OtiOtfam 

(siS 2 cri(T2 \Sf(jf) (SfLifCTfXf \jfrrif) (sj x cr x \ x \j x m x ) 
{t a t h v a v h \Tih>i) (txhvxvz \TfVf) (T f t x v f v x \T{Ui) 
Uf.l,i»j<», \JM) C n n; )t ('//,) y.v (?) 

where Li, Si and J are respectively, the two initial nucleon orbital angular momentum, 
spin and total angular momentum, and U aia} is a matrix element depending on the 
quantum numbers in both initial state a, = {Lj, Si, J and T{\ and final state atf = 
{Sf, Lf, l x , jf, j x and Tj}. This matrix element is also a function of the energy sharing 
parameter e. 

The individual nucleon spin quantum numbers are s a and s&, with the corresponding 
magnetic quantum numbers a a and a^. In the final state the nucleon spin quantum 
numbers are represented by S\ and s 2 with spin projections <j\ and oi, and the iViV total 
spin quantum number is Sf. The AW relative orbital angular momentum is Lf and the 
iViV total angular momentum is jf. 

For the produced meson, the spin is represented by quantum numbers s x and a x , the 
angular momentum relative to the iViV centre of mass is l x and the vector sum of s x and l x 
is j x . The isospin quantum numbers are represented by t and v. Ki is a kinematic constant 
arising from phase-space (as defined in Eq. (|G.12JI ). The sum in (|5.18j) extends over the 
initial and final state quantum numbers a.\ and a.2, and also the magnetic quantum 
numbers a m = {a iy Of, M, nif, m x , A/, X x , z/j and vf\. These notations are summarised 
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Spin s a and s& Si s± and S2 5/ s x 

projection a a and cr;, <7; <7i and <r 2 a x 

Orbital mom. Lj Lf l x 

projection A/- A x 

Total ang. mom. J jf j x 

projection M rrif m x 

Isospin t a and t& 7$ ti and t 2 Tf t x 

projection v a and v\> Vi V\ and v 2 Vf v x 



Table 5.2: Summary of the notations for the nucleon and pion spin, orbital momentum, total 
angular momentum and isospin. 

in Table O 

Note that as in the two particle case, the complexity of the angular pattern is deter- 
mined by the final state orbital angular momentum values Lf and l x . Since s x = a x = 
=>- j x = l x and m x = \ x , ()5.18|) simplifies to 

Ol«2am 

(siS 2 CTi(T2 \Sf<Tf) (SfLfCTfXf \jf1Tlf) (jfl x m f X x \ JM) 

U aia2 Y x L ; (q k )Yl : (q) 

with 

In = {t a t h v a v h \TiVi) (tit 2 viv 2 \T f Uf) (Tft x u f u x \TiUi) (5.20) 

If one follows the obvious procedure of including only initial and final angular momentum 
states that are consistent with the anti-symmetrisation requirements, then the resulting 
formula will satisfy all the necessary symmetries for identical particles in the initial and 
final states. The only requirement is that when the initial state has identical particles, 
an extra factor of 2 must be included in the cross section to reproduce the appropriate 
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incident flux normalisation [98| 2 . This extra factor of 2 is already included in Eq. ()5.19jl 
via the definition 1)5.21)1 of I s . 



for 7T° 



V2 

\f2 for 7r + production 
1 for 7T~ 

Expansion of the cross section 

The simplest observable is the unpolarised cross section, which reads jHHl EH]: 

1 



(5.21) 



(5.22) 



(2s a + i; 

To obtain a formula for the partial-wave expansion of the cross section which shows ex- 
plicitly the various allowed angular dependencies, one substitutes Eq. (J5.19)) in Eq. ()5.22|) 
and then combine the spherical harmonics and carry out an angular momentum reduction. 
More precisely, using (jTHfl} . (pT88|) . (frOsjl and (jHT2j), one obtains, 



\K, 



(2s a + 1) (2s b 
with the coefficients! 



T) (h) ^ E <L qL U*U:,y™ Lq g) (5.23) 

/ \ / r r. . r. 



= (W (-l) Lf+l '* +L '< 6 SiSi 5 Sf s>V2Li + 1 V / 2X[+1 V / 2T7TT V /2L' / + 1 (5.24) 



and 



+ ly/2l' x + lv/2j f + ly^j) + lv / 27TTv / 2J / + 1 
(L/L^OOlLpO) (^^00|L ? 0) (LiL-00|L0) 

jf lx J 

W (j f S f L p L' f ;L f f f )W {L^W; JL' % ) \ f f \> x J' 

L p L q L, 



for 7T C 



(Us 



\(\\vxv 2 IT/0) (7)101 1 1 1 >] 2 for 7T+ production 
^ | [<||^6 1^0) (111 - 1 iTfi)] 2 for 7T- 



(5.25) 



2 The corresponding requirement for the particles in the final state is included in Eq. 1)5.19(1 through 
the definition of K, (Eq. {£33). 
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(Usf 7T° 7T+ 71- 

T=l 1 11 
2 4 

Table 5.3: (I n I s ) 2 values for the several pion production reactions considered. 

In Eq. (I5.23j) and Eq. (|5.24jl . a = {Li, Si, J, Sf, Lf, l x , jf, Oi, Of, M,rrif, A/, \ x } is a short- 
hand for summation over all the indices. The values for (I n I s ) 2 are listed in Table 15.31 

The various cross section terms in ()5.23j) arise from the interference between the two 
matrix elements depending on the quantum numbers a and a'. Moreover, from Eq. (J5.24D 
one sees that two partial waves may interfere only if the NN spin quantum numbers 
match in both the initial and final states. Since the 8s f s' f is & l so present in the general 
formula for the polarisation observables, one can divide the partial waves in groups, where 
interference is only permitted within each group. For instance, for tt° production, these 
first two groups are the Sf = group and the Sf — 1 group, respectively {Ss, Sd, Ds} 
and {Ps, Pp}. Since Ss is dominant at energies close to threshold, one concludes that the 
Sd and Ds terms may be more important at low energies that one might otherwise have 
thought. Note also that the expansion of the cross section involves only angular functions 
yipL q w hh A = 0. This results from the angular momentum algebra and is a consequence 
of the fact that the cross section must be invariant under rotations about the z-axis. 

5.2.4 Selection rules for NN -> NNx 

The selection rules for iVW — > NNx reactions are based on the symmetries of strong 
interaction that imply the conservation of parity, total angular momentum and isospin. 

For the two-nucleon system, its parity (in other words, the behaviour of its wave 
function under a reflection of the coordinate system through the origin), determined solely 
by the relative angular momentum of the two particles, is (— 1) L . Also, the symmetry of 
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the spin and isospin wave function is (—1) and (—1) . Then, because the nucleons 
are fermions, 

(-1) = (_ 1 )5 l +T l +L I = ^Sf+Tf+Lf ^ 2 Q) 

where Lj (Lf), Si (Sf) and (Tf) denote the angular momentum, total spin and total 
isospin of the initial (final) two-nucleon system. In addition, for a reaction of the type 
NN — > NNx, one finds from parity conservation [H] 

(-1)^ = 7 r a .(_i)( L /+ 1 -) (5.27) 

where ir x and l x denote, respectively, the intrinsic parity and the angular momentum (with 
respect to the outgoing two-nucleon system) of the particle x. The two criteria (|5.26|) and 
()5.27j) can now be combined to give 

(-1) (A5+Ar) =7r a! (-1) 1 " (5.28) 

where AS (AT) is the change in total (iso)spin when going from the initial to the final 
NN system. We can now analyse in detail the channels contributing to the different 
charge reactions: 

• pp — > ppTT 

The pp — > ppn° reaction corresponds to a (T^ N — 1J — > (Tf N = l) transition. The 
possible partial waves for pp — > ppn result from applying Eq. (|5.27jl and Eq. (J5.28|) . 
They are listed in Table up to J = 3, and restricted to pion s- and p-wave states 
relatively to the final iV./V" system. The channels are grouped in sets of equal J. 
Within each set, the different channels were ordered according to the total angular 
momentum jf of the final iViV pair. 

• pp — > pmr + 

The partial waves considered for pp — > pmr + are both in Table 15.41 and Table 15. 5^ 
since pn can be either in a isospin isotriplet state T NN = 1 (Table 15. 4j) or in an 
isosinglet state T NN = (Table l5~5|) . The selection rules for this last case follow 
from that (Lj + Si) must be even and (Lf + Sf) must be odd, since T^ N = 1 and 
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Table 5.4: The lowest partial waves for (NN T= i) — > (NN T= i) ir, ordered by increasing values 
of J. The notation is Si(Sf), Li(Lf) and J(jf) for the total spin, the orbital momentum and 
the total angular momentum of the initial(final) two-nucleon pair, respectively. The orbital 
angular momentum of the pion relative to the final iViV pair is l x . The different channels are 
grouped in sets with the same total angular momentum J. Within each set, the order of the 
channels is ruled by the total angular momentum jf of the final NN pair. 
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Table 5.5: The same of Table El but (NN T=1 ) -> (NN T=0 ) tt. 
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• pn — > ppn 

The lowest partial waves for pn — > ppir~ are in Table 15.41 and Table 15.61 which 
correspond to the case when pn is in an isotriplet and isosinglet state, respectively. 
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Table 5.6: The same of Table El but (NN T=Q ) -> (NN T=1 ) ix. 

It is clear from Table 15.41 and Table 15.61 that partial waves with Ss final states can 
only contribute to a"n and partial waves with Sp final states only to cr i- 

Reactions proceeding dominantly via the A resonance are expected to be significantly 
larger than others where the excitation of the A resonance is suppressed by selection 
rules|2]. Since the A has isospin | and the nucleon has isospin |, the iVA system can 
couple to total isospin equal to 2 or 1. On the other hand, the nucleon-nucleon system 
can couple to total isospin equal to 1 or 0. Thus, the iVA states do not contribute to the 
nucleon-nucleon channels with isospin equal to and only the iVA channels with isospin 
equal to 1 will couple to the nucleon-nucleon systemjH]. Table IB~7I lists the lowest partial 
wave states for the iVA system. For instance, for o"oi, intervening in 7r _ production, A- 
excitation in the initial state is impossible because a pion-nucleon state with isospin | 
cannot couple with a nucleon with isospin zero. A A excitation is possible in the final 
state but then the intermediate AiV system needs to be in a L& N = 2 state (see Table I5~7jl . 
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Table 5.7: The lowest iVA partial wave states(21 IE]- 



The necessary rotational energy in the A A" system can be estimated through 

Lan (Lan + 1) 



(5.29) 



2/xr 2 

with /i the reduced mass. Assuming a distance r of lfm between the A and the nucleon, 
one obtains 

73MeV for L AN = 1 (5.30) 
219MeV for L AN = 2, (5.31) 

which means that the excitation of the A-resonance will be strongly suppressed at thresh- 
old in this particular isospin channel. The cross section an can proceed through the 
A but here the AA" system, like the initial A^A" system, has at least an orbital angular 
momentum of 1. For the cross section o"io, and consequently for the n + production reac- 
tion, however, there is one partial wave where the AA" system is in an s-wave state and 
therefore, the A excitation contribution is expected to be important given the weight of 
the ( 1 D 2 ) NN -> ( 5 S 2 ) NA ( 3 S 1 )p transition^. 

So far, most of the investigations considered only the lowest-partial waves in the 
outgoing channel with the pion in an s-wave state relative to the nucleon pair. Such 
calculations permit only conclusions on the absolute magnitude of the production cross 
section near threshold. The inclusion of higher partial waves, in the A^A" as well as 
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in the irN sector, allows predictions for differential cross sections and, in particular, 
spin dependent observables. Therefore, it is possible to examine whether the considered 
production mechanisms lead to the proper onset of higher partial waves, as suggested by 
the data. Results of the latter kind are especially interesting because they reflect the spin- 
dependence of the production processes and thus should be very useful in discriminating 
between different mechanisms |6()j. 

Very close to threshold, n° production is naturally dominated by s-wave production. 
To go beyond threshold, we made convergence studies of the partial wave series by in- 
cluding partial waves up to a total angular momentum of the system J = 3. In our 
calculations, the coupled channels were properly taken into account in the convergence 
studies. 

5.3 Convergence of the partial wave series 

In Fig. 15.31 we show the results for the convergence of the partial wave series for the 
cross section of 7r° production (with all the mechanisms discussed in Sec. 15.11 included). 
As expected, the importance of the contributions from the channels with higher J's in- 
creases with the laboratory energy, 7] a &. Till ~ 30MeV above threshold, the contributions 
from the several sets of J's are immaterial compared to the J = result, which alone 
describes the data in that energy region. The same conclusions hold when the NN in- 
teraction is described through the Bonn B potential [H2|> due to a smooth variation of the 
nucleon phase-shifts in that region. We compare in Fig. 15.41 the cross sections with both 
interactions. 

The cross sections with the Ohio NN model for the nucleon-nucleon interaction (dot- 
ted, dashed-dotted, dashed and solid line with +'s) are slightly lower than the ones with 
the Bonn B potential (dotted, dashed-dotted, dashed and solid lines). Also, naturally, 
the deviations between the results obtained with the two interactions increase with the 
energy, and also with J. 

For 7r~ production, the deviation between the cross section with all the contributions 
up to J = 1 (dashed-dotted line in Fig. I5.5J) . J = 2 (dashed line), and up to J = 3 (solid 
line) is very small. However, unlikely to what happens for n°, the J = channels (dotted 
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Figure 5.3: Convergence of the partial wave series for pp — > ppir . The dotted, dashed-dotted, 
dashed and solid line correspond, respectively, to the cross sections with all the contributions 
up to J = 0, J = 1, J = 2 and J = 3, as listed in Table l5.4l The nucleon-nucleon interaction 
is described through the Ohio NN model. The data points are from Ref. |23j. 



line) are not enough to describe the data in the region of 30MeV above threshold. The 
same is found when the NN interaction is described using the Bonn B potential, as it is 
shown in Fig. 15.61 

We show in Fig. 15. 71 the convergence of the partial wave series for the ir + production. 
For this reaction, although the convergence is slower, the experimental data near threshold 
are also described. As already found for 7r~ production, since the channels with J = 
(dotted line) are not enough to describe the data, one needs to include at least the channels 
up to J = 2 (dashed line). 
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pp -> ppjt (J=0) 



pp -> pprc°(J<1) 
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Figure 5.4: Effect of the nucleon-nucleon interaction in the cross section for pp — > ppir . 
The dotted, dashed-dotted, dashed and solid line correspond, respectively, to the cross section 
considering the Bonn B potential for the nucleon-nucleon interaction and all the contributions 
up to J = 0, J = 1, J = 2 and J = 3. The corresponding lines with +'s consider the 
nucleon-nucleon interaction to be described by the Ohio NN model. The data points are 
from Ref. [23]. 



Moreover, the NN potential used to describe the NN final- and initial-state interaction 
can be decisive for the case of 7r + production. We compare in Fig. l5.8l the calculated cross 
sections with the Bonn B and Ohio interactions, for all sets of J's considered. As expected, 
the deviation of both results increase with Ti a b and with J, but for this reaction the Bonn 
B potential fails in describing the data in a much wider region than the Ohio NN model 
(solid line vs. solid line with +'s in the bottom-right panel of Fig. 15 .8j) . 
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pn -> pprc (J=0) 
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Figure 5.6: Cross section for pn — > ppn for the Ohio NN model and Bonn B potential for 
the NN interaction. The meaning of the lines is the same of Fig. I5.4I The data points are 
from Ref. 
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Convergence of PW series for pp -> prm' 
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Figure 5.7: Convergence of the partial wave series for pp — > pnn + . The meaning of the lines 
is the same of Fig. I5.3I The data points are from Ref. |28J. 
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pp -> prm + (J=0) 



pp -> prm + (J < 1) 
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Figure 5.8: Cross section for pp — > pnir + for the Ohio NN model and Bonn B potential for 
the NN interaction. The data points are from Ref. [28J. 
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Convergence of the PW series for pp -> ppji 
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Figure 5.9: Comparison between the contributions for ir° production of the mechanisms of 
Fig. I5.ll The dotted line corresponds to the direct-production mechanism alone (diagram a)) 
and the dashed-dotted line considers the additional contribution of the re-scattering mecha- 
nism. In the legend, A# refers to the subtraction of the infinitely heavy limit of Eq. ()5.4|) 
from the q values of Eq. (I5.2jl . The dashed line further includes the explicit A-resonance 
excitation. The solid line represents all the contributions (diagrams a)-d)). The data points 
are from Ref. 



5.4 The role of the different production mechanisms 



We compare in Fig. I5.9I the contributions of the pion production mechanisms that 
are illustrated on Fig. I5.ll The direct-production and re-scattering mechanisms alone 
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(dashed-dotted line), highly suppressed, are clearly not enough to describe the ir° pro- 
duction data. Also, the additional mechanism considering the explicit excitation of a A 
(dashed line) was found to be small (this is natural, given the discussion in Sec. I5.2.4|) . 
As already discussed on Ref. [SZJ, t ne Z-diagrams (solid line) are decisive to remove the 
discrepancy with the experimental data. 

For 7r~ production, the re-scattering operator comprises also the Weinberg- Tomozawa 
term, which does not contribute to 7r° production, and its importance relative to the 
single nucleon emission term is therefore enhanced (dotted line vs. dashed-dotted line in 
Fig. I5.1()j) . For this case, the A-isobar mechanism is also found to play an important role 
(dashed line in Fig. I5.10j) . 

The relative importance of the several mechanisms of Fig. 15. 1| for ir + production, is 
shown in Fig. 15. 1 1| where the re-scattering mechanism is found to contribute significantly 
(dashed-dotted line). 

Since the A-resonance is included explicitly in all our calculations, its contribution was 
subtracted from the values of the low energy constants q's of the irN scattering amplitude 
in Eq. (|5.2|) . which was done taking the infinitely heavy or static limit of Eq. (|5.4|) (this was 
indicated by the superscript A H in the legends of Fig. 15.91 Fig. 15.101 and Fig. 15. lip . We 
show in Fig. 15.121 the effect of not including explicitly the mechanism of pion production 
through the A-excitation. In that figure, in contrast to previous ones, the re-scattering 
amplitude (in the dotted line) is taken with the low energy constant C3 including the A 
contribution. The solid lines are the full calculation for ir°, tt~ and 7r + . They coincide 
with the full lines in Fig. 15. 9| Fig. 15.101 and Fig. 15. 1H respectively. 

For all the cases, the inclusion of the A improves the description of the experimental 
data. As expected, its importance increases with the energy. Also, both for 7r° and ir + 
production, the inclusion of the A increases the cross section. These general trends were 
already found for the Julich modeljHniinniEI], presented in Sec. 12. 21 (see Fig. l2.8p . However, 
contrarily to this phenomenological model, in which all the short range mechanisms are 
included through cj-exchange (diagram c) of Fig. 15. ip and adjusted to reproduce the 
total cross section of the reaction pp — > ppn° close to threshold, in our calculations no 
adjustment is made, as the parameters for the Z-diagrams are taken consistently from 
the NN interaction employed. 
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Figure 5.10: Comparison between the contributions for n~ production of the mechanisms 
of Fig. 15.11 The meaning of the lines is the same of Fig. 15.91 The re-scattering term now 
comprises also the Weinberg- Tomozawa mechanism. The data points are from Ref. 
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Figure 5.11: The same of Fig. 15.91 but for pp — > pnir + . The data points are from Ref. (2 



5.5 Importance of the different orbital contributions 
in the three-body final states 

For all the production reactions, the Ss and Sp (NN) 7r final states were found to 
be the dominant ones (solid line in Fig. I5.1H| Fig. 15.141 and Fig. 15.151 and dotted line in 
Fig. 15.141 and Fig. 15.151 respectively). 

In particular, for ir° production, the cross section with only Ss final states (solid line 
in Fig. I5.1HJ) is enough to describe the data in a wide energy range. For charged-pion 
production reactions, the cross section results from an interplay between the Ss final 
states (which dominate for 7r~ production) and the Sp final states (which dominate for 
7r + production) (solid and dotted line in Fig. 15.141 and Fig. I5.15|) . 
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Figure 5.12: Effect on the cross section of not considering the A as an explicit degree of 
freedom. The solid line is the full calculation of Fig. l5.9l Fig. r5.11l considering the mechanisms 
of direct-production, re-scattering, A-excitation and Z-diagrams (Imp. + Rescatt Aff . + A + 
Z). The dotted line corresponds the calculation with only the direct-production, re-scattering 
terms and Z diagrams (Imp. + Rescatt. + Z). 
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pp -> ppji 




280 285 290 295 300 305 310 315 320 325 



Figure 5.13: Relative importance of the (NN) n final states in pp — > ppn°. The cross sections 
including only the Ss, Ps and Pp (NN) n final states are ergs (solid line), a Ps (dashed-dotted 
line) and a Pp (dotted line), respectively. The data points are from Ref. |23j. 
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Figure 5.14: Relative importance of the (NN) tx final states in pn — > ppix~ . The cross sections 
including only the Ss, Sp ,Ps and Pp (NN)ir final states are as s (solid line), as p (dotted 
line), <jp s (dashed-dotted line) and ap p (dashed line), respectively. The data points are from 

Ref. jig. 



119 



120 



Charged and neutral pion production reactions 



pp -> prui + 
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Figure 5.15: Relative importance of the (NN) n final states in pp — > pnii + . The meaning of 
the lines is the same of Fig. 15.141 The data points are from Ref. [28J. 
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pp -> ppt 



pn -> pp7i 
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Figure 5.16: Ratio between the cross section with different energy prescriptions and the 
TOPT cross section. The dotted, dashed-dotted and dashed line corresponds on-shell 
approximation, fixed-kinematics and static approximation, respectively. The solid line is 
the calculation within the S-matrix approach. The Ohio interaction is used for the NN 
distortions. 



5.6 Approximations for the energy of the exchanged 
pion 

In this section we start by showing that for all charged pion reaction channels the S- 
matrix approach for the effective operators, which describe the physical meson production 
mechanisms, is a very good approximation to the DWBA result obtained directly from 
TOPT. 

We compare in Fig. l5.1(j1 the result from time-ordered perturbation theory to the total 
cross sections of the several charge production reactions obtained with the on-shell, fixed- 
kinematics, static and S-matrix approximations. The validity of these approximations has 
been investigated previously (HH] but only for the re-scattering operator and 7r° production 
(see Chapter Hj). 

For all the production reactions, the deviation of the S-matrix approach result from 
the TOPT cross section is much less than the deviations obtained with the other approxi- 
mations. The maximum deviation, on the other hand, occurs for the static approximation 
(as checked in Chapter 01 at variance, below threshold, as expected, the static approxima- 
tion performs better, giving a basis for the traditional static interactions at low energies.) 



We compare in Fig. 15.171 Fig. 15.181 and Fig. 15.191 the effects of the on-shell, fixed- 
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Charged and neutral pion production reactions 



kinematics and static approximations on the total cross section to the result within the 
S-matrix approach. As expected, these results are in agreement with those shown in 

Fig. EH 



pp -> ppn (Energy prescriptions) 
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Figure 5.17: Ration between the cross section forpp — > ppn with different energy prescriptions 
for the production operator and the S-matrix (coinciding with the "exact" TOPT result, as 
shown in Fig. I5.16| cross section. The dotted, dashed-dotted and dashed line corresponds 
on-shell approximation, fixed-kinematics and static approximation, respectively. 
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pn -> ppjr (Energy prescriptions) 
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Figure 5.18: The same of Fig. 15.171 but for pn 



ppir~ 



The deviation from the cross section calculated within the S-matrix approach is max- 
imum for the static approximation, which considers the energy of the exchanged pion to 
be zero (dashed line). For energies close to threshold, the static approximation overesti- 
mates the cross section by a factor larger than 2. The on-shell approximation ("E — E' n 
approximation, which considers the energy of the exchanged pion to be the difference 
between the nucleon on-shell energies before and after pion emission) deviates also from 
the reference result, but it may be off 20% for all the charge channels, even for energies 
close to threshold. 

For the 7r° and n~ production cases (Fig. 15.171 and Fig. I5.18| respectively), about the 
same happens for the fixed kinematics approximation (energy of the exchanged pion set 
to its threshold value, m 7r /2) as seen in the dashed-dotted lines. As for n + production, 
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pp -> prm + (Energy prescriptions) 
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Figure 5.19: The same of Fig. 15.171 but for pp — > pnn + . 



represented in Fig. 15.191 the cross section within the fixed-kinematics approximation is 
underestimated by a factor of 1.5 — 2.5 in the region near threshold, and, for higher 
energies, by a factor of ~ 5. 

Our results for 7r + production are consistent with discrepancy with the experimental 
data reported in Ref. |SE], where the calculation also used the 7rN amplitude, but 
the fixed kinematics approximation was assumed 3 , leading to a calculated cross section 
which underestimated the data by a factor of 5 or by a factor of 2, depending on the 
A contribution being included or not includedjOEl- We verified these reduction factors, 
when the same channels were introduced in our calculation. We furthermore concluded 



3 Also, only two NN channels were considered, namely the 3 Pi — ^ 3 Si and 3 Pn - So- The Argonne 
V18 and the Rcid93 potentials were employed for the NN interaction. 
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that the deviations observed for the ir + production reaction are originated by the p-waves 
contributions to the ttN amplitude, as illustrated in Fig. 15.201 

In other words, the crucial role of the A in ir + production enhances the importance of 
the 7rN p-wave states resulting in a slower convergence and in a higher sensitivity to the 
fixed kinematics approximation. 
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Charged and neutral pion production reactions 



pp -> prur + (J=0) 
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Figure 5.20: Dependence of the effect of the fixed kinematics approximation for the cross 
section for pp — > pnir + on the orbital it (NN) angular momentum. The solid line is the full 
calculation. The dashed line is the cross section considering only s-wave contributions 
to the ttN amplitude. The dotted line with +'s further assumes the fixed kinematics 
approximation. The Ohio interaction is used for the NN distortions. 
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5.7 Conclusions 

This Chapter investigated charged and neutral pion production reactions within the 
S-matrix approach. The main conclusions are: 

1) the S-matrix approach for the description of the pion production operators repro- 
duces numerically the DWBA (TOPT) result 4 . We note that the S-matrix approach 
simplifies tremendously the numerical effort demanded in TOPT by the presence of log- 
arithmic singularities in the pion propagator for ISI, 

2) the contributions from the impulse, re-scattering, A-isobar and Z-diagrams, within 
the S-matrix approach successfully describe the cross sections not only for neutral but 
also for charged pion production, near threshold, 

3) the newly developed realistic Ohio NN interaction tuned well above pion produc- 
tion threshold was found to improve the description of the experimental data. For 7r + 
production, its effect is crucial to reproduce the data in the near-threshold region, 

4) the reason for previously reported (BE] failures to describe the pp — > pnir + reaction 
with the %PT 7rN amplitude (giving strengths for the cross section from 1/5 to 1/2 of 
the experimental one) stems from the approximation taken then for the energy of the 
exchanged pion, 

5) convergence of the partial wave series for neutral pion production, and for energies 
up to 30MeV above threshold, does not require channels with total angular momentum 
higher than J = 0. However, for the charged 7r~ and tt + production cases and at the same 
energies, channels with J = 1 and J = 2 at least are needed, due to the importance of p 
waves in the irN and NN interactions, and to the admixtures of both T = and T = 1 
NN isospin states. 

We add that a more detailed description of the experimental data at higher energies 
may further require a complete coupled-channel NA calculation. The S-matrix approach, 
shown here to be successful near threshold, should naturally be tested also by the more 

4 As this thesis was being concluded, an investigation (arXiv:nucl-th/0511054) identified within xPT 
the irreducible contributions of loop diagrams to be included consistently with the distorted NN wave 
functions. That reference provides the general recipe of calculating the energy dependence of the vertices 
by using its corresponding on-shell value, in agreement with the optimal prescription that we found 
numerically in this work. 
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sensitive polarisation observables and, above all, in high precision calculations, as for 
instance in charge symmetry breaking studies. 
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Chapter 6 
Conclusions 



Understanding pion production, in all its possible charge channels, in nucleon-nucleon 
collisions near threshold has been a challenge over the last decades. With the experimental 
data from the new cooler rings, the field entered in a new domain of precision. These 
unprecedented high quality data on unpolarised, as well as on polarised observables, 
triggered a plethora of theoretical investigations. 

The physics of meson production in nucleon-nucleon collisions is very rich. As the 
first strong inelasticity for the NN system, the phenomenology of the NN — > NNw 
reaction is closely linked to that of elastic NN scattering. The fact that low- and medium- 
energy strong interactions are controlled by chiral symmetry led the hope that chiral 
effective theories could be used to analyse these processes, and achieve a fundamental 
understanding of the production processes. This provides the opportunity to improve 
the phenomenological approaches via matching the chiral expansion and inversely, to 
constrain the chiral contact terms via resonance saturation. This early excitement was 
quickly abated by the realisation that proper evaluation within x?T involves surmounting 
several severe difficulties, which are caused by the high-momentum transfer nature of this 
process. 

The cross section for the reaction pp — > ppn near threshold was especially elucidative, 
since the main pion exchange contribution is there ruled out by isospin conservation. It 
clearly showed that the transition amplitude results from a delicate interference between 
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the single-nucleon term, the pion re-scattering term, and various additional and individ- 
ually important contributions from shorter range mechanisms ( "Z-diagrams" from heavy 
meson exchanges). To establish the importance of the latter in an overall description of 
the different charge channels of pion production, was the main motivation for this work. 

DWBA approaches apply a three-dimensional quantum-mechanical formulation for 
the NN distortion, which is not obtained from the Feynman diagrams describing the 
production mechanisms. As a consequence, the energy dependence of the pion production 
operator has been treated approximately and under different prescriptions in calculations 
performed so far. A clarification of these formal issues was thus needed to enable sound 
conclusions about the physics of the pion production processes, and it was the aim of 
Chapters El and EJ 

Moreover, advance in understanding pion production in NN collisions must follow 
not from the exclusive study of pp — > ppn , but also from the global understanding of the 
effects that afflict all channels. Also, to pin down the different production mechanisms, 
one needs to go beyond the lowest partial waves. These were the main motivations for 
the work of Chapters |U and 

In short, this thesis aimed to address the problem of charged and neutral pion pro- 
duction in nucleon-nucleon collisions. The main conclusions are: 

1. From field theory to DWBA 

The effective pion re-scattering operator was obtained starting from the corre- 
sponding four- dimensional Feynman diagram. By integrating over the energy of the 
exchanged pion, the amplitudes were transformed into those following from time- 
ordered perturbation theory and could be identified with six time-ordered diagrams. 

From these diagrams, the "stretched box diagrams" (i.e., those with more than 
one meson in flight in the intermediate states), were seen quantitatively to give 
a very small contribution, and therefore were neglected. The remaining diagrams 
were identified with the distorted-wave Born approximation (DWBA) amplitude, 
allowing to extract an effective production operator. The calculations considered a 
physical realistic model for pions and nucleons and their interaction vertices. 
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2. Energy prescriptions for the production operator 

The DWBA amplitude originated by the decomposition of the Feynman diagram 
into the corresponding time-ordered perturbation theory (TOPT) terms, was used 
as a reference result to study the commonly used approximations (on-shell, fixed- 
kinematics and static) for the energy of the exchanged pion. This energy, not fixed 
in the quantum-mechanical three-dimensional, non-relativistic formalism, was fixed 
as ou n , the on-shell energy of the exchanged pion. 

Near threshold, the retardation effects in the pion re-scattering mechanism were 
seen to be not very decisive. On the other hand, the effect on the re-scattering 
vertex alone was found to be significant, largely overestimating the cross section, 
as unnaturally obtained in previous calculations with ad-hoc choices for the pion 
energy. 

3. The S-matrix approach 

When applying TOPT, each of the re-scattering diagrams for the initial- and final- 
state distortion defines a different off-energy shell extension of the pion re-scattering 
amplitude, since energy is not conserved at individual vertices. This imposes the 
evaluation of two different matrix elements between the quantum-mechanical NN 
wave functions, for such operators. 

Within the S-matrix technique, however, after the on-shell approximation is made 
consistently, the pion re-scattering parts of FSI and ISI diagrams coincide and one 
can identify them with a single effective re-scattering operator. 

The amplitudes and cross sections obtained with the S-matrix effective operator 
were verified to be very close to those calculated with the "exact" TOPT result, up 
to 150MeV above threshold. 

4. Charged and neutral pion production 

If one considers the contributions from the single-nucleon emission term, re- 
scattering, A-isobar and Z-diagrams, the S-matrix approach successfully describes 
the cross sections for both neutral and charged pion production, near threshold. 
The S-matrix approach reproduces numerically the results from TOPT. 
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For neutral (charged) pion production, the summation of all channels up to total 
angular momentum J = ( J = 2) is enough up to 30MeV. 

The NN interaction was found to be decisive to reproduce the data in the near- 
threshold region (in particular for ir + production). In this channel we found also 
that the correct treatment of the exchanged pion energy in the re-scattering term 
beyond the fixed kinematics approximation, is essential to remove the previously 
reported too low cross section strength. 

In the overall, even for energies close to threshold, the cross sections with the ap- 
proximations for the energy of the exchanged pion other than the S-matrix, deviate 
at least 20% from the full calculation, being larger for the 7r + channel. 

Our findings for the cross section are in agreement with those of the Julich model. 
However, contrarily to this approach, for which all the short range mechanisms are in- 
cluded only through ^-exchange and adjusted to reproduce the total cross section of the 
reaction pp — > ppn close to threshold, in our calculations no adjustment is made, as the 
parameters for the Z-diagrams are taken consistently from the NN interaction employed, 
providing a good description of all charge channels. 

There are several issues which would undoubtedly be a step further in the investiga- 
tions of this thesis. Since the importance of the A contribution increases with energy, 
a detailed description of the experimental data in regions of higher energy may further 
require a complete coupled-channel NA calculation. Also, although the effect of other 
higher mass nucleon resonances was found to be small, these mechanisms should also be 
considered in future calculations. 

Moreover, the S-matrix approach, which successfully describes the near-threshold cross 
sections for neutral and charged pion production, should be submitted to the most strin- 
gent test of describing pion (and other mesons) polarisation observables. Finally, since 
the S-matrix approach passed so well the test of pion production, it is natural to extend 
the method to the rj meson sector. 

In recent years, the improvement of experimental facilities, detectors and accelerators 
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permitted the acquisition of a great deal of high-precision data, leading to the improve- 
ment of its theoretical understanding (both in x?T and in other phenomenological ap- 
proaches), and are a promise of very exciting times to come in the near future of meson 
production physics, where for instance, rj and rj production still pose interesting problems. 
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Appendix A 



General remarks on kinematics 



Fig. rO illustrat 




Figure A. 1: Illustration of the choice of variables for NN —> NNn. The solid lines and dashed 
lines represent the nucleon and pion, respectively. The initial(final) nucleon four-momenta are 
Pi and P2{U\ and [7 2 ). The four-momentum of the produced pion is Q. 



Final and initial momenta 

The initial nucleon four-momenta in the overall centre of mass system are given by 

Pi = (E u p) (A.l) 
P 2 = (E 2 ,-p), (A.2) 

with E lj2 = \/M 2 + p 2 = E. For energies close to the pion production threshold, the 
non-relativistic expressions for the energy of the final nucleons can be used. Thus, in the 
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General remarks on kinematics 



centre-of-mass of the final three-body system, 



U 2 



J? a W 

fi, Qu — y 



F 2 , ~q u ~ 



(A.3) 
(A.4) 



where q u is the relative momentum of the final nucleons and q n is the centre-of-mass 
momentum of the produced pion. From energy conservation it follows that 



2E = F 1 + F 2 + E T 

2 2 

~ 9M + — + -l 
2/i AM 



(A.5) 
(A.6) 



Laboratory energy and maximum pion momentum 

For the kinematics of the pion production reactions investigated in this work, the nu- 
cleons were considered as having equal masses (M). The laboratory energy can then be 
determined from 



Tr 



lab 



s-AM 2 
2M ' 



with s = (Pi + P2) 2 ■ The excess energy Q is given by 

Q = - (2M + m^) . 



(A.7) 



(A.8) 



The maximum pion momentum, which occurs when q u = 0, can be calculated (in units 
of m n ) using 

1 r 1 1 : 1 

(A.9) 



V 



s -ml- AM 2 ) - 16M 2 ml 



Eqs. (lA~7j) . fOjl and (lA~9j) can be easily generalised for the case when the nucleons have 
different masses 



T, 



lab 



S-M&-M&- 2M a M l2 
2M a 



(A.10) 



Q = y/s - (Mfi + M }2 + m n ) 



(A.11) 



and 



1 

77 = 



(A.12) 



J^^J (s - mi - M 2 ) 2 - AMjml 

where the masses of the initial (final) nucleons are Mn and M i2 {Mfi and Mf 2 ), and 
M f = M/i + M /2 . 
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It is important to remark that the relativistic and non relativistic calculations of the 
threshold laboratory energy T t hr. differ significantly, 

p „, 4M 2 + 4Mm n + ml- AM 2 
Tgf; = (A. 13) 

= T ™'+(^), (A.15) 

since (j^j ~ lOMeV. This difference can have decisive effects on the results, especially 
at threshold. 
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Appendix B 

Kinematic definitions for the 
diagrams 



Re-scattering diagram 



The kinematic conventions adopted for the re-scattering diagram are represented on 
Fig. IB. II The relations resulting from four-momentum conservation are listed on Ta- 
ble EU 
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* 1- 
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a) 



b) 



Figure B.l: Illustration of the choice of variables for the re-scattering diagram for a) FSI and 
b) ISI. The solid lines and dashed lines represent the nucleon and pion, respectively. The 
initial(final) nucleon four-momenta are Pi and P2(U\ and {7 2 )- The four-momentum of the 
produced pion is Q w . The four-momentum of the exchanged pion(a meson) is Q'(K). The 
four-momenta of the intermediate nucleons for the FSI(ISI) case are Qi and (^(^i and S 2 ). 
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Kinematic definitions for the diagrams 



FSI 

Pi + Q' - Qn - Qi = 
P2 - Q' - Q2 = 
Qt -K-U 1 = 
Q 2 + K - U 2 = 



ISI 

Pi - K - Sx = 
P 2 + K-S 2 = 
S 1 + Q'-Q„-U 1 = 
S 2 -Q'-U 2 = 



Table B.l: Four- momentum conservation for the re-scattering diagram. 



For the final state NN interaction (diagram a) of Fig. IB.lj) . the expressions for the 
on-shell energies of the particles are: 



OJ-n 



ml + 



-P + y + 9fc 



^2,1 



M 2 




ml + \Qk ~ Qu I 



F 2 ,i = M + 



2M 



<?7T 



exchanged 7r 

intermediate nucleons 

a exchanged 
final nucleons 



and, for the initial NN interaction, the corresponding expressions are 

exchanged 7r 



mi 



Ul,2 



M 2 + i^r 

m l + \~P - Qk\ 



F 2 i = M + 



2M 



intermediate nucleons 
a exchanged 

final nucleons 



(B.l) 

(B.2) 

(B.3) 
(B.4) 



(B.5) 

(B.6) 
(B.7) 

(B.8) 



where p and are the initial nucleon three-momentum and the emitted pion three- 
momentum, respectively. The relative three- momentum of the two intermediate (final) 



nucleons is qk(q u ), an d the energy of the emitted pion is E n = ym^ + \ q n \ . The 
nucleon, the pion and the a mass are M, m, and m a , respectively. All quantities are 
referred to the three-body centre-of-mass frame of the irNN final state. 
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FSI 




ISI 




Pi 


- Qn ~ Ql = 


Pi 


- K - Si = 





Qi 


-K-U 1 = 


Si 


-Q«-U! = 
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P2 


+ K -U 2 = 


P2 


+ K-U 2 = 






Table B.2: Four-momentum conservation for the direct-production diagram. 

Direct-production diagram 

The kinematic conventions adopted for the direct-production diagram are represented 
on Fig. IB. 21 The relations resulting from four-momentum conservation are listed on 
Table El 
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Figure B.2: Illustration of the choice of variables for the direct-production diagram for a) FSI 
and b) ISI. The meaning of the variables is the same of Fig. IB. II 



The expressions for the on-shell energies can be easily obtained from Eq. (JBJJ)- 
Eq. (JR4J) and Eq. (|B~5l)-Ea. (IB~8l) . 
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Appendix C 

Partial fraction decomposition and 
TOPT 



For the loop diagrams considered in this work, we have in general to evaluate integrals 
like 

M = J ^-f (Q' Q ) x U X L X x U 2 L 2 x U 3 L 3 (C.l) 

where C/j, Lj are terms with the poles in the upper-half plan and in the lower half-plan, 
respectively, which can be genericaily written as: 

Qo ^ Qo ~\~ i& 

and / (Qq) is a generic vertex function which depends on the momentum of one of the 
exchanged particles. 

Direct integration of Eq. (|C.1J) by closing the contour in one of the half-plans and 
picking up the poles enclosed, give as a result an expression where one has to evaluate 
the vertex / for three different energies (the pole contributions) of the exchanged particle 
Q' Q . For instance, if we close the contour in the upper half-plan, / must be calculated for 
Q' = «i, Q'q = u 2 and Q' = n 3 . We thus gain little insight on the adequate choice to 
be done for Q' in the traditional three-dimensional formalisms, where this energy is not 
fixed. 

However, by performing a partial fraction decomposition of the integrand before in- 
tegrating in Q' , provides us a visible connection with time-ordered perturbation theory. 
The way the partial fraction decomposition is done ensures that we have only contribu- 

143 



144 



Partial fraction decomposition and TOPT 



tions from terms with one pole on the on-shell energy of the particle interacting through 
/ (Qo) and the remaining one or two poles in the other half-plan. Closing the contour in 
the half-plan where there is only one pole, we can relate directly Q' at the vertex with 
the on-shell energy of the considered particle. Of course, the final result must not depend 
on the method adopted for integration. 

Let u 3 and l 3 be the poles we want to isolate. We start by writing the product U 3 L 3 
difference of two terms 

U 3 L 3 = a 33 (U 3 - L 3 ) with a 33 = — — - , (C.3) 

u 3 - h 

and doing the same for the terms concerning the other exchanged particle, U 2 L 2 : 

U 2 L 2 = a 2 2 (U 2 - L 2 ) with a 2 2 = — — r (C4) 

u 2 - h 

Then Eq. flC.l|) reads 

M = J ^rf(Q' )a 22 a 33 U 1 L 1 U 2 U 3 - J (Q' ) a 2 2a 33 U 1 L 1 U 2 L 3 (C.5) 

-J j^f(Q , )<X22<* 3 aU 1 L 1 L 2 U 3 + J j^fiQ^^assU^LzL, (C.6) 

Now, in ()C.5|) . we apart L1U2 and in (IC.6I1 . we apart UiL 2 , 



M = J ^rf(Q l Q )a 2 2a 33 a2iU 1 U 2 U 3 (C.7) 
' dQ' 



(27T) 

dQ'o 



f (Qo) a 2 2a 33 a2iUiLiU 3 (C.8) 
f(Q'o)a 22 a 33 a 21 U 1 U 2 L 3 (C.9) 



(27T) 

+ J ^rf(Q' Q )a 22 a 33 a 21 U 1 L 1 L 3 (C.IO) 

- J ^f(Q' G )^ 22 a 33 a l2 U 1 L l U 3 (C.ll) 

+ J j^f '(Q^^a^a^L.U, (C.12) 

-f (Q' )a22<y 33 a 12 U 1 L 1 L 3 (C.13) 



dQ'p 

(27T) 

dQ'p 

(27T) 
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/ (Qo) ol 22 ol 33 ol\ 2 L 2 L\L 3 (C.14) 
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and finally, in (jTTSf . (^TToj) . (f(TTTj) and (jDHj) , we apart ET^: 

^ = / ^/Wo)«22«33«2lt/lf/ 2 t/ 3 (C15) 

4-4/ (Qd) Q;2 2 q; 3 3Q;2iQ:ii?7iC/ 3 + / tt^t/ (<5d) a22«33a2iaii-^iE/"3(C16) 

(27T) J (27T) 

777% (Q'o) ^a^U^L, (C.17) 

TTj^j/ (Q'o) «22«33tt2lttllE r iL3 - ^ TTj^v/ (Qd) a 22 a 33 a 21 a u LxL 3 (C. 18) 
/(<5d)«22a33«12ail^lf/ 3 + / TTT^f (Q'o) «22«33«12ail^l^ 3 (C.19) 

J ( 27r ) 



(27T) 

f(Q' )a 22 a 33 a 12 L 1 L 2 U 3 (C.20) 



(2tt) 

+ y (Qo)oi 22 a 33 a 12 a n UiL 3 - J j^pjf (Q'o) a 2 2a 33 ai 2 a 11 LiL 3 (C.21) 

- / ^/(Qd)«22«33«12^2^1^3 (C.22) 

where, accordingly to our conventions, an = — 



A4 is written as a sum of eight types of terms: terms with three-poles in the upper- 
half plan (JC.15j) or in the lower half-plan (jC.22jl : with two poles in the upper half-plan 
(first term of ()C.16Jl and ()C19Jl ) or in the lower half-plan (second term of (jC.18j) and 
(IC.21|) ): with one pole in each half-plan (second term of (|C16|) and (jC.19|) and first 
term of (jC.18|) and ()C.21)l ): and finally, terms with two poles in the upper half-plan and 
one in the lower half-plan (jC.17j) or vice- versa ()C.20j) . The corresponding schematic is 
represented on Fig. 13.21 

Assuming that the integrals with all the poles in the same half-plan vanish (which 
happens if / (Q' ) is linear in Q' ), M. is reduced to only six terms: 

M = J J^rf (Q' )oi 22 a 33 a 21 a n L 1 U 3 + j J^rf (Qo)oi220i330i 12 auLxU 3 (C.23) 

+ / (Q'o) «22a33«21«llE r i^3 + / (Qd) a 22«33ai2«ll^l-^3 (C.24) 

J (27TJ J (27T) 

+ J j^f (Q'o) a^auL^U, - J ^rf(Q' Q )a 22 a 33 a 2l U l U 2 L 3 (C.25) 

These terms can be directly identified with the ones resulting from TOPT. More precisely, 
(jn2HJ)- (f(T24|) can be identified with the DWBA terms and (f(T25jl with the stretched 
boxes, as is illustrated on Fig. 13 .41 for the FSI case. 

The next section exemplifies this method for the FSI case. 
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Partial fraction decomposition and TOPT 



Re-scattering diagram with final state interaction 



From Eq. ()3.5|) we have 



Ui = E 2 — uj 2 + is 

u 2 = E 2 - F 2 - u) a + ie 

u 3 = —to n + ie 



and 



l\ = E % + ui — Ei — ie 
l 2 = E 2 - F 2 + uj a - ie 
h = - is 



(C.26) 



and consequently, 



1 

a u 
1 

1 

"33 
1 

1 

«21 



ui - k = E tot - E w - u) X - lo 2 + 2ie 



u 2 -l 2 = —2u a + 2ie 



= u 3 - l 3 = —2u n + 2ie 



u 2 -h = E tot - F 2 - E w - Ui - u a + 2ie 



(C.27) 
(C.28) 
(C.29) 



ui - l 2 = F 2 - lj 2 - u a = E tot - F 1 - E w - lu 2 - Ld a + 2ie (C.30) 



(C.31) 



with E to t — E\ + E 2 = F 1 + F 2 + E n . The residues of Ui, U 2 , L x and L 2 are 



Res U\\ QI . 



1 



1 



1 



r ° h h - ui E 2 -u 2 - u n E tot ~ E 1 - u 2 - co n 
Res f/ 2 | Q , = , 3 = 



li 2 -^2 — F 2 — LOo- — U n E to t — Ei — F 2 — LU a — u n 



(C.32) 
(C.33) 



and 



Res LiL,_„ = = = (C.34) 

w o-«3 U3 _\ x Ex-E^-uji- lo^ E tot - E 2 - E n - ui - u n 

11 1 

Res L 2 \ Q , o=U3 ~ U3 _^~ F ^_ E ^_ u ^_ u ^- Em F\ E n E 2 uj a LO n 
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Therefore, the amplitude Eq. (|3.5J) reads: 

M™? = -gl [ x (C.3f 

J (27r) <kOaUn 

/K) 



. (E to t - 


- E^ - 


-LOx - 


- U} 2 ) 


(Etot - 


-Ex- 


- L0 2 - L0 n ) (E tot - 

fM 


-Fx- 


- E w - 


-LU 2 - 


" Wo) 




(E to t 


— E n 


- LOx 


- L0 2 


) (Etot 


-Ex 


- u 2 - Uk) (Etot 

f (-^7r) 


-F 2 


— E w 


- LOx 


~ Wo) 




(Etot 


— E n 


- LOx 


- L0 2 


) (Etot 


— E 2 


— E n — LOx — UJ ni 
/(-^tt) 


(Etot 


-Fx 


— E n 


— L0 2 


- u a ) 


(Etot 


— E n 


- LOx 


- L0 2 


) (Etot 


— E 2 


— E n — LOx — UJn, 

/K) 


(Etot 


-F a 


— E w 


- LOx 


- u a ) 


(Etot 


-E x 


- L0 2 




) (Etot 


-F 2 


-E^ -L0 X - Uo) 
f (-^tt) 


(Etot 


-Ex 


-F 2 


~ W n 


-Uo) 


(Etot 


-E 2 


— E n 


- LOx 




(Etot 


— Fx — E n — lo 2 


- u a ) 


(Etot 


-F 2 


-Fx 


— E n 
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Appendix D 

T-matrix equations and phase shifts 



NN propagator 

The NN propagator, 



GnN ~ 2 (Eo -E) + is (D ' i: 



is given for non-relativistic nucleons by 



Z I 2M 2M I + 16 



where E = The nucleoli- nucleon relative momentum (on-shell momentum) is k (ho). 
Eq. (|D.2|) formally coincides with the Blankenbecler-Sugar propagator |1UU| HUlj consid- 
ered with minimal relativistic definitions of the amplitudes. The relativistic nucleons can 
be described by the Thompson propagator fl()2j . 

r R - 1 E ° + — - 1 E ° + — rn %\ 

^ NN ~ 2 {El - EP) + is ~ 2 (k 2 - k 2 ) + ie 1 ' 

since E = V 1 M 2 + k 2 . In a short-hand notation, 

c Hk) mi) 

Gnn - m - ^) + is { ] 

with j3(k) — M for the propagator of Eq. (jD.2|) and j3 (k) = E+E ° for the propagator of 
Eq. 
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T-matrix equations and phase shifts 



Uncoupled channels 

Equations for T 

For uncoupled channels, the equation for T is calculated from the potential V using 

T (k, k ) = V (k, k ) + J V (k, k') G NN (k') T (k', k ) k' 2 dk' (D.5) 

where ko is the on-shell momentum. The ie term of the denominator of Gnn in Eq. (|D.5j) 
may be evaluated by applying the Cauchy principal value theorem |lU3j . Numerically 
we used also the subtraction method with regularisation (see Appendix |FJ) to calculate 
T(k,k ): 

'g(k,k') 



T(k,k ) = V (k, ko) + V [ 9 (fc ' k 'l 9 - i fc ' ko) k' 2 dk' - tirSRes 

J K~ k! 1 



h 2 — h' 2 

A^Q rv 



(D.6) 

k'=k Q 



Here, g (k, k') = V (k, k') (3 {k') T (k', k ) k' 2 is the regular part of the integrand. Using 
Eq. (|D.4j) . Eq. (|D.6|) transforms into: 

,/2 



T (k, k Q ) = V (k, k ) + J V (k, k') (3 (k') T (k', k ) ^J^dk' (D.7) 

- J V (k, k ) (3 (k ) T (k , k ) ^Aj^dk' (D.8) 

- iir^-V(k,ko)f3(k )T(k ,ko) (D.9) 
The corresponding equation in the discrete form is 

T % o = Vi + J2 [VyTjokfijWj] - V l0 k 2 (3 (k ) ( V t^-^w, ) T 00 (D.10) 

+ i<* Woo (D.ll) 

with a = — z^pa for the non-relativistic case (Eq. ()D.2Jl ). and a = -ir^- for the rela- 
tivistic one (Eq. (|D.3ft ). Defining My = 8y — My, with 

VijGjtfwj j ^ 

Ma = I ( \ , (D.12) 

V i0 k 2 P (k ) E ra^' I + i^Vio j = 



the T-matrix can be calculated from: 



MijTjo = V i0 => T j0 = M~ V fe0 . (D.13) 



3 
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Phase shifts 



-e sm6=>tg(26) = z 



(D.14) 



Coupled channels 

Equations for T 

For each value of J ^ there are three possible values for the orbital angular momentum, 
L = J and the coupled L = J ± 1. The equations for the coupled case are 



T = V + V G T~ 
T+- = V + ~ + V + ~G Q T- 
T-+ = V~ + + V~ + G T- 
T++ = V++ + 1/ ++ G T- 



V ++ G T+- 
V"G T~ + 
V + -G T-+ 



(D.15) 



which reads in the matrix form 



-M] 



M + ~ 



T 



jo 



T 



jo 



7 1 — ! 

J j0 









V + 



(D.16) 



Phase-shifts 

Stapp, Ypsilantis and Metropolis parametrisation 



In the Stapp, Ypsilantis and Metropolis |1U4| parametrisation, the S'-matrix is written as 
a product of the following three matrices 



S 



From 



and 



e ld - \ / cos2e z sin 2s \ / e ld ~ 
e fs+ J \ isinZe cos2e I I e f&+ 

cos(2e)e^- zsin (2e) e^ 4 "*") 
is\n(2e)e i ^ ++ ~ 5 -) cos (2e) e i2 ^+ 

S~+ + S + ~ = 2z sin (2e) e*^^-) 



S++S" =cos 2 (25)e i2 P ++;f - 
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cos (2e) e 



(D.17) 
(D.18) 

(D.19) 
(D.20) 
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T-matrix equations and phase shifts 



it follows that 



and thus 



tg (2e) 



1 s- + + s + - 

2 y/S++S— 



1 _ » s-++s^ 
-lln / 2 ^s+^s 1 



5 + and <5_ are now easily calculated from 



5. 
7h 



1 

2 aig 



2 ar § 



1 g~ ++5-> 

2 VS++S- 



s— 

cos (2e) 
cos (2e) 



Blatt and Biedenharn parametrisation 



In the also used Blatt and Biedenharn parametrisation , 5 is written as 



-i 



S 



cos e sin £ 
— sin e cos £ 



e i2S + 



COS £ sin £ 
— sin £ cos £ 



cos 2 (e) e i2<5 - + sin 2 (e) e i2S + sin e cos £ (e i2<5 - - e i25+ ) 
sin £ cos £ (e m ~ - e i25+ ) sin 2 (e) e i2S ~ + cos 2 (e) e i2<5 + 



Defining 



and 



+ 5 ++ e l2S - + e l26 + 



^ppmm 



Sj \-S—— ~\~ Sj S- 



3 j2(<5++<5_) 



(D.21) 



(D.22) 



(D.23) 
(D.24) 



(D.25) 
(D.26) 

(D.27) 



(D.28) 

(D.29) 
(D.30) 



it is straightforward to obtain 



Crf 



i28- 



m ^ppmm ^mx 



(D.31) 



Cdp ^ppmm '-mi ^ j (D.32) 
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and thus giving for the phase-shifts 5- and 5+ the following expressions: 



2<J + ) = i [ f (D.34) 

2 - Im [ C dp] 



For both the uncoupled and coupled cases, S = 1 + iaT. 
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Appendix E 

Partial wave decomposition of the 
amplitudes 



Re-scattering diagram 

The re-scattering diagram amplitude involves terms like (a^' ■ q) . Since q is the three- 
momentum of the exchanged particle, it can be expressed as 

q = Aq af - Bq bf + Cq ej (E. 1 ) 

where q b f (q a f) is the relative momentum of the two nucleons before (after) pion emission 
and q e j is the momentum of the produced pion. A, B and C are numerical coefficients 
which depend on the reaction considered 1 . 

With xi = cos Z (q e j, q a f — q b f) and x 2 = cos Z (q a f, the amplitudes can be written 

as 

Mr = F(x 1 ,x 2 ) (a {i) -q) (E.2) 

where T is a function depending on the kinematic variables of the diagram considered. 
Using Eq. (lFL43|) and Eq. (ITT591) . Eq. (jQ reads 

Mr = -V4^[a^®Y 1 (q)]° Q \q\^(x 1 ,x 2 ) (E.3) 
= -(47r)^(xi,x 2 ) p^^^/^ft/^)]] (E-4) 

h+U=l 

F(l;l 3 k)\Aq af - Bq bf \ h (Cq ej ) h (E.5) 

1 Amplitudes should be defined such as A = 1. For the cases of A ^ 1. then x\ — 
cos Z(q e j,Aq a f ~ Bq hf ). 
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156 Partial wave decomposition of the amplitudes 



where F (I; lxl 2 ) = y ( 2 ^i+i)K^2+i)! ' ^ e f unc ^ on -F ( x u x "i) 1S now decomposed using the 
Legendre polynomials 2 



T x 2 ) = f h M P h M | with h EE y/2h + 1, (E.6) 



leading to 



Z 3 +U=l h 



F(l;y 4 ) |a« ® [^ 3/4 (Aq^Bq bf ,q^j ®y*> h (Aq af - Bfaq^ 



o 



where Eq. (IH.17J) was used. Coupling the spherical harmonics with the same arguments 
Eq. (IIL66|) and then using Eq. (jlTBIj) . Eq. (lR7l) simplifies to 



M r = -^v 7 ^ £ |Ag a/ -5g 6/ |' 3 (C ?ei ) k F(l ; y 4 )^ (E.8) 

Hi *1*2 /l+/2=*l 



|Ag o/ - Sgl,/!* 1 I / 3 ; 4 / x 



o 



The Legendre polynomials are now used to decompose f l± (x 2 ) \Aq a f — Bq b f\ 



h-ti 



G hhh (x 2 ) = f h {x 2 )\Aq af -Bq bf \ h - u (E.9) 

2 



E^ lW3tl ^ 2 (^)| (E.10) 



2 The coefficients of the expansion are / !l (#2) = / -^(^l,^) Ph ( x i) 
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and Eq. (|E.8|) reads 3 



M r = - { -^V^ Yl {Cqe j ) h F{VMk)hhY, C ^C l ^t (E.13) 



^3+^4 = 1 



tlt 2 



h t 2 1 



£ F ( tl ; hf 2 ) (Aqaf)* (-Bq bf )^ £ (-l)^ 2 

^4 ^3 ^1 I /i+/ 2 =ti ilia 



(7 



(') 



® { K/ 2 «v) ® F * 2 1 ® % 2 (&/> 9V) 



n o 



Using Eq. (lrL66|) and F* 2 



4?r 3 ; t 2 2 o (ft* qbf), one obtains 



(4tt) S 



^ (,Cq e j) k F (1; / 3 / 4 ) Z 3 Z 4 ^ ^ooo^C'ooo 1 ' 2 



(E.14) 



i 3 +i 4 =i 



tl*2 



^4 h h I /i+/ 2 =ti ilia 



9 2-^ 000 0( 



2hh2 
000 



hi h 2 h 
h fi k 



a 



(i) ® Mo (ft* 9b?) ® 3^ fa ft?)] 1 



Orbital matrix elements 



With Eq. (j!T67j) and Eq. ()H.57jl . the orbital part of A4 r reads 

X = V hhF(l;hh)(Cq ej ) h YTiCSS 

7T 

tl*2 



4 a/Itt 



'Wlt 2 

000 



(E.15) 



z 3 +z 4 =i 



tx * 2 1 f £ /i/ 2 ^(ti;/i/ 2 )(Ag a/ ) /l (-5ft/) /2 E(- 1 ) 

^4 ^3 ^1 | /i+/ 2 =*i iii 2 



il+«2 



h+t 2 +i 



h\h,2 
ti 1 t 2 
gi hi h 2 



h ft h | gi 



3 Since 



(*3+U=l) | | C'^^h+h+U even lZi 



h h 1 

h h o 

ti ta i 

Hi l/l Z 4 ^3 il 




1 



— V^(i 3 Mi* 2 ;i/i) 



(E.il) 



(E.12) 
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Partial wave decomposition of the amplitudes 



Thus, the orbital matrix elements given by 

O = (q af q ej £'M'\M°\q hf LiM) (E.16) 

= J dn qbf dn qaf dn qej yt f Z M ' (£, £) (£) M° (e.it) 

= | dn^dn^dn^D^^' (£?. 8) [(-i) M ^ (S/)] (e.is) 



are 



O = -i!^--L(-l) A/ £ ll£F(l;/ 8 /4)(C7gei) u (E.19) 

4 V47T 



-2^2 

i z 2 < 



^1 lx 1 
^4 ^3 ^1 



/l+/2=*l 

/ 2 /1 h J I £' i/ U J 

Finally, applying the Wigner-Eckart theorem |l()tij 

(q af q ej C'M'\O klX \q bf LiM) = \\O hX \\Li) (E.20) 

= f(-l)^ +M C^ A <£'||0*,a||L*>, (E.21) 

one obtains the reduced matrix elements 

{c'WMiWh) = -M^(-if +l ~ x i E WlF(i;kU)(c qej ) u (E.22) 

h+h=l 

h l x 1 



E m)^ 2 s 2 s 2 css'' ^tiav* 

E /i/ 2 F(ti;/i/ 2 )(Ag a/ ) /l (-i?g f ,/) /; 

/l+/2=tl 



^4 ^3 ^1 



r hhL Srh i 2 u \ L f Li h \ \ h 1 Z x 
O ooo °ooo S ( \ 

h h h a Lf Li 



Spin matrix elements 
From Ref. [MI, it follows that 



S f Si 1 



[S f \\*M\\ = (-i)«+«+ 8 f 2S i si^/s 1 ( Sl + 1) I 1 1 } (E.23) 

[ Sl Si s 2 
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(S f \\a^\\ Si) = 2SiS 2y /s 2 (s 2 + l) 



Sf Si 1 



S 2 S 2 Si 



(E.24) 



Direct-production diagram 

The direct-production diagram involves terms like ■ (q^ + pi + q\) S (q 2 — p 2 ), which 
have the form of 



■ {Aq ej + Bq af + Cq bf )5 



S (&/ - 9a/) - -f 



(E.25) 



where <fb/ (<7a/) is an external(loop) momentum and q e j = q n . With x = cosZ (q a f,Qbf) 
the amplitudes can be written as 



Mi = ■ (Aq ej + Bq af + Cq bf ) 5 
Using Eq. (ITT431 . Eq. (lE~26l) reads 



S (96/ - 9a/) - 



Qej 



Mi = -V4:ir\Aq ej + Bq af + Cq bf \5 



s (96/ - 9a/) - -y 



'4tt |A<f e .,- + B<f a/ + C9&f 1 S 
and thus, by Eq. (jH.59|) one has, 



S (96/ - 9aJ 



Qej 



Mi 



(4vr)5 



(96/ - Qaf) 



Qej 



J /l+/2=l 

^(i;A/2){a»®^ l/2 (9Tf,9a7)}o 



(E.26) 



(E.27) 



(E.28) 



^ (5'g a/ ) /2 {C'q bf ) h (E.29) 



Orbital matrix elements 



The orbital matrix elements are calculated from 



U = (q a fq ej C'M' \O kl x\q b fLiM) 

dn qbf dn qaf dn qe .y%f {q7f ^ €j) Y u {qTf) M . 
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(E.30) 
(E.31) 
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Partial wave decomposition of the amplitudes 



Using Eq. ()H.65j) and 



dn qej 8 



s {qbf - q a f) - y 



s=±l 



5 [\s (q bf - qaf)\ - q ~f] 

Wbf - Oaf I 2 



U reads 



(E.32) 



u = - (4tt) J dn qbf dn qaf 

91+92=1 

5 [\s (q bf ~ q af )\ ~ q ~f] 



m) 91 {sqafT 



Wbf-qafV 

£ {B'q af ) h (C'q hf ) h F{l;hf 2 ) 



F(l x , gi g 2 ) (E.33) 



\qbf-q a f\ h+f 2 =i 



Now, the part depending on a; is decomposed in Legendre polynomials 



A 



(sqbf) (sq a f) 



92 



(B'q af ) h (C'q bf ) h 5 [\s (q bf - q af ) \ - q ej ] (E.34) 

(E.35) 



1 9ft/ - Qa 

(4n)J2f hl ' h9 H-l) hl p? lh (*f,Q) 



where 4 



leading to 



U 



f hl ' fl91 = J dxAPi 



h (x) 



(4tt) 5 



/l+/2 = l 



gi+92=lx 



{Y* L f (Q) ® y;[; 2 (fa, q7f)} c M , Ytt m yf lh (§67, Q) . 



'Since Isl = 1. 



hi'figi 



(B' qi ) J2 (C'q fx ) Jl P h (x) 

QlQfx 



f 



with xo being the solution of J — 2q b fq a fXo + q\f — ^j- = 0. 
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(E.37) 



J dn qbf dn qaf F &hh) E F (!*>9i92) (e.3 



(E.36) 
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Using Eq. ()H.63j) . one obtains 

u = - { ^r f f x J dn qbf dn qaf fihF^hh) Yl 92 (e.39) 

/i+/a=l gi+92=lx 



'000 

h c 2 92 lx Lf I dld2 



d\ di 1 

h h h 



® %\ d2 yfi£* Q) Y& . 



Finally, the result of the angular integration (of the orbital part) is 5 



O = -yi/i.U E /i/ 2 F(l;/i/ 2 ) ^ 92F(l x ; 9l g 2 ) (E.43) 

/l+/2 = l 9l+92=lx 



/2 /l *1 C 2 C 1 



'2'lC2 

000 



and, therefore, the reduced matrix elements are 



O = £ TihF^hh) 92F(l x , gi g 2 ) (E.44) 

/l+/2 = l 9l+g2=lx 

^f^k^rl^ 02 91 Ie^oT^ 

h c 2 [ 92 lx Lf J d 

dlC2l \ c&» ( U d > 91 > (-1)* 



'2^1C2 

000 



5 Note that 



and also by Eq. i|H.57|) 



® yZ 2 (<s,?o}° = E ( E - 4 °) 



B = ^(^.^^^'(S.®)^ 4 ^) (E.41) 

= E (S) ^ («) (E.42) 

77117712 

E (6c) m (s) ■ 
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Partial wave decomposition of the amplitudes 



For the simplest case of a term like ■ (pi + <fi), Eq. f)E.44|) reduces to 

/ A'jr ^ 

LfLil x i 92F(l x ; gi g 2 )J2f hlx ° 91 

gi+92=lx h 



o = 

^2 x -v 



L f9 2C2 I £ ' C 2 #1 I rflhctrtLihgi J L * ^ & 

ooo S ( °ooo °ooo S 

[ 92 l x Lf J ^ c 2 £ 1 

Spin matrix elements 

The spin matrix elements are given by Eq. (lE~23l) and by Eq. (IR241 



(E.45) 



Z-diagrams 

The computation of the amplitude for the Z-diagrams involves the calculation of terms 
like ■ [a^ x <f] . Since 



a 



[a^ x g] = [<r« x c?< 2 >] ■ g 



(E.46) 



and 



•g = iv / 2v / 4^|gl{[a w ®(T {j ' ) ] 1 ®F 1 (g)}° (E.47) 
the matrix elements of the orbital part are thus i\/2 times the matrix elements of Eq. (jE.22|) . 



Spin matrix elements 

The spin matrix elements are 



SaSbSf 



(E.48) 





Sf 


Si j 


V S 1 ( S l + 1) V S 2 («2 + 1) < 


Sl 


Sl 1 > 




S2 


S2 1 



and 



SaSbSf 



{aw®*®}' 



sis 2 Si 



\j+2( Sl +S2-Si) 



4JSijsis 2 



(E.49) 



y/s 1 (si + l)y/s 2 (s 2 + 1) < 



S f Si j 

S2 S 2 1 
Sl Si 1 
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A resonance contribution 

The A contribution diagram orbital and spin matrix elements are similar to the corre- 
sponding ones for the re-scattering diagram (Eq. (jE.22|) and Eq. (|E.23[) . respectively). 



Isospin matrix elements 



The isospin matrix elements involved in the calculations are 



with t± = ^4p, and 



PP 
PP 
np 



.(i) 



-(i) 



r 



(i) 



pp 
np 
pp 



PP 
PP 
pn 



T- 



(2) 



r 



(2) 



.(2) 



PP) = 1 

pnj = 2 
pp) = 2 



pp 
np 



e i r (1) r (2) 

t abc If, 1 c 

e , r (1) r (2) 

tabc I h 1 r. 



np 

PP 



-2i 



2i 



(E.50) 
(E.51) 
(E.52) 

(E.53) 
(E.54) 



Coefficients A, B and C for the diagrams considered 
Re-scattering diagram 

The coefficients A, B and C of the partial wave decomposition of the re-scattering dia- 
gram, Eq. ([E.ip . are listed on Table |E~T1 The coefficients for the A contribution diagram 
can be easily inferred from Table IE. 11 

Direct-production diagram 

The coefficients of the partial wave decomposition of the direct-production term of Eq. (jE.25|) 
are listed in Table IE. 21 

Z-diagrams 

The coefficients of the partial wave decomposition of the Z-diagrams are in Table lE~3l (for 
a-exchange) and in Table lE~4l Table lE.61 (for u;-exchange) . 
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Partial wave decomposition of the amplitudes 





b bl ' 


Fbl ' 


Ibl ; 


Ibl ' 


q 




f+p-qk 


f - % + q u 


f + qk-q u 


q a f 


Ik 




qu 


n 

qu 


Qbf 


P 


V 


% 


qk 


A 


1 


-1 


1 


-1 


B 


1 


-1 


1 


-1 


C 


1 


l 


i 


i 


2 


2 


2 


2 



Table E.l: Coefficients of the partial wave decomposition of the re-scattering diagram. The 
b) superscript refers to the 1 <-> 2 diagram. 



qbf 

A 
B 

c 

s 



fsf 



v 



q f =p 



Ett 



2{E+M) 

Ett 



E+M 

Ett 

Ei+M 
1 



FSI b ) 



V 

q* f =p 
1 



E % 



2{F+M) 
E„ 



F+M 

E n 
E 2 +M 

-1 



isi a > 

= q u + f 



1 + 



Ett 



2{F 1 +M) 
E„ 



Ett 

E'+M 
1 



isi b ) 
& = q u - % 

1 + 



Ett 



2(F 2 +M) 

Ett 



F 2 +M 

Ett 

F'+M 
1 



Table E.2: Coefficients of the partial wave decomposition of the direct-production diagram. 
The b) superscript refers to the 1^2 diagram. 
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FSP) 


FSI b > 


ISP) 


isi b ) 


Qi 


qbf q-K 


qbf y-7r 


Haf \ Hit 


— n £ A- n 

qaf ~ q-K 


k 


qbf <4af 2 


<Jbf qaf i 2 


qbf q a f 2 


fLt — n t A- ^ L 
qbf qaf i 2 


A 


1 


1 


1 


1 


Fj+M 


F 2 +M 


M-E' 


M-F' 


B 


1 


1 


1 


1 


M-E 


M-F 


Ei+M 


E 2 +M 


C 




1 1 1 






1 1 1 




1 


1 




M-E 1 2(Fi+M) 






M-F 1 2(F 2 +M) 




2(M-F') 


2(M-F') 



Table E.3: Coefficients of the partial wave decomposition of the cr-exchange part of the Z- 
diagrams Eq. flEj}. Qbf = V and Qaf = Qu- 




Table E.4: Coefficients of the partial wave decomposition of the cu-exchange part of the 
Z-diagrams for the FSI case Eq. (IE. 1ft . 



A 
B 

C 



ISP) 



i 

M-E' 
1 

Fi+M 
1 



2M 



+ 



1 

F 2 +M 

E 2 +M ' 2M 
1 



2(M-E') 2(F 2 +M) 4M 



isi b > 



1 



1 



Fi+M M—F' 2M 
1 1 



Ei+M E 2 +M 

1 1 
2(M-F') 2(Ji+M) 



2M 



4M 



Table E.5: The same of Table lE~4l for the ISI case 



A 
B 

C 



FSP 



i 



F 2 +M 
1 

' E 2 +M 
1 



2M 



2M 



2(F 2 +M) ^ 4M 



FSI b ) 



l 



Fi+M 
1 

Ei+M 
1 



2M 



2M 



2(Fi+M) 



4M 



ISP» 



1 



F 2 +M 
1 



2M 



E 2 +M " r 2M 



1 



2(F 2 +M) 4M 



FSI b > 



i 



Fi+M 
1 



2M ) * 



Fi+M 
1 



2M 
f - 

2(Fi+M) r 4M 



Table E.6: Coefficients of the partial wave decomposition of the Lj-exchange part of the 
Z-diagrams Eq. (IE. 46ft . 
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Appendix F 

Numerical evaluation of integrals 
with pole singularities 

The integrals appearing in usual T-matrix equations are of the form (see Appendix ID|) : 

where a pole exists at k' = k. g (k') is a regular function of k' (i. e., it has no poles). The 
ie term of the denominator of Eq. (jF.ljl can be evaluated by using the Cauchy Principal 
value theorem |lU3j . 

r dk'k' \, 9 f') = PV r dk'k ,2 -^--m r dk'k ,2 6 (k') g (k 2 - k' 2 ) , (F.2) 

Jo k 2 - k' 2 + ie J k 2 - k' 2 J 



where 

r°° . n(h'\ r r k ~ e . n(V\ r°° . n (b>\ 1 

(F.3) 



PV / dk'k' 2 ^ > = lim 
A; 2 — /c' 2 e^o 



£.2 _ U/2 U2 _ U/2 

_J0 * * Jk+e * * 

The pole at k' — k in Eq. (|F.2jl has to be treated carefully. This can be done numerically 
by using the subtraction method with regularisation [TDK] to calculate the integral. The 
basis of the subtraction method is to treat the pole in the numerical integral by subtracting 
(and adding separately) an analytically defined integral 



oo 



9(h') _ mr r.rjfrn 9(k>) g{k) 



o 



-pv / dk'k 2 9 {k) 



-m 



n 



k 2 - k' 2 
dk'k ,2 5 (k 2 - k' 2 ) g (k') 
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Numerical evaluation of integrals with pole singularities 



The first integration on the right hand side of Eq. (|F.4|) can be performed numerically 
because the pole no longer exists (when k' — > k, it goes to jj). The second integration 
may be solved analytically as 



h = PV 
= PV 



lim 



g (k) 



dk'k 



'7„2 



k 2 - k' 2 



dk'k 2 g (k) 



1 



1 



1 



In 



k + k' 



k'-k 



k—e 



2k \ k - k' k + k 1 
k + k' 



+ 



In 



k — k' 



~| fcmax 

k+e 



;9(k) 



2k -e 


-In 


k 


+ In 


kmax k 


-In 


2k + e 








e 




k 


k—k 




e 



lim ( In 



-L(k)\n kmax k 
2 k max + k 



- 2 9(k) 



(F.5) 
(F.6) 
(F.7) 

(F.8) 
(F.9) 



In the limit k n 

reads 



dk'k 



oo the integral of Eq. (|F.9|) reduces to zero, and in that case, Eq. (jF.4 

g(k') 



nj2 



k 2 — k' 2 + is 



PV 



PV 



dk' k' 2 



g(k') 

k 2 - k' 2 



k 



dk' k' 2 



2 9(k) 
k 2 - k' 2 



g(k') _ k2 g{k) 



i^kg(k). 



(F.10) 

dk'k' 2 5(k 2 -k' 2 )g(k') (F. 11) 

(F.12) 



The integral of Eq. §FJM is no longer a numerical problem. To evaluate it numerically, 
one simply needs to map a variable x in a finite interval to k' |1()9[ IllOj . For instance, 
using 

fc' = ctan(|x) withO<x<l, (F.13) 

where c is an appropriate constant, or 

x+1 



k' 



1 — x 



with — 1 < x < 1. 



(F.14) 
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Appendix G 

Cross section for pion production 



The unpolarised cross section for pion production in nucleon-nucleon collisions can be 
calculated from |1 11] 

/ (2Mi)(2M 2 ) (27r)454 ( ^ + p 2 _ Qi _ Q2 _ Qn) (G1) 

J - '(Pi ■ P2) 2 - M 2 M 2 



(2Mj) (2M 2 ) d 3 qi d 3 q 2 d 3 q n s 
(2tt) 9 gi °g 2 2g0 

where 5 is the symmetry factor which accounts for different particles in the final state. 
Pi (Qi) ar e the nucleons initial(final) four-momenta. Note that a is built up from three 
factors, namely the flux factor, 



^ " ^ M> M," = MA ' (G ' 2) 



the phase-space, 



Mxfqi M 2 d i q 2 d 3 q T 4 4 

9i(2tt) $( 27r ) (2?r) 2g0 

and the amplitude for the process (averaged over the initial spins and summed in the final 
spins) , 

spins 

Choosing the overall centre-of-mass system, Eq. (|G.1|) reads, 
M 2 1 f M 2 1 2 
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Cross section for pion production 



where 2E = y/s is the initial CM energy, E k = V M 2 + k 2 is the energy of each final 



nucleon relatively to the CM of two final NN system. E pp = 2E — E n and E n = \Jq 2 + m 2 
are the energy of the final NN system and the energy of the produced pion relatively to 
the CM, respectively. 

The integral in d 3 q pp can be easily calculated using the 5 3 -function, 

da M 2 1 r M 2 1 
dn k dn n ~ ApE (2tt) 5 J ~E~~E~E k 

The remaining 5 function can then be used to perform the dk integration. From 



r m 2 1 2 

■ S / — — — \M fi \ 5 (2E - Epp - E v ) k 2 qldkdq n . (G.6) 

J tipp^k 



^ 2 ^s-2 v rsE 7T + E* = 4(M 2 + e)+ql, (G.7) 

one concludes that the argument of the 5-function vanishes for 

klro =\{ s ~ 2>/S>A£+£ + 4M 2 ) . (G.8) 

Since 



E pp 



2VM 2 + k 2 



2 



ql 2EppdE P p = Skdk, (G.9) 



the final expression for the cross section is 1 

da M 2 1 ATM 2 1 



dQkdQ^ ApE (2tt) 5 4 J E n E k 



i r m i 2 

S- / — \M fi \ k zero qldq n (G.10) 



K 2 \M fl \ dq wi (G.li; 



where 



M 2 1 1 [M 2 1 2 , 2 



K t = — I ST ~B fceroS- (G.12) 



4pE(2vr) 5 4 J E n E 
Eq. (|G.10|) can also be written in terms of an integral over the pion energy as 

da M 2 1 „ f E ™ M 



dn k dtt w ApE (2tt) 5 



M 2 2 , 

S / \M fi \ k zero ^/El - mldE v , (G.13) 



1 To express the cross section in barn, Eq. (|G.10|I must be multiplied by a factor of {he) 2 x 10 2 . 
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Appendix H 

Clebsch-Gordan coefficients, Six-J 
and Nine-J symbols 

This Appendix lists a series of useful relations of tensorial calculus [lU6t I1U7] employed in 
the partial wave decomposition of the amplitudes (see Appendix IE]). 



Legendre Polynomials 



«<*>=sd? (*•-!)' (h.i) 

P, (x) P v (x) dx = — ^—fc (H.2) 
00 91 + 1 

—^-Pi (*) = ^ (x - x') (H.3) 



-1 

d(a-b)f(a-b) Pi (a ■ S) (H.5) 
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Clebsch-Gordan coefficients, Six-J and Nine-J symbols 



Spherical Harmonics 

Y^(-a) = (-) l Yi(a) (H.6) 

Yi{*-0,* + <p) = {-?Yi>(6,<p) (H.7) 

Y^(a) = (-) m YL m (a) (H.8) 

Y l m (0) = (H.9) 

J dqY^(q)Yi,(q) = S U 'Smm' (H.10) 

(m 

(p) >i* (p) = (-i) m2 -4= v/^TTv/^TT (H.12) 



Ym (p) m ? ft^i - m 2 |LM) (Z^OOlLO) 



E ij ((9, ^ (6K, = 5 (cos 9 - cos 0') 5 (cp - V >) 



(H.13) 



tin 



E y ™ «) y ™ (? ) = \/ ( R) = ^P (g; s') 



47T 



47T 



(H.14) 



27 —I— 1 



m=—l 



4tt 



with a = £[(9,ip);(0',ip')} (H.15) 



Legendre polynomials related to spherical harmonics 

-1)' 



^(«-&)=-j=4^°(a,&) (H.17) 
3&(M) = \l {2k l l ^l\^ l) (hW\lO)Yi(a) (H.18) 
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Coupling of angular momenta 

Clebsch-Gordon coefficients 

Closure relations / — 

( _ 1) -2i 1 +2, 2 -m-w /^±1 ( jlhmim2 \ jm ) ( Jlj2mim2 |jW) = 6 jf 5 mm , (H.19) 



J2(-l)- 2n+2j2+2m (Jij2m im2 \j - m) {h32m' x m' 2 \j - m) = <W/ m2m > (H.20) 

jm 

Symmetries 

(jxj 2 m 1 m 2 \jm) = (-1) J1+J2 ~ J (j x j 2 - m x - m 2 \j - m) (H.21) 

(ji32rn 1 m 2 \jm) = (-1) J1+J2 ~ J (j 2 j 1 m 2 mi \jm) (H.22) 



(jij 2 m 1 m 2 \jm) = l/ g^ + \ ^ (jiJ m i ~ 171 \h ~ m i) (H.23) 



2 J + 1 m 

(j\j 2 m x m 2 \jm) = x \ — - - (_i) j2+m2 (jj 2 _ mm2 | ^ - mi ) (H.24) 



Special cases 



(Jihmim2 |00> = (-l) J1 - mi v^j7Tl5 JU2 5 mi ,_ m2 (H.25) 
(jiOmiO |jm) = 5 hj 5 mim (H.26) 



3j-symbols 



(-l) jW2+m j| Jl J2 J 
mi r«2 — m 



6j-symbols and Racah coefficients 

Symmetry relations 

jl 32 h\_l 32 33 Jl ^ _ J J3 Jl J2 J2 Jl J3 



J4 J5 J6 J5 J6 J4 J6 J4 J5 J5 J4 J6 



J2 


33 


jl 


35 


36 


J I 


jl 


35 


33 


jl 


j2 


36 
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Clebsch-Gordan coefficients, Six-J and Nine-J symbols 



Orthogonality relation 

{ 31 32 f ) I * 32 f \ - Sj & 



.1 



33 34 3 \ \ 33 34 3 3 



Relation to Racah coefficient 



( 33 J J23 J 
Relation to 3j symbols 



lymw w {jij2 j j3 . Ji2 j 23) (H . 30) 



W ( 3l j 2 Jj 3 ; J 12 J 23 ) = (2J+l)(-l) J23 ~ Jl2 ^ I 31 32 JU ] (H.31) 

33 J Jl2 \ [ 32 33 J23 \ I J jl ^23 

^7 e I \ j3 5 —(p I \ 7 a <p 

Sum over components 

WUAJJ*JM\ 3 31 ' ) = Y,(-lY*- h *- a -' ( 31 32 JU ) (H-32) 

\ 7 a <p J /36s y a (3 -e 

33 J Jl2 \ [ J2 J3 ^23 

5 7 e / I /3 S -(f) 



I 1 \J23~ Jl2~S 

W V&Jfa JuJn) (JjijalJ* - <f>) = Jj ZTi (H.33) 

p5e (2J 12 + 1) 

(jij2a/3\J 12 s) {jsJSil J12 - e) 

{3233PS\J 23 (f)) 



9j-symbols 

Relation to 6j-symbols 

3i 32 J 12 

33 34 J34 
Jl3 J24 J 



J\2 J34 J 

J' 3i 32 



jl 33 Jl3 
J24 J J' 



32 J4 J24 

33 J' J34 



X 



(H.34) 
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Relation, to 3i-sVmbols 

7i I2 Jw. " 



Jl 32 -J 12 
h fa J 34 
Jl3 J24 J 



> = (2ii + l)E 



jl 32 J12 

a (3 7 



32 fa -hi 
per] 



J12 J34 J 

7 (j) V 



jl 33 Jl3 

a 5 p 



33 J4 J34 

5 e d> 



F elatior^to julet >sch-Gordan coefficients 

(2jl + 1) ( — \^~ 2 h-Ji-2+J3A+2jA-Jl3 + J2A 



33 34 J34 > 
J 13 J24 J 



V2Ji 2 + V2J34 + V2J13 + V2J 2 4 + 1 (2 J + 1) 



(H.35) 



Jl3 


<^24 


J 


P 


V 


V 



(H.36) 



E(- 1 )" P "^ 2 " 7 "" 0'U2«/?|J 12 - 7 ) ^24 - rf) 

(J12J341MJ - v) (jij 3 aS\J 13 - p) {j 3 faSe\ J34 ~ <P) (Jl3J24PV\J - v) 



Sum over components . T 

\ T 7i .72 J12 

jl 32 J 12 



a (3 7 



Ji J2 J12 

J3 J4 J34 
J\3 J24 J 



> 



J2 J4 ^24 

Per] 



(H.37) 



jl J3 ^13 

^(-l)»-"+<^ (j 2j4 /3 £ | J 24 - (J l2 Jul<P\J- v) x 

(jlj3tt5|^13 - p) (j334&e\ J34 - 0) (J 13 J 2 4pV\ J ~ V) = 



Jl3 


J 24 


J 


P 


V 


V 



J12 J34 J 

7 1/ 

33 34 J34 

5 e 6 



(H.38) 

( -^^2,32+2^4+^12+7— J34+J13— J24 

yg J34 + lyg Ji 3 + 1 y/2 J 24 + 1 (2 J + 1) 
V2J12 + 1 

jl 32 J12 
(3i32aP\Ji2~l) { j 3 j 4 J34 

Jl3 J24 J 



Syfhmetry rel&tion£ 

1 Jl J2 J J3 



< J4 J5 J6 
J7 J8 J9 



Jl J4 J7 

32 fa fa 

33 fa fa 



37 3s 39 

jl J2 fa 
34 35 fa 



34 J5 J6 
37 js J9 
jl fa 33 



(H.39) 



jl +.72 +J3 + J4 +36 +36 +37 +3S +39 



34 J5 J6 
jl J2 J3 

37 js 39 
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Clebsch-Gordan coefficients, Six-J and Nine-J symbols 



Orthogonality 

„ ^_2— -2 i 
/ J Jl3 J2A < 

J13J24 



relation 



3 1 


J2 


■ J \2 




Jl 


J2 


J' 

J 12 


h 


Jl 




> < 


J3 


Jl 


J' 


■hi 


J24 


J 




Jxz 


J24 


J 



- — -2- — 2 
J\2 J34 



;H.40) 



Rules for tensorial algebra 



■x/3{a«®6 (1) }i 



a ■ b 



a-b = -bVArr la 1 &Y 1 [b 



with 

(o) 



-1,0,1 



(o) 



a x b 



(i) 



(i) 
1" 



(o) 




a- (b-c) = --LX:(-l//{[aW®6W] (/) ®c«}| 

6xc) =zv / 6{a (1) ® ® c (1) ] (1) }^ 
x tj ■ c = iVQ { [a« ® &«] (1) ® C W}^ 

(a • ft) (c • d) = 3 { [a« ® (0) ® [c« ® d«] (0) }^ 
x &) (c ■ d) = iV6 { [a« ® 6«] (1) ® [c« ® (0) } (1) 
x fe) ■ (c x d) = 2v/3 { [a« ® &«] W ® [c« ® eZ«] (1) } ^ 



(H.41) 
(H.42) 
(H.43) 
(H.44) 

(H.45) 

(H.46) 
(H.47) 

(H.48) 
(H.49) 
(H.50) 



Relations for tensorial operators 

{ A (a) g, B (6) = (-l) 2c3 ^Cg 3C3 A oa3 5 663 (H.51) 

13^3 

^3^3 = Cf 3 C b3C3 {A {a) ® 5 W }^ (H.52) 

cc 3 

{^°><8>B< 6 >}^ = (-l) a+6+c {B (fe) ® if [A, 5] = (H.53) 
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(d) i ( e ) 



/ | e a / 



(H.54) 



a) 



a )®£( 6 )] (/) <»c (c) 



(e) 



"I) 



a+6+c+e 



b a f 



d e c d 

(d) i ( e ) 



;h.55) 



a ) ® i?W] (c) ® [£>«> s = ^c/ji 



a 6 c 
d e / 
j k g 



(H.56) 



Spherical (tensor) harmonics 



mim2 

dadby 3 ™ 2 (a,b)y^f (a,bj = S jf 8 mm >8 jlj{ 8 hj > 2 



YL S+b 



^ ahbh 

h+h=l \a + b 



4tt (2/ + 1)! 



'V (2/ 1 + l)!(2/ 2 + l) 



Pi(d-b 



-I) 1 



47T 



V2Z + IP/ (a ■ b 



V21 + 1 



' 4tt^° ( 



(2Zi + 1) (2Z 2 + 1) 
4vr (2/ + 1) 



'nnn 1 » 



000 - 1 m 



(H.57) 
(H.58) 

(H.59) 

(H.60) 
(H.61) 

(H.62) 
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Clebsch-Gordan coefficients, Six-J and Nine-J symbols 



yf g (a + b,c + d 



a 



H? 2 c 9l d 92 / 4tt(2/+1)! 



/i+/ 2 =/ \a + b c + d 

91+92=3 I 



(2/ 1 + l)!(2/ 2 + l)! 



4tt (20 + 1)! 



(2^ + 1)! (2(fe + l)! 



! 1^/1/2 



(0 



(H.63) 



3^1 a + &,c 



h+h=f \a + b 
yf if (a,b)®Y*(c) 



4tt(2/ + 1)! 



/V (2/ 1 + l)!(2/ 2 + l)! 

(0 



y}™ (a,c + d 



91+92=9 



c 9l d 92 / 4tt(2c/ + 1)! 



c + eZ 



a) ® y 9 9192 ( a, d 



' (2^ + 1)! (20 2 + l)! 

(0 



{y L LlL2 (P, q) ® 3^ 2 (p, 9)}^ = E ^LlL^LzhC^C^X 



hh 



x < 



Li L2 £ 
Zi Z 2 / 

h h £ 



c 

J M 



X 



Ji 32 3 
k\ k 2 k 

h h £ 



hh 



> x 



.VI 



(H.64) 



(H.65) 



(H.66) 



(H.67) 



Wigner-Eckart theorem for the matrix elements of 
spherical tensors 



(J>M> \T KQ \ JM) = C^ M , (J 1 \\T K \\ J) (H.68) 
178 



179 



with 1 

{J' \\T K \\ jf cL m = (-1) 2K (-1) J '-V+V {J > || Tjc || jjHet EE] (R6g) 



^ince (J'M> \T KQ \ JM) Rct ^ = (-l) 2K CfMg, (J' \\T K \\ J) Ret mi . 



y MQW 
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